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ABSTRACT 

Advances  in  x-ray  imaging  technology  and  x-ray  sources  are  such  that  a  new 
technology  can  be  brought  to  commercialization  enabling  the  three-dimensional  (3- 
D)  microvisualization  of  hydrated  biological  specimens. 


Statement  A  per  telecon 
Dr.  Paul  Kepple  NRL/Code  l\121 
Washington,  DC  20375-5000 
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I.  INTRODUCTION 

The  Company  is  engaged  in  a  program  whose  main  goai  is  the  development  of 
a  new  technology  for  direct  three  dimensionai  (3-D)  x-ray  holographic  imaging.  It 
is  beiieved  that  this  technology  will  have  a  wide  range  of  important  applications 
in  the  defense,  medical,  and  scientific  sectors.  For  example,  in  the  medical  area, 
it  is  expected  that  biomedical  science  will  constitute  a  very  active  and  substantial 
market.  The  basis  of  this  view,  is  represented  directly  below. 

The  application  of  physical  technologies  for  the  direct  visuali2ation  of  biological 
entities  has  had  a  long  and  extremely  fruitful  history.  The  invention  of  the  light 
microscope  in  the  17th  century  and  the  development  of  the  electron  microscope  shortly 
before  World  War  II,  have  obviously  been  enormously  successful  scientifically. 
Equally  significantly,  these  two  landmark  advances,  in  addition  to  revealing  radically 
new  physical  features  of  the  human  environment,  have  also  had  a  profound  and 
unexpected  influence  on  man's  spiritual  perception  of  his  world.  The  light  microscope 
opened  up  an  unseen  universe,  not  only  of  strange  plant  and  animal  life,  but  also 
one  embodying  new  shapes  and  forms,  serving  to  challenge  and  stimulate  the 
mind.  The  electron  microscope,  by  greatly  enhancing  the  spatial  resolution  achievable, 
led  to  further  seminal  findings,  such  as  the  first  views  of  viral  particles  and  the 
complex  cytoskeletal  structure  of  cells.  As  a  consequence  of  the  broad  biological 
and  medical  applications  of  these  two  technologies,  both  now  represent  valuable 
basic  and  widely  used  tools  in  biological  research  and  clinical  medicine. 

Presently,  the  opportunity  exists  to  take  another  step  in  this  historical 
process  with  the  development  of  a  powerful  new  means  of  visualization,  three- 
dimensional  ^-ray  microholography.  ^  Importantly,  unlike  electron  microscopy, 
which  generally  requires  stained  and  desiccated  specimens,  this  new  technology 
will  uniquely  combine  a  main  feature  of  the  light  microscope,  the  ability  to  observe 
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hydrated  living  matter,  with  the  principal  capability  of  the  electron  microscope, 
the  power  to  resolve  very  small  spatial  elements,  such  as  those  defining  the 
cytoskeletal  architecture. 

This  microholographic  technology  is  composed  of  two  basic  and  separate 
components.  They  are  (A)  an  x-ray  laser  and  (B)  an  x~ray  holographic  camera 
designed  specifically  for  biological  specimens.  The  suitability  of  x-ray  laser 
technology  for  this  type  of  imaging  arises  very  naturally,  since  x-rays  (i)  can 
readily  penetrate  hydrated  matter,  (ii)  intrinsically  embody  a  short  spatial  scale- 
length  permitting  high  resolution,  a'i  (iii)  can  attain  extraordinarily  high  exposure 
rates  as  a  direct  consequence  of  the  very  high  source  brightness 

The  holographic  camera  has  been  sufficiently  well  developed  by  th'  Company 

\ 

that  a  prototype  instrument^"'^  is  currently  under  development.  Indeed,  since  the 

capability  of  this  technology  is  quite  general,  it  has  been  shown  by  the  Company^ 

that  (t  is  possible  to  develop  a  technique  of  direct  imaging  for  the  sequencing  and 
mapping  of  the  human  genome  based  on  a  modification  of  this  prototype.  However, 
as  noted  above,  the  full  potential  of  holography  fo.  hydrated  specimens,  can  only 
be  realized  when  it  incorporates  a  suirable  x-ray  laser,  The  resulting  3-D  images, 
on  account  of  the  rate  of  radiation  damage,  require  a  source  with  a  brightness 

that  can  only  be  achieved  with  an  x-ray  laser.  It  is  important  to  note  that  the 
use  of  such  a  laser  uniquely  provides  two  important  new  functions,  atomic  specificity 
and  "stop-action”  imaging.  By  appropriately  selecting  the  laser  wavelength,  it  is 
possible  to  selectively  enhance  the  imaging  of  a  particular  atomic  species  of 
biological  interest  such  as  nitrogen,  phosphorus,  sulfur,  calcium,  and  potassium, 
among  others.  Finally,  the  exposure  time,  predicted  to  be  ~  10'^’  seconds,  is 

faster  than  most  biological  or  chemical  processes. 

Special  conditions  are  required  to  produce  the  amplification  needed  for  iasei 
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action  at  x-ray  wavelengths.  Amplification  In  the  x-ray  region  mainly  requires 
prodigious  energy  deposition  rates  spatially  organized  in  a  high-aspect-ratio  volume 
of  material.  The  fundamental  requirement,  therefore,  reduces  to  the  efficient 
deposition  of  energy  at  very  high  specific  powers  in  a  spatially  controlled  and 
nontherrnal  manner.  Canonically,  a  deposition  rate  in  the  10^^  -  10‘®  W/cm^ 
range  is  necessary. 

To  meet  this  requirement,  we  have  developed  an  approach  that  uniquely 
combines  three  different  components;  (1)  an  energy  deposition  mechanism  based 
on  a  highly  nonlinear  coupling®'^  to  the  source  of  energy  used  for  excitation  of  the 
x-ray  amplifying  medium,  (2)  a  mode  of  channeled  propagation®"^®  capable  of  guiding 
both  the  ultraviolet  and  x-ray  radiation,  and  (3)  a  now  extremely  high  brightness 
ultraviolet  pulsed  laser  technology^ to  provide  the  source  of  energy  for  excitation. 

It  can  be  seen  that  these  three  elements,  which  appear  capable  of  producing 
and  controlling  the  most  extreme  conditions  of  power  density  yet  realized  v/ith 
laboratory-scale  apparatus,  fit  together  in  a  remarkably  congenial  way,  as  outlined 
in  Appendix  A,  Importantly,  the  conditions  needed  to  produce  the  (1)  strong 
nonlinear  coupling  observed  experimentally  are  exactly  (2)  those  required  for 
channeled  propagation  and  (3)  that  these  conditions  can  be  achieved  by  present 
laser  technology.  The  confluence  of  these  three  conditions  is  the  key  feature  of 
this  method  for  the  efficient  attainment  of  the  very  high-energy-density  states  of 
matter  needed  for  x-ray  amplification. 

This  report  deals  with  two  central  aspects  of  the  developing  technology 
which  cirectly  affect  the  feasibility  of  commercial  application.  They  are  (l)  the 
technology  for  holographic  image  reconstruction  and  (2)  the  x-ray  source  technology 
in  particular,  the  issue  associated  with  the  optimum  spatial  control  of  the  deposition 
of  energy.  The  effort  on  both  of  these  questions  is  aimed  principally  at  reducing 
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the  complexity,  cost,  and  physical  scale  of  the  instrumentation  needed  to  serve 
the  targeted  applications.  Most  significantly,  the  results  show  indeed  that  a  small 
laboratory  instrument  is  possible.  Full  satisfaction  of  the  Statement  of  Work 
represented  in  the  Phase  II  proposal  has  been  accomplished. 

II.  PHASE  II  RESEARCH 

The  Statement  of  Work  for  the  Phase  II  effort  is  presented  for  reference  in 
Appendix  B. 

A.  DEVELOPMENT  OF  ALGORITHMS  FOR  HOLOGRAPHIC  IMAGE  RECON¬ 
STRUCTION 

1.  3-D  Algorithms 

Algorithms  suitable  for  the  reconstruction  of  3-D  images  in  the  x-ray  range 
have  been  developed  and  tested  by  appropriate  simulations  at  visible  wavelengths. 
Appendix  C  gives  a  description  of  the  important  aspects  of  this  work. 

2.  Images  of  Drosophila  Embryo 

Other  subsequent  experiments  have  shown  that  a  properly  prepared  Drosophila 
embryo  can  serve  as  a  very  convenient  and  informative  model  system  for  examining 
the  3-D  imaging  properties  of  the  microscope.  Specifically,  since  we  desire  images 
with  (1)  high  spatial  resolution,  (2)  biological  structural  specificity,  (3)  3-D  informa¬ 
tion,  and  (4)  high  contrast,  a  staining  procedure  involving  silver-enhanced  immuno- 
gold  complexes^ ^  was  used.  This  stain  was  used  to  visualize  the  segmentation  of 
the  Drosophila  embryo  as  displayed  by  using  an  antibody  against  the  engrailed 
protein. 

A  reconstructed  image  of  a  Drosophila  embryo,  at  three  successive  depths, 
(z),  is  shown  in  Fig.  (1)  taken  under  dark-field  illumination.  The  study  of  such 
images  is  informing  on  the  operational  characteristics  of  the  system  such  as  (l) 
ti'.e  sensitivity  of  the  image  to  the  conditions  of  exposure,  (2)  the  staining  procedure, 
and  (3)  the  numerical  procedures  used  in  the  reconstruction.  These  characteristics 
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represent  simulations  of  their  counterparts  In  the  x-ray  regime.  Of  course,  the 
high  contrast  that  can  be  naturally  obtained  in  the  x-ray  region  will  eliminate  the 
need  for  staining. 

Several  important  aspects  of  the  imaging  system  are  revealed  by  examining 
the  reconstructions  shown  in  Fig.  (i).  They  are  (1)  the  high  contrast  shown  by 
the  banding,  (2)  the  sharp  edges,  demonstrating  the  high  spatial  resolution,  shown 
most  clearly  at  the  edge  of  the  embryo.  (3)  the  absence  of  speckle,  and  (4)  the 
3-D  character  of  the  images.  The  latter  point  is  apparent  through  examination  of 
the  region  of  invagination  toward  the  bottom  right  of  the  embryo. 


Fig.  (1);  Reconstructed  images  of  Drosophila  embryo  stained  with  silver-enhanced 
Immuno-gold  complexes.  The  image  reveals  the  banded  pattern  governed 
by  the  seg.mentation  genes.  The  three  panels  represent  reconstructions 
corresponding  to  three  different  depth  (z)  planes  in  the  specimen.  In 

the  left  panel,  z  =  i-54  pm;  in  the  center  panel,  z  =  0  pm;  in  the  right 

panel  z  =  -54  pm.  The  full  size  of  the  reconstructed  image  is  -  256 

pm,  as  indicated  by  the  black  edges  at  the  left  and  right  of  the  embryo. 
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B.  HOLOGRAPHIC  MICROSCOPE 


Two  basic  types  of  holographic  microscopes  have  been  constructed  using  the 
Fourier  configuration.  One  version  utiiizes  a  giyceroi  microdrop  as  the  reference 
wave  scatterer.  This  instrument  is  described  in  Appendix  D. 

A  second  basic  system  we  have  used  for  the  studies  conducted  with  visible 
radiation  is  shown  in  Fig.  (2).  A  convenient  configuration  is  a  modification  of  a 
Mach-Zehnder  (MZ)  interferometer  geomefy.  In  this  configuration  have  we 
replaced  the  spherical  reference  scatterer  with  a  suitable  microscope  objective, since, 
for  visible  wavelengths,  this  permits  facile  and  flexible  operation.  The  primary 

beam  of  coherent  radiation  is  split  into  two  beams  with  a  beam  splitter.  Advantage 
is  taken  of  the  inherent  polarization  of  the  coherent  laser  radiation,  and  a  polarizing 
beam  splitter  is  used.  For  these  studies,  the  polarization  is  utilized  to  control 

the  relative  intensity  of  the  two  subbeams.  In  the  x-ray  case,  this  relative 

intensity  would  be  governed  by  the  radius  of  the  spherical  reference  scatterer 

One  subbeam  of  the  M2  interferometric  section  of  the  microscope  is  used  to 
illuminate  the  biological  specimen,  while  the  other  subbeam  is  brought  to  a  focus 
by  illumination  of  a  compound  microscope  (CM)  objective.  Both  subbeams  are 
then  recombined  in  a  beam-combiner  cube  and  these  recombined  beams,  which 
form  the  holographic  interference  of  the  specimen,  diverge  to  fill  the  aperture 
subtended  by  the  detector  (CCD).  The  fourth  direction  of  the  beam  combiner  is 
occupied  by  a  CM  to  permit  alignment  of  the  specimen  and  the  reference  beam 
This  configuration  of  the  microscope  functions  as  a  dark-field  microscope  and  has 
performed  exceptionally  well  at  NA  ~  0.50  with  highly  scattering  biological  specimens 
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Fig.  (2):  Schematic  of  Mach-Zehnder  holographic  microscope. 


A  photograph  of  the  existing  instrument,  which  corresponds  to  the  schematic 
presented  in  Fig.  (2),  which  was  used  for  the  studies  at  visible  wavelengths  is 
given  in  Fig.  (3a).  The  specific  design  incoiporating  the  same  concepts  that  has 
already  been  built  for  work  at  x-rav  wavelengths  is  shown  in  Fig.  (3b).  As 
described  in  Appendix  E,  this  x-ray  camera  is  sufficiently  flexible  that  it  can  be 
used  with  various  sources,  such  as  a  free-electron  laser. 

Studies  of  the  characteristics  of  operation  of  this  instrument  have  been 
performed  on  several  biological  materials  These  include  (1)  parascaris  univalens  larva 
(2)  a  fish  melanoma  histopathology  specimen  (Xiphophorus  maculalus  x  helleri),  (3) 
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a  section  through  duodenum  including  a  nerve  gangiion.  and,  as  noted  above.  (4)  a 
Crosophiia  embryo.  The  latter  was  used  to  examine  specifically  the  three-dimensional 
character  of  the  imaging. 


Fig.  (3);  (a)  Photograph  of  the  visible  holographic  microscope.  (b)  Photograph 

of  x-ray  holographic  system. 


C.  ELECTROMAGNETIC  PROPAGATION 


As  outlined  in  Section  I,  the  spatial  control  of  the  energy  deposited  for 
excitation  of  an  x-ray  amplifier  plays  an  important  role  in  the  fundamental  scaling 
relationship  between  the  required  energy,  the  gain  and  the  wavelength.  The 
critical  governing  issue,  which  determines  the  scaling  relationship  between  the 
required  excitation  energy  (E)  and  the  amplifier  gain  (G)  of  x-ray  lasers,  is  the 
spatial  control  of  the  deposited  energy.  The  information  presented  in  Fig.  (4) 
shows  that  optimizing  the  gain  (G)  per  unit  energy  (E)  calls  for  the  guided  mode 
of  propagation  in  order  to  optimally  control  the  deposition  of  the  energy.  Overall, 
in  comparison  to  traditional  forms  of  excitation,  for  a  fixed  x-ray  energy  output 
(Ey)  and  wavelength  (Xy),  a  reduction  of  several  orders  of  magnitude  in  the 
necessary  energy  (E)  results,  as  shown  in  Fig.  (5),  if  this  form  of  confined  (channeled) 
propagation  can  be  achieved.  Therefore,  if  this  scaling  holds,  a  relatively  small 
and  useful  laboratory-scale  technology  becomes  feasible.  One  joule  of  energy  is 
seen  as  sufficient  up  to  a  wavelength  uf  several  kilovolts  (~  I.e  A). 

At  least  two  types  of  conditions  can  lead  to  the  desired  confined  propagatio  .. 
They  involve  (1)  the  use  of  guiding  strucfires,  such  as  static  waveguides,  or  (2) 
propagation  that  is  dynamically  confined  Roth  aspects  have  been  con.sidered. 
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Fig.  (4):  Spatial  distribution  of  energy  of  excitation  (E)  for  an  x-ray  (hoj^) 
amplifier.  Parameters  are  the  same  as  in  Fig.  (5)  with  X  the  wavelength 
of  the  excitation  energy,  assumed  longitudinally  delivered,  and  with  p, 
p*,  and  Pg  representing  the  particle,  inversion  and  electron  densities, 
respectively.  The  analysis  shows  that  optimization  of  G/E  requires  a 
guided  mode  of  propagation  so  that  high  concentrations  of  power  can 
be  organized  into  high-aspect-ratio  spatial  volumes. 
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Energy  (E) 


Energy  (EVWavelength  (X)  Scaling 


Fig.  (5):  Scaling  relationship  between  required  excitation  energy  (E)  and  quantum 
energy  (fiojx)  characteristic  of  the  amplifier.  Parameters.  total  gam 
exponent  G  »  100,  energy  efficiency  rij^  =  10'*,  channel  diameter  6  =  3 
pm,  x-ray  (hto^)  cross  section  for  stimulated  emission  0,^,  x-ray  cross 
section  for  stimulated  emission  for  radiatively-broadened  transition  a^p. 
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1.  Capillary  Guiding  Structures 

A  static  guiding  structure,  assuming  that  the  wave  can  be  properly  launched 
and  that  the  losses  are  sufficiently  low,  could  potentially  serve  to  provide  the 
desired  confined  propagation.  This  possibility  was  experimentally  explored,  under 
subcontract  to  the  Company,  by  the  University  of  Illinois  at  Chicago.  Basically, 
it  was  found  that  the  use  of  static  guiding  structures  is  not  feasible  for  the 
range  of  parameters  required  for  x-ray  amplification.  The  report  of  the  subcontractor 
is  given  in  Appendix  F. 

2.  Dynamical  Guiding  of  Electromagnetic  Propagation 

Theoretical  indicates  that  dynamical  guiding  of  very  intense  pulsed 

radiation  may  be  possible  under  conditions  suitable  for  x-ray  amplification.  The 
essential  aspects  of  this  theoretical  work  are  given  in  Appendix  A  and  Appendices 
(G  -  I). 

Significantly,  experimental  evidence  now  exists,  as  described  in  Appendix  J, 
that  demonstrates  for  the  formation  of  such  a  mode  of  propagation.  Specifically, 
the  experimental  studies  have  examined  a  new  relativistic  regime  of  high-intensity 
short-pulse  propagation  in  plasmas  which  indicates  the  development  of  a  stable 
mode  of  spatially  confined  (channeled)  propagation.  For  an  electron  density  of  - 
1.35  X  10*‘  cm’’  and  a  power  of  ~  3  x  10**  W,  the  results  indicated  a  channel 
radius  <  1  pm  and  a  peak  intensity  ~  10**  W/cm’.  Comparison  of  these  findings 
with  a  dynamical  theory  yield  agreement  for  both  the  longitudinal  structure  and 
the  radial  extent  of  the  propagation  observed. 

The  most  important  outcome  of  these  findings  is  that  it  enables  the  scaling 
relationship  shown  in  Fig.  (5).  The  existence  of  this  mode  of  propagation  makes 
possible  a  small  laboratory-scale  x-ray  laser  technology  that  can  serve  a  multitude 
of  applications  including  the  importani  biomedical  ones  of  interest  to  the  Company. 
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III.  CONCLUSIONS 


The  implications  of  developments  achieved  under  the  support  of  this  work 
for  general  applications  to  x-ray  imaging  and  the  microcharacterization  of  condensed 
matter  are  extremely  important  and  propitious.  In  terms  of  the  x-ray  source, 
they  are  (1)  that  a  properly  controlled  energy  deposition  rate,  sufficient  for  the 
production  of  stimulated  x-ray  sources  up  to  a  few  kilovolts  in  quantum  energy, 
can  now  be  achieved  with  an  excitation  energy  of  ~  l  J,  (2)  that  an  x-ray  output 
energy  of  ~  1  mJ  per  pulse  is  achievable  with  laboratorv-scale  technology,  and  (3) 
that  an  x-ray  beam  diameter  (-2-3  pm)  arises  as  a  natural  consequence  of  the 
physics.  These  parameters  represent  an  exceptionally  high  peak  brightness  figure 
that  permits  a  new  and  completely  unexplored  range  of  physical  measurements  to 
be  made.  Indeed,  a  high-brightness  source  of  this  nature  is  ideal  for  the  microimaging 
of  condensed  matter.  In  particular,  an  x-ray  source  with  these  parameters  is 

perfectly  matched  to  the  requirements  for  holographic  imaging  of  biological  materials^" 
in  terms  of  all  its  relevant  properties,  specifically,  wavelength  (10  -  40  K),  pulse 
energy  (-  1  mJ),  pulse  length  (~  10'**  s),  beam  diameter  (-2-3  pm),  and  divergence 
(~  1  mrad).  Finally,  the  3-D  camera  technology  matched  to  this  source  is  available 
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Abstract 

The  concentration  of  energy  arising  from  charge- 
displacement  self-channeling  Is  discussed.  Since 
high  energy  deposition  rates  are  expected  to 
arise  from  multiphoton  coupling  and  the  channel 
can  also  serve  as  an  effective  waveguide  for 
secondary  radiation,  such  circumstances  are 
ideal  for  generating  coherent  short  wavelength 
radiation. 

Channeled  Propagation 

Ampllficallan  In  the  x-ray  region  requires  prodigious 
energy  deposition  rates  (1]  spatially  organlz^  in 
a  hIgh-aspect-ratIo  volume  of  material.  The 
fundamental  question,  therefore,  centers  on  the 
controlled  deposition  of  energy  at  high  specific 
powers.  We  show  that  the  use  of  extremely 
Intense  (10‘*  -  10”  W/cm*)  short  pulse  (-  100 
fs)  radiation  may  be  able  to  produce  both  the 
necessary  deposition  rates  and  spatial  control  by 
combining  (1)  the  energy  deposition  [2,3]  arising 
from  high-ordar  multiphoton  processes  with  (2)  a 
new  mode  of  channeled  propagation  Involving  a 
charge-displacement  mechanism.  A  significant 
point  that  will  emerge  Is  that  the  conditions 
needed  lor  the  strong  multiphoton  coupling  are 
Identical  to  those  found  required  for  the  confined 
propagation. 

It  appears  [4-d|  that  a  fundamantally  rww  regime 
of  electromagnetic  propagation  will  develop  In 
plasmas  for  subpicosecond  radiation  of  sufficient 
IntensHy.  In  the  high  Intensity  case  of  Interest. 
In  which  processes  generating  Ionization  dominate 
the  coupling,  multiphoton  Ionization  Is  expected 
to  produce  a  substantial  reduction  of  both  the 
linear  and  nonlinear  refractive  Indicea  of  the 
medium.  For  a  suitably  Ionized  plasma,  this 
combination  of  high  charge  state  Ions  with  the 
electron  density  produced  by  the  multiphoton 
coupling  appears  capable  of  producing  a  new 
mode  of  radiative  channeling. 


Indeed,  preliminary  analysis  [7]  supports  the 
presence  of  the  following  nwchanlsm.  For  a 
sufficiently  short  pulse  (-  100  fs),  the  massive 
Ions  remain  spatially  fixed  while  the  relatively 
mobile  electrons  are  expelled  by  the  ponderomotlve 
force  from  the  high  Intensity  zone.  Thereby,  a 
state  of  equilibrium  can  be  established  between 
the  ponderomotlve  and  the  electrostatic  force 
densities  owing  to  the  charge  displacement. 
Since  the  electrons,  which  embody  a  negative 
contribution  to  the  Index,  are  expelled,  an  on- 
axls  region  of  relatively  high  refractive  Index  Is 
formed  which  supports  the  channeling.  In  a 
limiting  case.  It  appears  that  the  behavior  approaches 
that  characteristic  of  a  metallic  waveguide. 
Interestingly,  a  focusing  mechanism  of  this  type 
Is  not  dcscribable  In  terms  of  the  conventional  (8) 
nonlinear  Index  parameter  n,.  A  further  unusual 
aspect  of  this  process,  since  there  Is  a  strong 
tendency  to  locally  reduce  the  plasma  frequency 
dtp  In  the  region  where  the  Intensity  Is  high.  Is 
that  it  may  enable  propagation  for  appreciable 
distances  In  plasmas  that  would  normally  be 
considered  as  overdense. 

The  essence  of  the  Idea  can  be  understood  In 
the  foltowing  manner.  The  steady-state  force 
balance  between  the  radially  outward  ponderomotlve 
force  and  the  oppositely  directed  electron-ion 
attraction  Is  represented.  In  cylindrical  symmetry 
for  a  completely  Ionized  tenuous  (u  >  >  Up) 
plasma,  by 


2X0* 

m<i)*c 


71  (r)  +  e*n. 


1 


p(r)  (r  -  r‘) 
Tr-rf 


d*r' 


0.  (1) 


For  a  cylindrical  gausslan  Intensity  distribution 
l(r),  the  total  charge  density  p(r)  can  be  written  • 
In  exact  form  as 

p(r)  -  2B  {1-r*/r5  «tp(-r*/rS  (2) 

In  which  B  »  I,(ma)*cnprJ)'*  and  m,  c,  n^,  and 
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r,  d«no(*  tfM  •iMtron  mau.  tpsAd  of  light, 
quiascant  plasma  dantity,  and  radlua  (HWHM)  of 
tha  aaaumad  gautalan  IntonaKy  proflla,  raapaetivaly. 
Tha  form  of  this  solution  Is  shown  In  Fig.  (1)  for 
a  valua  of  B  ■  0.4.  Tha  axpulslon  of  tha  alactrons 
from  tha  cantral  raglon  Is  apparent. 

A  charge  distribution  of  this  natura  can  laad 
to  channeled  propagation.  An  astimate  of  the 
condition  necessary  ::an  now  be  made  by  describing 
the  charge  displacement  by  two  regions,  such 
that  p(r  <  r«ii^-p(0)  and  p(r  >  r^t/S) -p(r«wS), 
and  equating  the  an^e  of  total  internal  reflection 
to  the  angle  corresponding  to  the  first  minimum 
of  diffraction.  This  procedure  yields  a  critical 
Intensity 


64(1  +  e'*)rg 


nHii  u 

=  (4.3  X  10  *)  - ,  (3) 


In  which  r^  denotes  tha  classical  electron  radius. 
Significantly,  this  Intensity  Is  Independent  uf  the 
electron  density  n,  and  contrasts  with  the  power 
threshold  (6]  normally  arising  from  induced  Index 
charges  In  transparent  dielectrics.  Analytic 
approximations  extending  to  high  electron  density 
ah^  a  significant  dependence  of  on  n^  only 
quite  near  tha  critic^  density.  Furthermore,  It 
is  easily  demonstrated  that  Iq  normalized  to  the 
Compton  Intensity  (9)  Is  a  constant  (0.S4)  Independent 
of  fraquertcy,  a  flrxflng  which  shows  that  the 
charge-displacement  mechanism  Is  associated  with 
relativistic  conditions.  For  a  wavelength  of  248 
nm,  Iq  -  2.4  X  to‘*  W/cm*.  Finally,  a  striking 
aspect  of  the  analysis  (7]  Is  the  extremely  strong 
frequency  dependence  favoring  the  use  of  ultraviolet 
wavelengths  to  establish  the  conditions  for 
channeling.  Scaling  laws  are  derived  (7]  Indicating 
variations  as  rapid  as  u*. 


Charge -DIaplacement  Focuaing 

fl.ow  Oensity  Plasma) 


Figure  1.  Illustration  of  the  total  charge  density 
p  as  a  function  of  normalized  radial  distance  r/r, 
for  B  -  0.4. 


Controlled  Energy  Conceniratlon  -  Summary 

In  tarms  of  the  bealc  question  of  the  controlled 
deposition  of  energy  at  high  specific  powers,  the 
answer  that  emergea  from  these  considerations 
involves  three  separate  components.  They  are 

(1)  a  naw  ultravlolst  pulsad  p^sr  technology,  (10] 

(2)  an  energy  deposition  mechanism  bas^  on 
highly  nonlinear  coupling,  [it]  and  (3)  a  condition 
for  channeled  propagation.  [7]  It  can  be  seen  that 
thasa  threa  atamenta,  which  appear  capable  of 
producing  conditions  comparable  to,  or  possibly 
excaedlng  those  of  a  thermonuclear  environment, 
(10]  fit  together  In  a  ramarkabiy  congenial  way. 
Figure  (2)  illuatretee  theta  relationships.  A 
principal  finding  la  that  tha  radlativa  conditions 
naadad  for  tha  atrong  multiphoton  coupling 
governing  the  energy  transfer  rate  (it]  are  essentially 
Identical  to  those  required  for  the  channeled 
propagation  (7]  gag  that  the  later  technology  (10) 
can  readily  produce  the  regime  of  Irradiation 
necessary.  The  compatibility  of  these  three 
factors  Is  a  key  feature  of  this  method  for  the 
attainment  of  very  high  energy  density  states  of 
matter. 


Multiphoton  Energy  Tranafar  Croaa  Section 
Nty  A  - O  everything  (ionization) 


Figure  2.  Illustration  of  how  (1)  tha  energy  transfer 
cross  section  a,  (2)  the  channeling  condition,  and 
(3)  the  pulse  power  technology  all  fit  together. 
The  significant  fact  Is  that  the  conditions  needed 
for  the  strong  multiphoton  coupling  governing  the 
energy  transfer  (I  s  1^^  at  tO**  W/cm*)  are 
Identical  to  those  required  for  tha  channeled 
propagation  (If  =  I^^  *  10**  W/cm*),  and  that  the 
developing  large  aperture  ultraviolet  laser  technology 
moves  us  Into  the  desired  range. 
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Appendix  B; 


Statement  of  Work 


NRL 

STATEMENT  OF  WORK 
FOR 

MCR  TECHNOLOGY  CORPORATION 
PHASE  II  SBIR 


1.  Background 

Advances  in  x-ray  imaging  technology  and  x-ray  sources  are  such  that 
a  new  technoiogy  enabling  the  microvisualization  of  condensed  materials 
can  be  brought  to  commercialization.  In  order  to  accomplish  this  goal, 
further  development  of  the  computational  technology  of  holographic 
image  reconstruction  and  corresponding  aspects  of  source  technology, 
particularly  concerning  the  spatial  control  of  the  energy  deposition,  are 
needed.  This  Phase  II  SBIR  is  directed  at  these  objectives. 

2.  Scope 

This  is  an  SBIR  Phase  II  contract,  therefore  it  is  anticipated  that  a 
commercially  viable  product  will  be  developed  under  this  contract. 

3.  Applicable  Documents 

None. 

4.  Requirements 

Development  of  a  high  brightness  x-ray  source  for  directed  energy  and 

holographic  imaging  applications. 

4.1. 1A  The  contractor  shall  develop  advanced  algorithms  for  multidimensional 

holographic  reconstruction  from  a  holographic  microscope  operating  in 
the  visible,  but  able  to  simulate  the  parameters  of  an  x-ray  holographic 

microscope.  This  task  to  be  performed  in  conjunction  with  a  subcontract 
to  Robotronix,  Inc. 

4.1. IB  The  contractor  shall  explore  techniques  for  low  resolution  optical 

reconstruction  of  selected  portions  of  the  hologram  as  a  means  for 
rapid  screening  of  data  prior  to  numerical  reconstruction. 

4.1. 1C  The  contractor  shall  begin  the  technological  developments  that  will 

enable  the  performance  of  experimental  studies  to  examine  channeled 
propagation  in  capillaries  with  Nj  serving  as  the  medium.  In  these 
studies,  the  contractor  shall  place  emphasis  on  (1)  the  launching  of  the 
wave  in  the  capillary  and  (2)  determining  the  relationship  governing  the 
length  of  propagation  that  can  be  achieved  and  the  den.sity  of  the 
medium.  This  task  to  be  performed  on  a  subcontract  basis  with  the 
University  of  Illinois  at  Chicago. 

4.1, ID  In  conjunction  with  a  subcontract  to  Robotronix,  the  contractor  shall 

develop  supporting  theoretical  capability  for  the  analysis  of  the  experiments 
being  performed  for  in  task  4i  iC.  The  contractor  shall  place  particular 
emphasis  on  (1)  the  dynamics  of  channel  formation  and  (2)  channel  stability, 

4.1.1  Deliverable  products 
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Preprints  of  any  articles  submitted  to  Journals  for  Publication,  which 
cover  any  of  the  work  under  this  contract.  A  letter  type  or  formal 
progress  report. 

4.1.2  Schedule 

Tasks  4.1.1A  -  4.1.1D  are  to  be  completed  twelve  months  after  start  of 
work. 


4.2.1A  The  contractor  shall  design,  build,  test  and  operate  a  holographic 

microscope  operating  in  the  visible  region  suitable  for  simulalion  of  the 
parameters  governing  the  x-ray  holographic  application. 

4.2.1B  The  contractor  shall  continue  development  of  the  algorithms  of  task  4.1. 1A, 

4.2. 1C  The  contractor  shall  evaluate  techniques  for  microfabrication  of  target 

systems  for  both  visible  wavelength  simulation  and  x-ray  exposures. 

4.2.1  D  The  contractor  shall  continue  the  development  of  the  low  resolution 

optical  reconstruction  techniques  started  in  task  4.1. IB. 

4.2.1  E  The  contractor  shall  continue  the  experiments  noted  in  task  4.1. 1C. 

4.2. IF  The  contractor  shall  experimentally  explore  the  possibility  of  self 

channeling  in  Nj.  The  contractor  shall  evaluate  the  suitability  of  other 
materials,  both  atomic  and  molecular,  for  the  generation  of  considerably 
shorter  wavelengths.  This  task  to  be  performed  on  a  subcontract  basis 
with  the  University  of  Illinois  at  Chicago. 

4.2.1G  The  contractor  shall  continue  the  theoretical  developments  of  4.1, ID. 

4.2. 1H  The  contractor  shall  develop  a  prototype  of  the  target  conversion 

system  based  on  the  work  of  tasks  4.1. 1C,  4.1. ID,  4. 2. IE,  and  4,2.iF. 

4.2.1  Deliverable  products 

Preprints  of  any  articles  submitted  to  Journals  for  Publication,  which 
cover  any  of  the  work  under  this  contract.  A  final  report  containing 

designs,  calculations,  source  code  of  and  users  guide  for  any  software 
developed  wholly  under  this  contract,  drawings,  cost  estimates,  test 
plans  and  conclusions  resulting  wholly  from  tasks  4.1,1  A  -  4.1. 1D  and 

4.2.1A  -  4.2. 1H. 

4.2.2  Schedule 

Task  4.2.1A  -  4.2.1H  to  be  completed  twenty-four  months  after  start  of 
work. 

5.  Progress  reports 

The  contractor  shall  submit  a  letter  type  or  formal  progress  report  at 
the  end  of  iv/elve  months  and  a  final  report  by  the  end  of  twenty-four 
months,  in  addition,  as  directed  by  the  COTR,  the  contractor  shail 
submit  vicwgraphs  and  other  material  for  use  in  program  reviews. 
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Contract  No.  N00014~89-C-2274 


STATEMENT  OF  WORK  REVISIONS 

Tasks  4.2. IB,  4. 2. ID,  4.2. 1E  and  4.2. 1G  are  revised  to  read  as  follows: 

4.2.1  B  The  contractor  shall  develop  the  algorithms  r.r  task  4.1.1  A. 

4.2. ID  The  contractor  shall  develop  the  low  resolution  optical  reconstruction 
techniques  started  in  task  4.1, IB. 

4. 2. IE  The  contractor  shall  do  the  experiments  noted  in  task  4,1. 1C. 

4.2.1G  The  contractor  shall  do  the  theoretical  developments  of  4.1.1D. 
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FOURIER  TRANSFORM  HOLOGRAPHIC  MICROSCOPE 
Waleed  S.  Haddad,  David  Cullen,  Johndale  C.  Solem,  James  W.  Longworth, 
Armon  McPherson,  Keith  Boyer  and  Charles  K.  Rhodes 

MCR  Technology  Corporation 
P.  0.  Box  10084,  Chicago.  IL  60610-0084 

ABSTRACT 

We  describe  a  holographic  microscope  with  spatial  resolution  approaching  the 
diffraction  limit.  T.he  instrument  uses  a  tiny  drop  of  glycerol  as  a  lens  to  create 
the  spherically  diverging  reference  illumination  necessary  for  Fourier-transform 
holography.  Measurement  of  the  point-spread  function,  obtained  by  imaging  a 
knife  edge  in  dark-field  illumination,  indicates  a  transverse  resolution  of  1.4 
with  wavelength  X  =  514.5  nm.  Longitudinal  resolution  was  obtained  from  the 
holograms  by  the  numerical  equivalent  of  optical  sectioning.  We  describe  the 
method  of  reconstruction  and  demonstrate  the  microscope's  capability  v/ith  selected 
biological  specimens.  The  instrument  offers  two  unique  capabilities:  (1)  it  can 
collect  three-dimensional  information  in  a  single  pulse  of  light,  avoiding  specimen 
damage  and  bleaching;  and  (2)  it  can  record  three-dimensional  motion  pictures 
from  a  series  of  lig.'it  pulses.  The  conceptual  design  is  applicable  to  a  broan 

range  of  wavelengths  and  vvj  di.scuss  extension  to  the  x-ray  regime.  • 

Key  Words;  Microholography,  f  ourirr-Trarisform  Microholography,  Biological 

Samples,  Glycerol  Microarop,  ftecclution,  Numerical  Reconstruction 
Reference  Scaitt^rer 
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I.  INTRODUCTION 


Although  Gabor  invented  holography^  with  microscopy  applications  in  mind^, 
the  problem  of  recording  medium  resolution^  was  not  solved  until  the  introduction 
of  the  Fourier  transform^  technique.  Shortly  after  the  technique  was  proposed, 
Leith  and  Upatnieks'^  attempted  Fourier-transform  microholography  and  were  able 
to  resolve  a  bar  pattern  spaced  at  7  -  10  pm  using  a  632.8  nm  laser.  Although 
their  report  is  sketchy.  VanLigten  and  Osterberg^  claimed  an  optimum  resolution 
of  detectability  of  12  pm  and  very  serious  difficulties  in  obtaining  a  resolution  of 
~  5  pm.  Their  holographic  microscope  formed  a  hologram  from  a  pre-magnified 
image  rather  than  using  the  Fourier  method.  Apparently  because  initial  results 
were  disappointing,  visible-light  Fourier-transform  microholography  was  generally 
abandoned.  Most  papers®  of  the  70's  and  80's  are  concerned  with  its  application 
to  x-ray  microscopy. 

In  the  present  work,  we  describe  an  instrument  that  is  capable  of  resolution 
approaching  the  diffraction  limit.  The  camera  design  discussed  herein  represents 
a  new  technology  for  three-dimensional  visualization  that  is  applicable  to  a  spectral 
range  spanning  from  visible  (~  700  nm)  to  x-ray  (~  0.6  nm)  wavelengths.  Descriptions 
of  the  Fourier-transform  technique,  biological  target  design  and  related  numerical 

reconstruction  simulations  have  been  given  in  previous  work.^  The  new  camera 
system  has  been  operated  with  visible  radiation  (X  =  514.5  nm)  from  an  argon-ion 

laser.  The  experiments  in  the  visible  region  have  enabled  the  imaging  properties 

of  the  camera  to  be  studied  and  optimized  as  well  as  demonstrating  that  this 

camera  design  can  be  employed  for  applications  in  embryology  where  three-dimensional 
information  is  of  particular  value. 

^  The  original  concept  was  for  lensless  electron  microscopy. 
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II.  MICROHOLOGRAPHY  WITH  VISIBLE  RADIATION 


A.  Camera  Design 

Figure  (1)  is  a  schematic  of  the  microhoiographic  camera.  The  target  holder 
is  positioned  horizontally  so  gravity  will  hold  the  sample  and  a  glycerol  droplet  in 
place.  A  side  view  schematic  of  the  target  design  used  for  visible  light  exposures 
is  shown  in  Fig.  (2).  The  glycerol  droplet  proved  to  be  a  superb  lens  for  generating 
the  reference  illumination.  It  can  spread  the  illumination  to  the  large  angles  necessary 
for  good  resolution  and  can  be  placed  very  close  to  the  specimen.  Furthermore, 
its  size  can  be  varied  to  match  the  reference  intensity  to  the  intensity  of  radiation 
scattered  from  the  specimen. 

The  basic  functions  of  the  camera,  such  as  shutter  timing  and  triggering,  are 
controlled  by  a  personal  computer  (PC).  The  holograms  are  recorded  electronically, 
oigitized  and  stored  by  the  PC.  The  holographic  data  are  then  transferred  by  an 
ethernet  link  to  a  Stardent  GS  2000  supergraphic  work  station  where  the  numerical 
image  reconstruction  is  performed.  The  reconstructed  image  can  then  be  manipulated, 
scaled,  and  checked  for  three-dimensional  information  and  displayed  on  a  high 
resolution  graphics  monitor. 

The  holograms  were  registered  on  a  customized  charge-coupled  device  (CCD) 
having  a  sensitive  area  of  18.43  mm  x  18.43  mm  containing  2048  x  2048  square 
pixels  each  with  a  dimension  of  9  pm.  To  protect  the  CCD  from  the  direct  laser 
beam,  we  installed  a  cover  glass  with  a  central  opaque  square. 

B.  Liquid  Microdrop  Lenses 

The  drop  of  glycerol  is  placed  on  the  microscope  cover  slide,  close  to  the 
edge  of  the  sample  and  nearest  the  region  of  interest  in  the  sample,  as  shown  in 
Fig,  (2).  The  droplet  acts  as  a  lens  mtorcopting  a  portion  of  the  illuminating 
plane  wave  front  and  converting  it  into  a  diverging,  nearly  spherical  wave.  The 
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focal  point  of  this  lens  acts  as  the  point  source  of  reference  illumination.  Glycerol 
is  a  particularly  suitable  liquid  for  the  microdrops  because  of  its  low  vapor  pressure, 
high  surface  tension,  and  polarity.  To  improve  the  characteristics  of  the  microdrops. 
It  is  often  necessary  to  pre-coat  the  glass  cover  slide  with  a  thin  layer  of  silicone 
Dricote,  a  standard  anti-wetting  agent  for  PH  electrodes,  proved  to  be  adequate 
for  this  purpose.  Both  the  reduced  wetting  and  the  surface  tension  help  the 

glycerol  to  form  a  round  lens  with  a  relatively  short  focal  length.  This  conveniently 

produces  a  system  with  a  low  f  number.  The  microdrop  of  glycerol  is  positioned 
on  top  of  the  sample  cover  slide  so  that  the  illumi.''  ting  beam  first  impinges  on 
the  object,  then  on  the  microdrop  -  200  pm  closer  to  the  detector.  Since  the 
microdrop  is  a  positive  lens,  the  reference  point  source  is  formed  even  closer  to 

the  detector.  The  transverse  and  longitudinal  resolution  of  the  camera  system  is 

dependent  on  the  angular  aperture  subtended  by  the  hologram.  Consequently,  it  is 
necessary  that  the  cone  angle  of  the  diverging  reference  and  object  waves  be 
large  enough  to  produce  Interference  over  an  aperture  which  is  sufficient  to  give 
the  desired  resolution  in  the  image.  Hence,  the  f  number  of  the  microdrop  lens 

determines  the  usable  f  number  of  the  camera 

The  microdrops  are  formed  using  a  glass  micropipette  which  is  mounted  on  an 
X-Y-Z  micropositioner.  The  sample  and  micropipette  lip  are  viewed  through  a 
dissecting  microscope  while  manipulating  the  micropipette  into  position  with  the 

micropositionor.  When  the  lip  of  the  micropipette  has  been  positioned  properly,  n 
is  brought  almost  in  contact  with  the  glass  slide.  A  slight  pressure  at  the  bacK 
of  the  pipette  dispenses  the  desired  quantity  of  glycerol.  The  tip  is  then  slowly 

lifted  away  from  the  slide,  allowing  the  droplet  to  form. 

The  droplet  size  is  an  important  factor  in  the  formation  of  a  good  hologram 

because  it  is  primarily  the  diameter  of  the  drop  that  determines  the  brightness  of 
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the  reference  wave.  The  larger  the  diameter  of  the  drop,  the  larger  the  portion 
of  the  illuminating  wavefront  which  is  intercepted,  and  consequently,  the  brighter 
the  reference  wave.  In  order  to  form  a  hologram  with  a  good  signal-to-noise 
ratio^,  it  is  necessary  to  have  the  reference  wave  brightness  approximately  equal 
to  the  object  wave  brightness.  Depending  on  the  type  of  object  being  imaged,  we 
have  used  microdrops  ranging  between  20  and  200  pm  in  diameter. 

C.  Numerical  Reconstruction 

For  the  sake  of  clarity,  it  is  necessary  to  define  a  series  of  parallel  planes: 
(1)  the  hologram  plane  is  the  sensitive  surface  of  the  CCD;  (2)  the  reference 
plane  contains  the  focal  point  of  the  droplet;  and  (3)  an  object  plane  is  any  plane 
we  wish  to  reconstruct,  generally  a  plane  intersecting  the  specimen. 

The  basic  method  of  image  reconstruction  used  for  all  results  presented  in 
this  article  is  a  fast  Fourier  transform  (FFT).  However,  some  pre-processing  of 
the  holographic  data  is  performed  in  order  to  form  the  best  image.  Since  the 

hologram  contains  three-dimensional  data,  a  specific  object  plane  must  be  chosen 
as  the  focal  plane  for  the  reconstructed  image.  If  an  FFT  is  performed  on  the 
raw  holographic  data,  the  object  plane  will  be  the  same  as  the  reference  plane. 

Therefore,  the  hologram  formed  in  the  Fourier-transform  geometry  can  be  thought 
of  as  containing  an  implicit  lens  which  causes  it  to  form  a  focused  image  at  the 
reference  plane.  If  a  different  object  plane  is  chosen,  it  is  necessary  to  add  a 

second  lens  to  the  hologram  such  that  the  combined  action  of  the  implicit  lens 

and  the  additional  lens  cause  the  image  to  be  in  focus  at  the  desired  object 
plane.  This  second  lens  can  be  added  numerically  by  multiplying  the  hologram  by 
a  phase  factor  ({>  given  by 


4)  (x,y)  -  exp 


nr 


2  2 
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where  x  and  y  are  coordinates  in  the  hologram  plane  whose  origin  is  at  the 
center  of  the  hologram  and  fQ  is  the  desired  focal  length  of  the  additional  lens. 
The  focal  length  a  can  be  determined  from  knowledge  of  the  position  of  the 
reference  plane  and  the  position  of  the  desired  object  plane  by^ 
be 

fa  =  -  (2) 

b  -  c 

where  b  is  the  distance  from  the  reference  plane  to  the  hologram  plane  and  c  is 
the  distance  from  the  object  plane  to  the  hologram  plane. 

The  next  preprocessing  step  is  to  numerically  clip  the  corners  of  the  hologram 
and  remove  the  central  region,  thus  converting  the  originally  square  hologram  into 
an  annulus.  Removal  of  these  portions  of  the  data  reduces  noise  in  the  reconstructed 
image  and  involves  a  negligible  loss  of  image  quality.  The  holographic  data  are 
also  scaled  and  converted  to  a  floating  point  format  prior  to  performing  the  FFT. 

Post-processing  of  the  reconstruction  can  also  improve  the  quality  of  the 
image.  First,  the  central  portion  of  the  image  is  eliminated,  since  it  contains 
only  a  reconstruction  of  unwanted  scattering  sources  on  or  near  the  reference 
point  source  and  is  overwhelmingly  bright.  Secondly,  the  image  brightness  can  be 
rescaled.  The  best  images  are  given  by  the  power  spectrum  of  the  holographic 
data,  namely,  the  natural  log  of  the  sum  of  squares  of  the  real  and  imaginary 
parts  of  ..I  Fourier  transform.  The  image  data  must  then  be  converted  from 
floating  point  to  integer  format  for  display. 

D.  Resolution 

We  measured  the  microholograph  transverse  resolution  by  direct  observation 
of  the  spread  in  the  image  of  a  very  fine  line,  i.e.,  the  transverse  point  spread 
function.  The  line  was  produced  by  side  illumination  of  a  sharp  edge,  specifically, 
the  edge  of  one  of  the  figures  in  a  microscope  test  reticle.  The  image  of  this 
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edge  appears  as  a  line  because,  for  the  camera  geometry  presented  here,  only 
light  scattered  by  the  object  is  collected  to  form  an  image  (a  form  of  dark-field 
illumination).  The  width  of  this  line  as  observed  at  high  magnification  with  the 
same  sort  of  dark-field  illumination  in  a  compound  microscope,  was  seen  to  be 
much  less  than  1  pm.  The  hologram  of  the  line  was  recorded  and  reconstructed 
in  the  usual  manner.  An  intensity  distribution  at  a  representative  cross  section 
perpendicular  to  the  line  was  then  generated  at  small  spatial  increments.  This 

intensity  distribution  is  shown  in  Fig.  (3).  Remarkably,  the  full  width  at  half 

maximum  is  about  1.40  pm,  less  than  three  times  the  wavelength. 

The  standard  expression  for  diffraction  limited  resolution  of  a  microscope  is; 

1.22  X 

6  =  -  (3) 

2n  sin  9 

where  (n  sin  0)  is  the  numerical  aperture  (N.A.)  of  the  imaging  system.  This 

expression  is  valid  since,  for  the  camera  geometry  described  herein,  the  hologram 
fringe  spacing  is  always  greater  than  the  pixel  size  of  the  detector,  and  hence, 
the  image  resolution  is  not  limited  by  the  detector  pixel  size.  In  this  case,  the 
refractive  index  is  n  =  i,  and  at  f/2  the  N.A,  =  0.24.  Thus,  using  Eq.  (3),  <5  = 

1.29  pm.  This  compares  favorably  with  our  measured  value  for  the  point  spread 
function  of  1.40  pm. 

E.  Image  Reconstructions 

Although  we  have  examined  a  wide  variety  of  specimens,  we  show  here  an 
image  of  a  Parascaris  univalens  (Ascaris)  larva  because  it  demonstrates  several  of 
the  capabilities  of  our  microholograph  simultaneously. 

Figure  (4)  shows  the  reconstructed  image  from  a  hologram  taken  at  f/2  of  a 
section  of  the  Ascaris  larvae  in  its  maturation  stage.  Since  this  sample  is  not  a 
strong  scatterer,  the  reference  microdrop  had  to  be  only  ~  25  pm  in  diameter  in 
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order  to  make  the  reference  wave  and  object  wave  brightness  comparable.  Visible 
in  the  image  Is  the  cuticle  of  the  Ascaris  as  well  as  individual  cells.  Cell  walls 
are  clearly  visible  and  In  some  cells  nuclear  structure  is  evident. 

A  detailed  section  of  the  Ascaris  containing  10*  spatial  resolution  elements 
(10®  X  10*  X  10®)  in  the  volume  is  shown  In  Fig.  (5)  for  four  different  positions 
in  depth  of  the  focal  plane.  The  set  of  100  slices  was  reconstructed  from  a 
single  hologram.  The  spatial  range  covered  by  the  full  set  of  100  slices  (depth- 
wise)  is  72  pm;  consequently,  each  slice  is  separated  by  0.72  pm.  Hov/ever,  the 
actual  thickness  of  the  sample  is  ~  15  -  16  pm  as  measured  with  a  research-grade 
compound  microscope.  The  two  slices  shown  in  Fig.  (5b)  and  (5c)  represent  the 
lower  and  upper  surfaces  of  the  sample,  respectively,  and  are  separated  by  16.56 
pm.  Figure  (5a)  and  (5,b)  correspond  to  focal  planes  well  below  and  above  the 
boundaries  of  the  object.  For  Fig.  (5a)  the  focal  plane  is  24.48  pm  below  Fig. 
(5b),  and  for  Fig.  (5d)  the  focal  plane  is  24.48  pm  above  Fig.  (5c).  This  is  because 
longitudinal  blurring  of  the  image  introduces  ambiguities  which  can  be  partially 
alleviated  by  focusing  well  outside  the  sample  voijme.  This  quality  of  longitudinal 
resolution  is  comparable  to  what  we  observe  with  a  com.nound  microscope  with  an 
0.25  N.A.  (-  f/2)  objective. 

Blurring  of  the  image  is  primarily  due  to  two  factors:  (1)  a  longitudinal 
spatial  resolution  which  is  lower  than  that  governing  the  transverse  coordinates 
and  (2)  aberration  from  coma  that  is  inherent  in  the  use  of  a  Fourier-transform 
for  image  reconstruction.  The  changing  character  of  the  image  is  clear.  The 
central  nuclear  region  is  visible  and  the  structure  in  the  double  wall  can  also  be 
seen.  These  features  vary  with  the  frame  shown,  since  they  change  with  depth, 
V/e  have  found  that  a  video  representation  assembled  from  the  separate  100  slices, 
which  shows  the  continuous  variation  of  these  images,  is  a  more  informative 
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method  to  demonstrate  the  three  dimensional  character  of  the  microholographic  camera. 

Those  famiiiar  with  the  lore  of  holographic  microscopy  will  note  that  our 
reconstructions  are  conspicuously  devoid  of  speckle.  Speckle  is  not  a  probiem 
with  dark-field  illumination  because  the  image  is  produced  by  scattered  light. 

III.  CONCLUSION 

A  microholographic  camera  capable  of  three-dimensional  imaging  has  been 
constructed,  tested  and  applied  to  the  imaging  of  biological  materials.  This  system 
has  recorded  holograms  of  test  patterns  and  biological  specimens  with  transverse 
resolutions  near  the  diffraction  limit  at  f/2  with  514.5  nm  light.  The  three 
dimensional  character  of  the  recorded  images  has  also  been  demonstrated.  This 
instrument  can  be  directly  applied  to  imaging  in  the  x-ray  range  when  suitable 
sources  become  available. 
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FIGURE  CAPTIONS 


Fig.  (1):  Schematic  of  Fourier-transform  microholographic  camera  configured  for 

operation  in  the  visible  range. 

Fig.  (2):  Side-view  schematic  of  the  target  system  prepared  for  visible  microholog¬ 
raphy. 

Fig.  (3):  Transverse  point  spread  function  (PSF)  measured  by  imaging  the  edge 

of  a  microscope  test  reticle.  The  curve  is  comprised  of  250  data  points 
reconstructed  with  a  transverse  spacing  of  10  nm. 

Fig.  (4):  Reconstructed  image  of  an  Ascaris  in  maturation  stage.  Visible  are  the 
cuticle  and  individual  cells.  Cell  diameter  is  ~  40  pm.  The  nuclear 
region  visible  within  some  ceils  is  ~  20  pm  in  diameter. 

Fig.  (5):  Local  region  of  Ascaris  image  involving  10*  resolution  elements  (io’  x 

10®  X  10®).  The  region  is  centered  around  one  cell  adjacent  to  the 

cuticle.  Features  arising  from  the  depth  resolution  are  visible  as  the 

focal  plane  is  positioned  at  different  depths  (a)  -  (d).  See  text  for 
discussion. 
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APPENDIX  D 


Fourier-transform  holographic  microscope* 

W.S.  Haddad,  D,  Cullen.  J.C.  Solem,  J.  W.  Longworth, 
A,  McPherson,  K.  Boyer,  and  C.K.  Rhodes 

MCR  Technology  Corporation 
PO.  Box  10004,  Chicago,  IL  60610-0084 


1.  ABSTRACT 

A  camera  system  suitable  for  mlcrohotography  has  been  constructed,  Tested,  and  applied  to 
the  imaging  of  biological  materials,  The  design  of  this  instrument  Is  compatible  with  operation 
over  a  very  wide  spectral  range  spanning  from  visible  to  x-ray  wavelengths.  In  order  to 
evaluate  its  properties,  visible  light  Fourier  transform  microholograms  of  biological  samples  and 
other  lest  targets  have  been  recorded  and  digitally  reconstructed  using  a  glycerol  microdrop  as 
a  reference  wave  scatterer.  Current  results  give  a  resolution  of  -  4  X  with  X  «  514.5  nm. 

2.  INTRODUCTION 

Advances  in  high  speed  computation  and  electronic  detection  have  acceleraled  the  development 
of  practical  instruments  for  mlcroholography.  The  concept  discussed  herein  promises  a  new 
technology  for  three-dimensional  visualization  of  condensed  matter  that  Is  applicable  to  a 
spectral  range  spanning  from  visible  (-700  nm)  to  x-ray  (-  0.6  nm)  wavelengths.  Several  Important 
results  in  this  area  have  already  been  achieved  by  Howels  et  al.\  D.  Joyeux  and  F.  Polack^, 
McNulty  et  al.^,  J.  C.  Solem  <?t  al.***®,  and  London  et  al.^. 

Descriptions  of  the  Fourier-transform  system,  target  design,  and  related  numerical  simulations 
have  been  given  in  previous  work®"^®.  Tests  of  the  operation  of  the  system,  which  hive 
enabled  the  details  of  its  imaging  properties  to  be  studied,  have  been  performed  with  vis'bie 
radiaiion,  The  experiments  in  the  visible  region,  In  addition  to  allowing  optimization  ol  ih* 
instrument,  have  independently  demonstrated  that  this  Is  an  attractive  technique  for  visualizations 
in  fields  such  as  embryology  where  three-dimensional  information  is  of  particular  Interest. 

3.  VISIBLE  LIGHT  MICROHOLOGRAPHY 

3.1.  Apparatus 

A  schomalic  lor  the  holographic  microscope,  configured  for  operation  with  visible  radWitio" 
is  shown  in  Fig,  (1),  The  beam  path  through  the  target  is  vertical  so  that  the  target  holder  s 
horizontal  and  gravity  can  be  used  to  hold  the  sample  and  reference  wave  scatterer  In  place 
Several  objects  have  been  used  as  reference  wave  scatterers.  These  Include  GRIN  rod  lenses 
glycerol  microdrops,  and  reflecting  microspheres.  The  mIcrospheres  represent  the  optimum  chr^ro 
for  use  v/ith  soft  ’'-rays'^®,  and  the  corresponding  x-ray  target  design  Is  described  elsewhere”^ 
Glycerol  microdrops  have  performed  very  well  in  experiments  with  visible  light,  A  magnii'eo 
side  view  schematic  of  the  target  used  for  visible  light  exposures  is  Shown  In  Fig.  (2). 

The  basic  functions  of  the  apparatus,  such  as  shutter  liming  and  triggering,  are  conif.jiio-j 
by  a  PC,  The  hoiogram.s  are  recorded  electronically,  digitized,  and  stored  by  the  PC  f-e 
data  are  then  transferred  by  an  elhernet  imk  to  a  Stardent  GS  2000  supergraphIc  worksiauo" 
v/hcre  Ihe  image  reconstruction  is  performed  numerically.  The  reconsirucled  Imago  can  then  i:* 
manipula'ed,  scaled,  and  displayed  on  a  high  resolution  graphics  monitor.  A  schematic  ol  mu 
computer  system  Ls  shov/n  in  Fig,  (3|. 


*  maiiiiscri  pi  yubmittud  (or  llic  SPI  s  r./mpos  lum  on  Electronic  Imaging:  Science  S 

Tuclinolof'y,  2/^  Fubru/iry  -  I  M.arch  1991,  Sm  ,  CA ,  to  be  publi.shed  in  the  rrm-rr.ii 

of  the  Syinpo.sium  on  Elect rmitc  Imaging:  Science  and  Technology,  SPIE,  Bellingham, 


Fig.  (I):  Schematic  of  apparatus  for  visible/UV  Fourler-transform  microholographic  microscope 
configured  for  operation  in  the  vi.$ible  range. 


Fig,  (2):  Magnified  side-view  schematic  of  the  target  system  prepared  for  visible  microholography 
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Fig.  (3):  Diagram  of  computer  system  and  electronics  for  acquisition  and  reconstruction  of 
microholograms. 

3.2.  Hologram  detector 

The  holograms  are  recorded  on  a  charge-coupled  device  (CCD).  The  chip  has  a  sensitive 
area  of  18  mm  x  18  mm  containing  -  4  x  10*  pixels  (2048  x  2048).  The  pixels  are  square  wuh 
a  dimension  of  9  ym.  For  these  experiments,  a  protective  window  Is  placed  over  the  CCD 
which  carries  a  small  central  opaque  square  serving  as  a  blocker  for  the  unscaltered  Incident 

beam.  Specially  modified  drive  electronics  were  designed  for  this  system  Including  a  separate 
headboard  which  carries  the  CCD  and  allows  for  remote  mounting  of  the  detector.  Access  through 
the  headboard  to  the  back  of  the  CCD  also  is  available  to  permit  cooling  of  the  chip.  The  pixel 
clock  rate  (data  readout  rate)  has  also  been  reduced  from  -  20  MHz  to  2  MHz  to  allow  the 

system  to  function  well  when  cooled  to  -  -40®  C.  Such  cooling  would  allow  direct  Integration 
of  photons  on  ifie  chip  for  periods  of  up  to  about  one  hour  without  a  significant  accumulation 
of  dark  counts, 

3.3.  Liquid  microdrop  lenses 

A  liquid  microdrop  of  glycerol  has  been  found  to  serve  well  as  a  reference  wave  scaue'er 
for  visible  radiation.  The  drop  is  placed  on  the  glass  slide,  close  to  the  edge  of  the  sarr,pi(i 
and  nearest  the  region  of  interest  in  the  sample,  as  shown  in  Fig.  (2).  The  droplet  acts  as  t 
lens  intercepting  a  portion  of  the  Illuminating  plane  wavefront  and  converting  It  Into  a  diverging, 
nearly  spherical  wave.  The  focal  point  of  this  Ions  acts  as  the  reference  point  source,  rrie 

imaginary  plane  parallel  to  ihe  detector  plane  (hologram  plane),  and  containing  the  reference 

point  source,  will  be  referred  to  as  Ihe  reference  plane.  The  object  can  be  divided  up  into  a 
series  of  planes  (object  planes)  which  lie  parallel  to  the  hologram  plane.  In  most  cases. 
reference  plane  is  located  a  significant  distance  away  Irom  the  object  planes.  The  usual  sUu-riion 
is  that  the  microdrop  is  placed  on  top  of  the  sample  coverslide,  so  that  the  Illuminating  beam 
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firsi  impinges  on  the  object,  then  a  short  distance  later  (-  200  closer  to  the  detector),  the 
microdrop.  Sine  the  microdrop  is  a  positive  lens,  the  reference  point  source  is  then  formed 
even  closer  to  the  detector. 


Glycerol  is  a  suitable  material  for  the  microdrops  because  of  its  low  vapor  pressure,  high 
surface  tension,  and  polarity.  To  improve  the  characteristics  of  the  microdrops.  H  is  nUen 
necessary  to  pre-coat  the  glass  coverslide  with  a  thin  layer  of  silicone.  Dricote,  a  standard 
anil-wettlng  agent  for  PH  electrodes,  has  been  found  useful  for  this  purpose.  The  rodijr.od 
welting  and  the  viscosity  help  the  glycerol  form  a  round  lens  with  a  relatively  short  focal 
length.  This  produces  a  system  with  a  low  F  number. 

The  microdrops  are  formed  by  the  use  of  a  glass  micropipette.  The  mlcroplpette  Is  mounted 
on  an  X-Y-Z  microposilioner  and  the  operator  views  the  sample  and  micropipette  tip  through  a 
dissecting  microscope  while  manipulating  the  micropositioner.  When  the  tip  has  been  posit’onod 
properly,  it  is  brought  almost  into  contact  with  the  glass  slide,  and  a  slight  pressure  is  put  it 
the  back  of  the  pipette,  thus  dispensing  the  desired  amount  of  liquid.  The  tip  is  then  slowly 
brought  vertically  away  from  the  slide,  allowing  the  drop  to  form. 

The  droplet  size  is  an  important  factor  in  the  lormallon  of  a  good  hologram,  Thrs  is 
because  it  is  primarily  the  diameter  of  the  drop  that  determines  the  brightness  of  the  toforanro 
wave.  The  larger  the  diameter  of  the  drop,  the  larger  the  portion  of  the  lilumlnaling  wavefront 
which  is  intercepted,  and  therefore,  the  brighter  the  reference  wave.  In  order  to  form  a 
hologram  with  a  good  signal-to-noise  ratio,  it  is  necessary  to  have  the  reference  wave  brighmoss 
approximately  equal  to  the  object  wave  brightness.  Depending  on  the  type  of  object,  we  hjv* 
used  microdrops  ranging  between  20  pm  and  200  pm  in  diameter.. 


3.4.  Reconstruction  technique 

The  method  of  image  reconstruction  used  for  all  results  presented  In  this  article  is  has-': « 
a  fast  Fourier  transform  (FFT)  However,  some  pre-processing  of  the  holographic  data  is  pnr'  •>-< 
in  order  to  form  the  best  image.  A  specific  plane  parallel  to  the  hologram  plane  must  be  c*"'- 
as  the  focal  plane  for  the  reconstructed  image.  This  plane  will  be  referred  lo  as  if’e  « 
plane.  If  an  FFT  Is  performed  on  the  raw  holographic  data,  the  image  plane  will  be  tf^e  ^  v- 
as  the  reference  plane.  Therefore,  the  hologram  formed  in  the  Fourier  transform  geomp'w 
be  thought  of  as  containing  a  lens  which  causes  it  to  form  a  focused  image  al  the  n 
plane.  If  a  different  image  plane  is  chosen,  it  Is  then  necessary  to  add  a  second  lens  'u  ’ 
hologram  such  that  the  combined  action  of  the  implicit  lens  and  the  additional  lens  cau^«  • 
image  to  be  in  focus  at  the  desired  plane.  To  do  this,  the  hologram  must  be  multiplied 
phase  factor  L(xj^,  y^,)  given  by 

Vhl  ■  ‘  ’'I'’. 

v/hfjre  x^,  and  y^  are  coordinates  in  the  hologram  plane  whose  origin  js  at  the  center  ■* 
hologram. 

The  next  preprocessing  step  is  lo  clip  the  corners  of  the  hologram  and  remove  the  ’ 
region,  thus  converting  the  originally  square  hologram  into  an  annulus.  Removal  oi 
portions  of  the  data  reduces  noise  in  the  reconstructed  image  without  loss  of  image  o  i 
The  holoaraphic  data  are  also  scaled  and  converted  to  a  floating  point  formal  prior  to  pertor-^ 
the  FFT, 

Post  processing  of  the  reconstruction  can  also  Improve  the  quality  of  the  image  '  • 
the  central  portion  of  the  image  Is  eliminated,  since  it  contains  only  a  reconstfucllon 
unscaitered  beam  and  would  be  overwhelmingly  bright.  Secondly,  the  Image  brightness  ri'’ 
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r^'caled.  Th®  bast  Images  are  given  by  the  "power  spectrum*  of  the  holographic  data,  namely, 
the  natural  log  of  the  sum  of  squares  of  the  real  and  Imaginary  parts  of  the  Fourier  transform. 
The  image  data  must  then  be  converted  from  floating  point  to  Integer  format  for  display. 

3.5.  Results  of  Image  reconstructions 

3.5,1.  Swordtail  melanoma 


All  the  holograms  described  below  were  produced  with  the  use  of  514.5  nm  radiation  from 
an  argon  ion  laser.  Figure  (4)  shows  a  hologram  of  swordtail  (Xlphophorus  maculatus  x  heiierl) 
malignant  melanoma  and  the  reconstructed  Image.  Also  shown  are  Images  taken  with  a  standard 
low  power  optical  microscope  with  both  white  light  and  laser  Illumination.  Clearly  visible  in 
Fig,  (4a)  is  the  -  130  pm  glycerol  droplet  used  to  produce  the  reference  wave.  This  sample 
was  provided  to  us  by  R.  8.  Sellow  of  Brookhaven  National  Laboratory. 


(c)  (d) 

Fig,  (4):  Fish  melanoma  histological  sample  with  glycerol  mlcrodrop  reference  wave  source  (a) 
Image  of  the  target  taken  with  a  standard  optical  microscope  and  white  light  lllumlnaiion 
(b)  Same  as  (a)  except  with  Ar‘  (514,5  nm)  laser  Illumination,  (c)  Digital  recording  of 
hologram,  (d)  Numerical  reconstruction  ol  the  Image;  resolution  s  4  pm. 

3.5,2.  Air  Force  target 

Figure  (5)  Is  a  reconstructed  Image  from  a  hologram  ol  a  U.S.  Air  Force  resolution  tarq'’' 
The  full  reconstruction  is  shown  In  Fig,  (5a),  and  a  zoomed  Image  of  the  region  containing 
smallest  features  is  shown  in  Fig.  (5b).  These  slots  have  width  equal  to  the  slot  spacing  and 
for  (he  smallest,  this  dimension  is  ~  2  pm.  Since  the  features  of  the  target  are  holes  m  a 
metallic  coating,  they  scatter  the  Incident  radiation  by  diffraction.  Therefore.  H  Is  only  me 
edges  of  the  features  that  show  up  In  the  reconstruction.  In  principle,  the  Information  about 
the  central  regions  of  the  target  features  is  present  at  the  center  of  the  hologram,  however 
this  information  would  not  be  available  because  of  the  overwhelming  brightness  of  the  unscaitered 
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laser  beam.  In  Fig,  (5b),  it  is  possible  to  resolve  the  edges  of  both  the  vertical  and  horizontal 
sets  of  slots  lor  the  smallest  group  (~  2  pm  separation).  Since  X  «  514.5  nm,  a  conservative 
estimate  for  the  resolution  of  this  Instrument  Is  -  4\.  The  theoretical  resolution  limit  for  such 
an  imaging  system  operated  at  f/2  Is  -  f  X,  a  value  which  compares  favorably  with  our  result. 


Fig,  (5);  Reconstructed  Image  from  a  hologram  of  a  U.S.  Air  Force  resolution  target  taken  at 
f/2  with  514,5  nm  light.  A  glycerol  microdrop  was  used  as  the  reference  wave 
scatterer,  (a)  Full-size  reconstructed  image,  (b)  I6x  zoom  from  (a)  of  the  smallest 
slots  (2  pm  separation  between  each  edges). 

3.5.3.  Ascarls 

Figure  (6)  shows  the  reconstructed  image  from  a  hologram  taken  of  an  Ascarls  In  Its 
maturation  stage,  also  at  f/2.  Since  this  sample  Is  not  a  strong  scatterer,  the  reference  microdrop 
had  to  be  only  -  25  pm  In  diameier.  Visible  In  the  Image  Is  the  cuticle  of  the  Ascarls  as  well 
as  Individual  celts.  Cell  walls  are  clearly  visible  and  In  some  cells  nuclear  structure  Is  evident. 

4.  CONCLUSIONS 

A  holographic  microscope  capable  has  been  constructed,  tested,  and  applied  to  the  Imaging 
of  biological  materials.  This  system  has  recorded  holograms  of  test  targets  and  biological 
specimens  with  resolutions  in  the  range  of  -  4  X  with  514.5  nm  light.  This  Instrument  can  be 
directly  applied  to  Imaging  in  the  x-ray  range  when  suitable  sources  become  available. 
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Fig,  (6):  Reconstructed  Image  of  an  Ascarls  In  maturation  stage.  Visible  Is  the  cuticle  and 
individual  cells.  Cell  diameter  Is  ~  40  pm.  Nuclear  region  visible  within  some  ceili  is 
~  20  pm  in  diameter. 
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APPENDIX  E 


Prospects  for  x-ray  holography  with  free  electron  lasers 
J.C.  Solem*.  K.  Boyer,  W.  Haddad,  and  C.  K.  Rhodes 

MCR  Technology  Corporation 
P.O.  Box  10084 
Chicago.  IL  6061CH0084 

1.  ABSTRACT 

We  review  the  technical  advantages  offered  by  x-ray  holographic  microscopy  for  imaging 
the  structure  of  living  biological  specimens.  We  discuss  the  wavelength,  coherence,  energy, 
and  pulse-length  requirements  and  conclude  that  these  could  be  met  by  free-electron  laser 
architectures  of  the  near  future.  We  also  show  that  Fourier-transform  holography  using  a 
reference  scattering  sphere  is  the  best  optical  configuration  for  a  practical  Instrument. 

2.  INTRODUCTION 


X-ray  microholography  will  ultimately  allow  three-dimensional  imaging  of  living  organisms 
with  a  spatial  resolution  far  beyond  the  reach  of  optical  microscopes.  In  the  wavelength 
region  between  the  K-edges  of  carbon  and  oxygen,  we  find  a  high  contrast  between  protein 
and  water  that  will  allow  imaging  of  the  structural  elements  of  cells  without  staining.^ 
Because  of  the  quantum  nature  of  any  probe  on  the  <  100-X  level  of  resolution,  the  specimen 
will  be  severely  damaged  in  the  imaging  process.^  If  the  x-ray  energy  car.  be  delivered  in  a 
vary  short  pulse,  however,  a  snapshot  view  of  the  specimen  in  the  living  state  can  be  obtained, 
even  though  the  specimen  is  obliterated  in  the  process.^  An  ultrashort  pulse  x-ray  laser  could 
provide  the  illumination  for  such  snapshot  microholography.  Another  option  is  the  free  electron 
laser,  if  the  coherence  and  energy  requirements  can  be  met. 


3.  WATER  WINDOW 


Fiaure  1  shows  the  absorption  length  (e-folding  distance)  in  normal  density  protein  and 
water.^  Between  the  K-edges  of  carbon  and  oxygen,  there  is  a  region  of  high  contrast;  water 
is  transparent  and  protein  is  opaque.  This  wavelength  interval  has  become  known  as  the 
water  window  and  the  maximum  contrast  appears  at  the  K-edge  of  nitrogen.  Most  of  the  structural 
elements  of  a  cell  are  made  of  protein.  Therefore,  within  the  water  window,  most  biologically 
interesting  features  are  visible  without  staining.  This  compelling  fact  has  provided  the  principal 
impetus  for  the  emergence  of  the  field  of  x-ray  microscopy. 

Determining  exposure  requirements  for  x-ray  microscopy  from  the  data  In  Fig.  (1)  's  fairly 
straightforward.  For  holography,  however,  the  problem  is  a  bit  more  subtle.  About  a  decade 
ago  one  of  the  present  aulhors^'^  calculated  the  exposure  requirements  from  the  cold  opacities 
alone.  An  approximation  to  the  Mie  scattering  problem  for  a  semitransparent  sphere  in  a 
vacuum  revealed  that  the  scattering  cross  section  could  be  estimated  by  the  geometric  cross 
section  times  the  square  of  the  beam  extinction. ^  Since  that  time,  quantum  dispersion  theory 
has  been  applied  to  the  soft  x-ray  regime  to  give  both  raal  and  imaginary  atomic  scattering 
factors  and  researchers  have  solved  the  Mie  problem  for  protein  sphere  immersed  in  water, 
including  both  absorption  and  dispersion.®  The  results  have  revealed  that  the  best  wavelength 
for  holography  is  not  at  the  nitrogen  K-edge,  as  in  microscopy,  but  just  below  the  carbon  K- 
edge,  a  region  actually  below  the  threshold  of  the  water  window. 

There  appears  to  be  no  impediment  to  free-electron  lasers  producing  coherent  radiation  within 
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Fig.  1.  X-Ray  absorption  lengths  for  protein  and  water  at  1  gm*cm’*. 

the  water  window.  The  shortest  wavelength  achieved  to  date  is  2,400  A  by  the  Soviet  storage¬ 
ring  FEL  at  Novosibirsk,^  but  many  of  the  requisite  technologies  to  permit  extension  below 
100  A  using  RF  linear  accelerators®  are  now  in  hand.®  The  foremost  of  these  is  the  photoelectron 
injector.® 

Projecting  from  current  rf  linac  and  x-ray  optical  technology,®  either  an  FEL  oscillator  or 
a  single-pass,  self-amplifled-spontaneous-emlssion  amplifier  may  qualify.®”^ ^  Both  conceptual 
designs  could  be  operated  at  the  third  harmonic,  for  which  high  gain  could  be  achieved  with 
sufficiently  large  undulator  parameter.  In  one  speculative  single-pass  design, tuning  to  132 
A,  and  using  the  third  harmonic,  gives  about  three  orders  of  magnitude  more  gain  than  tuning 
directly  to  44  A.  An  oscillator  can  produce  the  requisite  power  with  higher  emittance  than 
required  for  a  single-pass  amplifier.  For  example,  an  FEL  oscillator  using  a  750  period  undulator 
(1.6  cm*period‘‘)  and  35%-reflectance  multilayer  mirrors  could  produce  44  A  radiation  at  5  * 
10*  W  in  2.5  x  10*^’  s  pulses  with  peak  beam  current  of  400  A  and  normalized  electron  beam 
emittance  of  lOit  mm»mr.  A  single-pass  laser  with  the  same  characteristics  would  require  an 
emittance  of  4Tr  mm*mr  and  an  undulator  twice  as  long.^^ 

4.  HOLOGRAPHY 

Holography  uses  interference  with  a  reference  wave  to  record  the  amplitude  and  phase  of 
the  specimen-scattered  wave.  From  the  symmetry  properties  of  Maxwell's  equations  it  can  be 
shown  that,  '*  the  fields  are  known  on  a  surface,  they  are  also  known  in  three-dimensional 
space.  Thus,  it  is  possible  to  reconstruct  rather  complete  three-dimensional  information  about 
the  geometric  structure  of  the  specimen  from  its  scattered  wave.  Although  the  configurations 
of  holographic  apparatus  are  diverse  and  varied,  all  recording  techniques  divide  into  two  ba,sic 
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categories:  Fresnel  transform  and  Fourier  transform. 

4.1.  Freanel-tranaform  holography 

In  Fresnel-transform  holography,^®  the  reference  waves  are  planar.  The  specimen  is 
illuminated  by  the  same  beam  of  radiation  that  supplies  the  reference  waves:  the  specimen 
simply  stands  between  the  radiation  source  and  the  recording  surface.  TSte  wider  the  ang!e 
ftom  the  specimen  to  the  edge  of  the  reference  illuminated  area  on  the  recording  surface,  the 
better  the  resolution.  The  same  ccne-angle  dependence  applies  to  an  optical  .riicrcscope.  The 
cone~angle  dependence  further  implies  that  the  resolutiori  depends  on  the  coheience  of  the 
radiation.  Spatial  coherence  must  be  maintained  over  ths  extent  of  the  reference  illuminated 
area  of  the  recording  surface  and  the  temporal  coherence  length  must  be  larger  than  the 
difference  between  the  perpendicular  distance  frrm  the  recording  surface  to  the  spR::imen  an'j 
the  distance  fiom  the  edge  of  the  refeience  illuminated  area  to  the  specirrioii.^^ 

A  more  severe  limitation  derives  from  the  intrinsic  resolution  of  the  recording  .-wedium. 
The  resolution  of  a  Fresne!  transforn?  holograpf<  is  limited  to  twice  the  gain  size  of  the 
recording  Tiedium.^^  This  limitation  makes  '.v-ray  hr'lo.^rsphy  exceeding!*/  difficult,  oven  with 
the  finest  grain  commercial  films.  The  deus  ex  machira  is  photoresist.  C*iveloped  for  microHthog- 
raphy,  photoresists  are  photocenaitive  elchaole  polymers  with  no  significant  gniin  structure.^*^ 
X-rays  will  break  cross  links  In  these  polymers  making  them  mote  susceptinte  to  etching  in 
the  regions  ci  x-ray  exposure.  After  etching,  the  resist  rcsembies  a  relief  map,  with  lowest 
elevations  corresponding  to  greatest  exposuras.  The  spatial  resolution  of  photoresists  (more 
properly,  x-ray  resists)  Is  dependenc  on  the  wavelength  of  the  x-rays  to  which  it  is  exposed. 
For  long  wavelengths,  it  is  llmite:,!  by  diffraction  within  the  resist  Itself.  For  short  wavelengths, 
it  is  limited  by  the  lange  of  sacondary  electrons,  which  can  also  break  cross  links  and  thereby 
cause  blurring.  The  optimum  waveiength’^  Is  about  50  A.  serendipltously  close  to  the  lower 
edge  of  the  water  window,  The  best  spatial  resolution  of  r^sKst  at  the  optimum  wavelength  is 
50  to  100  A. 

Two  groups  are  presently  making  holoaram''^  of  bMogical  specimens  using  the  Fresnel 
transform  technique.  An  Arr.erican  group'®  working  at  the  National  Synchrotron  Light  Source 
of  Brookhaven  National  Laboratory  uses  a  spv.)c!ally  equipped  beamllne  with  an  undulator  lor 
high  brightness  ir?  the  water  window.  After  exposing  and  etching  the  resist  on  which  their 
hologram  is  recosded,  they  shadow  .ho  contours  with  gold.  They  then  read  out  the  hologram 
with  a  transmission  electron  microscooe,  digitize  it,  and  reconstruct  It  with  a  computer.  A  French 
group^®  using  a  synchrotron  at  Orsay  records  their  holograms  on  photoresist  in  a  simiia. 
manner.  They  rd-co.nstiuct  the  hologram  by  passing  a  He-N«?  (6,328  A)  laser  rough  tre 
etched  resist,  whose  contours  produce  locally  differing  phase  shifts.  This  produces  .  magnified 
image  of  the  specimen.  They  can  accomplish  the  same  phase  shifting  effect  by  aluminizing  the 
surface  of  the  resist  and  reflecting  tiie  laser  beam. 

4.2  Fourier  transform  holography 

Fourier  transform  holography^  ^  does  not  suffer  Ir.-'m  the  resolution  limits  of  the  recording 
medium  that  affect  Fresnel  transform  holography.  Tho  Fourier  technique  uses  curved  rather 
than  planar  reference  waves.  In  the  far  Held  using  ideally  spherical  reference  waves,  the 
technique  maps  increments  of  distance  In  the  specimen  into  spatial  frequencies  at  the  recording 
surface  —  thus  the  name  “Fourier  transform*.  In  the  far  field,  the  spatial  relationship  among 
the  fringes  simply  scales  with  distance  from  the  specimen.  Thus,  thn  hologram  can  bo  faithfully 
registered  by  recording  media  of  ativ  resolution  by  simply  adjusting  the  distance  from  the 
specimen  to  the  recording  surface.  A  charge-coupled  device  (CCD)  array  is  particularly  attractive 
in  this  application.  CCDs  can  be  made  to  work  with  high  efficiency  in  the  water  window,’® 
and  their  signal  can  be  digitized  directly  for  simplified  computer  reconstruchon. 
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Curlou&ly,  small  Increments  of  distance  (minute  features  of  the  specimen)  are  mapped  into 
low  spatial  frequencies,  white  large  Increments  of  distance  (gross  features)  are  mapped  into 
high  spatial  frequencies.  This  counterintuitive  result  can  be  reconciled  with  familiar  optics  by 
noting  that  the  diffractive  resolution  limit  of  the  holograph  is  again  determined  by  the  angle 
over  which  the  hologram  1$  recorded.  A  larger  angle  will  allow  a  better  measurement  of  lower 
spatial  frequencies. 

The  Fourier-transform  technique  has  less  stringent  coherence  requirements  than  the  Fresnel 
transform  technique.'^  Spatial  coherence  is  reduced  because  the  coherent  pa*"*  the  wavefront 
need  only  span  the  specimen  and  the  lens  that  is  creating  the  reference  wav'..>  and  these  may 
be  placed  quite  close  together.  Temporal  coherence  is  reduced  because  the  difference  in 
pathlength  from  the  reference  source  and  the  specimen  Is  always  small.  Figure  2  compares 
the  coherence  requirements  for  the  two  techniques  as  a  function  of  specimen  mass  over  1 1 
orders  of  magnitude.  A  free-electron  laser  will  have  no  difficulty  meeting  these  coherence 
requirements,  but  the  Fourier-transform  technique  Is  clearly  less  demanding. 
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A  principal  difficulty  in  Fourier  transform  holography  is  the  generation  of  spherically 
diverging  reference  waves.  A  straightforward  technique  Involves  the  use  of  a  lens.  The  water 
window  is  too  far  from  any  strong  resonances  for  a  refractive  lens  to  be  used,  so  wo  must 
use  a  diffractive  lens  in  the  form  of  a  Fresnel  zone  plate.  It  can  be  shown,  however,  that  the 
resolution  is  then  limited  to  the  same  order  as  the  finest  spacing  of  the  zone  plate.  In  other 
words,  we  need  to  bo  able  to  fabricate  x-ray  optics  to  the  same  precision  as  we  wish  to 
resolve  the  specimen.  Present  technology  is  limited  to  a  few  hundred  ar^gstroms.^^  If  a 
pinhole  with  a  diameter  less  than  the  finest  zone  spacing  is  placed  at  the  focus,  then  the 
resolution  is  on  the  order  of  the  pinhole  diameter.  It  may  be  possible  to  fabricate  pinholes 
with  diameters  somewhat  smaller  than  the  finest  zone-plate  spacing.  The  real  solution  to  the 
reference  wave  dilemma  Is  a  reflecting  sphere,  which  will  be  discussed  In  detail  later. 

Despite  these  difficulties,  the  Brookhaven  group  is  now  generating  Fresnel  transform  x-ray 
holograms  using  a  zone  plate  and  a  CCD  as  a  recording  surface. The  CCD  output  can  be 
directly  digitized  for  computer  reconstruction. 


5.  EXPOSURE 

Holography  results  from  the  symmetry  properties  of  Maxwell's  equations — classical  electro¬ 
dynamics.  VVe  know  this  is  only  an  approximation  to  quantum  theory.  While  the  wavefunction 
is  analytic  and  deterministic,  any  measurement  is  inherently  statistical.  We  must  count  particles 
just  as  in  nuclear  physics  experiments.  Since  the  energy  carried  by  a  single  photon  in  the 
water  window  Is  more  than  a  hundred  times  the  energy  carried  by  an  optical  photon,  the 
damage  to  the  specimen  is  enormous.  Furthermore,  we  are  trying  to  resolve  much  finer 
features  than  in  the  optical  range  and  each  of  these  features  must  scatter  a  statistically  significant 
number  of  photons.  Thus  for  diffraction  limited  resolution,  the  damage-  is  increasing  very 
roughly  as  the  fourth  power  of  wavelength.  For  an  electron  microscope,  the  damage  for  the 
same  resolution  is  even  greater.^ 

Table  I  lists  some  specimens  and  the  dosage  required  to  Image  them  to  a  resolution  of  too 
A  with  signal-to-noise  of  5,  using  x-rays  at  the  nitrogen  K-edge.  Somewhat  higher  dosages 
are  required  for  the  virons  because  there  is  less  contrast  between  DNA  and  protein.^  ^  Higher 
doses  are  required  for  the  larger  cells  because  of  their  physical  thickness.  For  most  cells,  a 
few  times  10*  rad  Is  the  required  dosage.  Using  the  optimal  wavelength  just  below  the  water 
window  would  reduce  these  exposures  about  an  order  of  magnitude,  but  the  enormity  of  the 
dose  is  very  clear.  To  put  this  in  perspective,  lo'*'  rad  would  make  y:>u  very  sick  and  to*  rad 
would  be  surely  fatal.  The  sturdiest  known  living  thing  is  the  bacterium  Micrococcus  radiodursns 
which  has  been  alleged^^  to  survive  3  x  10*  rad.  Surely  no  specimen  can  survive  the  exposure 
The  best  we  can  do  is  try  to  record  Its  true  image  when  it  was  in  the  living  state,  and  this 
is  best  accomplished  by  taking  a  snapshot. 

6.  SNAPSHOT  HOLOGRAPHY 

X-ray  flash  pictures  or  snapshots  seem  essential  to  imaging  living  specimens.  The  required 
dosages  are  so  high  that  we  can  expect  the  specimen  to  explode  owing  to  the  sudden  heating 

If  the  pulse  is  short  enough,  the  Image  can  be  captured  from  fhe  plasma  ghost  of  the  specimen 

before  there  is  sufficient  bulk  motion  to  blur  the  features  we  are  seeking  to  resolve.^  Analytic 
expressions  for  hydrodynamic  blurring  ov/ing  to  the  ex,  loslon  of  semiopaque  features  are  useful 
for  estimating  pulse  length  requirements.  To  achieve  a  resolution  6,  a  specified  number  of 

photons  must  be  coherently  scattered  from  a  volume  j*  during  the  exposure  time  At.  and  no 

dimension  of  the  specimen  can  be  allo'wed  to  change  by  more  than  6  during  At.  To  examine  a 
feature  with  resolution  6  <<  d,  where  d  is  the  transverse  dimension  of  that  feature  me 
maximum  exposure  time  is 
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Specimcii 

Mass(g) 

Oosage(rad) 

Specimen 

Mass(g) 

Oosage(rad) 

Echerichia  coll 
phase  4>x174 

10**’ 

2  X  10*’ 

Anthrax 

bacterium 

10'** 

4  X  10* 

Herpes 

viron 

10'** 

2  X  10** 

Red  blood  cell 
(human) 

10'** 

5  X  10* 

Vaccinia 

viron 

10**' 

6  X  10** 

White  blood  cell 
(human) 

10'* 

4  X  10* 

Mycoplasma 

pneumoniae 

io“- 

3  X  10* 

Amoeba 

(dysentery) 

10'* 

2  X  10* 

Escherichia  coli 
{immature) 

10'** 

3  X  10* 

Smooth  muscle 
cell 

10*’ 

4  X  10** 

Escherichia  coll 
(mature) 

10'*’ 

3  X  10* 

Paramecium 

(protozoan) 

10'* 

2  X  10** 

Table  I.  Dosage  to  Obtain  100  A  Resolution  with  S/N  -  5  at  N  K-Edge, 
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where  m  is  the  mass  of  the  feature  being  examined,  Og  is  its  absorption  cross  section.  I  s  <*'•; 
x-ray  intensity,  and  F  is  a  dimensionless  coefficient  accounting  for  equalion-of-state  parameters 
and  the  rate  at  which  transport  processes  remove  energy  from  the  exploding  feature  to 
surrounding  material.  If  6  >>  d  is  acceptable,  then  ^t  should  be  multiplied  by  6/d  4">*r 
substituting  a  different  dimensionless  coefficient  G. 


The  number  of  useful  photons  elastically  scattered  from  a  single  resolution  eiemt>'>t  s 
approximated  by 


eX 

Ns—  pIag6*Al 
fim 


where  p  is  the  density,  Og  is  the  elastic  scattering  cross  section  of  the  feature,  a  s  '*•; 
wavelength,  and  c  is  the  quantum  efficiency  of  the  recording  medium.  This  number  •- 

statistically  significant  compared  to  photon  noise.  Combining  Eqs,  (1)  and  (2)  gives  the  roj  ,  r-; 
intensity  as  a  function  of  resolution. 
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Intensity  increases  very  rapidly  as  better  resolution  is  sought.  For  typical  b<oi  ;  . 
specimens,  intensities  on  the  order  of  10*’  W-cm  '*  with  pulse  lengths  on  the  order  of  '  ; 
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are  required  for  a  resolution  of  100  A.  The  projected  designs  for  free-electron  lasers  will 
have  no  difficulty  meeting  these  requirements. 

7.  REFERENCE  SCATTERING  SPHERES 

The  resolution  dilemma  encountered  with  the  Fresnel  zone  plate  reference  source  in 
Fourier  transform  holography  Is  solved  by  using  a  reference  scattering  sphore.^^  The  zone- 
plate  holograph  is  limited  to  a  resolution  about  the  same  as  the  finest  zone-plate  spacing 
because  the  largest  angle  covered  by  the  reference  illumination  is  determined  by  that  spacing. 
A  tiny,  highly  reflective,  smooth  sphere,  which  can  be  used  to  replace  the  zone  plate  as  the 
reference  source,  will  not  suffer  this  limitation. 

Reference-sphere-Fourier-transform  holography  also  offers  a  host  of  advantages  when  an 
x-ray  laser  is  used:  (1)  the  sphere  can  provide  reference  illumination  over  a  sufficiently  large 
angle  to  realize  diffraction-limited  resolution:  (2)  the  sphere  is  appropriate  to  the  narrow 
beam  of  a  laser  in  that  it  is  small  and  can  be  placed  quite  close  to  the  specimen;  (3)  orientation 
of  the  sphere  is  unimportant;  and  (4)  the  sphere  is  rather  inexpensive  to  fabricate  in  large 
quantities,  so  unlike  the  zone  plate,  it  is  unimportant  whether  the  sphere  Is  damaged  or 
destroyed  on  each  exposure.  Figure  3  is  a  schematic  of  how  the  reference  sphere  might  be 
used  in  practice.  The  most  logical  recording  device  is  a  CCD,  the  output  of  which  can  be 
digitized  for  direct  computer  analysis.  The  scattering  sphere  will  not  produce  spherical  reference 
waves,  there  will  always  be  an  angle-dependent  phase  shift.  Thus  It  will  not  generate  a  true 
Fourier'  transform,  even  in  the  far  field.  However,  any  scattering  shape  can  be  used  as  long 
as  we  know  exactly  what  It  Is  and  how  it  is  oriented.  A  set  of  basis  functions,  analogous  to 
the  sinusoidal  basis  functions  of  Fourier  reconstruction  which  incorporates  the  appropriate 
correction,  can  be  generated  by  the  computer  and  used  to  recover  the  Image. The  image 
can  also  bo  reconstructed  optically  using  a  scaled  reflecting  sphere  and  a  computer  modified 
hologram  that  compensates  for  the  angular  dependence  of  reflectivity.  The  hologram  might 
even  be  an  array  of  liquid  crystal  optical  gates  for  direct  interfacing  with  the  computer  and 
‘real-time*  reconstruction. 


X-R»  V  Laser 


X-Ray 

CCD 

Mosaic 


Fig.  3.  Schematic  of  a  reference-sphere-Fourier-transform  holographic  microscope. 
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That  it  is  possibie  to  manufacture  a  sphere  with  the  requisite  surface  smoothness  was  not 
obvious  in  the  beginning  of  our  research  program.  We  found  that  many  of  the  spheres  manufac¬ 
tured  by  iasar  abiation  had  rather  crystalline  surfaces  resulting  from  slow  coding.  However, 

In  examining  soma  hdlow  spheres  that  were  manufactured  as  laser  fusion  targets,  we  found 
surfaces  so  smooth  that  no  structure  could  be  resdved  with  an  electron  microscope.  The 
smoothness  was  attributed  to  rapid  coding  of  the  thin  layer,  which  left  the  metal  in  a  glassy 
state  rather  than  forming  crystals.^^'^^  Smooth  spheres  have  also  been  manufactured  in  very 
small  sizes  at  Lawrence  Livermore  Laboratory.  Finally,  if  there  Is  a  problem  with  the  metallic 
spheres,  a  clear  alternative  exists.  For  example,  we  can  use  readily  available  glass  spheres, 
which  will  have  the  necessary  smoothness,  and  then  coat  them  by  metal  evaporation,  a  process 
which  Is  known  to  produce  amorphous  layers  with  a  smoothness  better  than  3  A.  Actually, 
recent  work^^  on  x-ray  mirrors  has  demonstrated  the  ability  to  achieve  a  microroughness  of  - 
2  A.  Processes  to  provide  uniform  coating  have  been  developed.  The  spheres  also  turn  out  to 
be  exceedingly  round  owing  to  rapid  damping  of  multipole  oscillations  on  the  time  scale  of 
cooling.  This  removes  completely  the  necessity  for  measuring  departure  from  sphericity  and 

for  orientation  of  the  sphere. 

We  surveyed  the  periodic  table  for  elements  of  optimum  specular  reflectivity  to  serve  as 

reference  scatterers  in  the  water  window.  Unless  a  compound  can  significantly  enhance  the 
density  of  an  element,  the  pure  element  will  be  the  best  reflector.  We  used  the  tables  of  real 
and  Imaginary  atomic  scattering  factors^^  calculated  from  quantum  dispersion  theory^^  to 
obtain  the  reflectivities  shown  in  Fig.  4.  Nickel  Is  found  to  be  the  best  reflector  at  the 
water  window  threshold  (carbon  K-edge).  If  we  are  able  to  precisely  select  our  wavelength,  it 
would  behoove  us  to  use  the  optimum  wavelength  just  below  the  threshold,  so  nickel  is  the 
element  of  choice.  Furthermore,  osmium  Is  only  slightly  superior  to  nickel  and  suffers  from  a 
variety  of  fabrication  and  handling  difficulties.  Figure  5  shows  the  reflectivity  of  nickel  as  a 
function  of  grazing  angle  at  44.7  A.  It  provides  adequately  uniform  reference  illumination  out 
to  an  angle  of  about  26°.  This  is  sufficient  to  obtain  good  resolution  in  both  the  transverse 
and  longitudinal  directions.  The  transverse  resolution  is  given  by  5t  s  X/26,  while  the  longitudinal 
resolutions^  Is  given  by  6i  =  X/4e^  Table  11  shows  the  resolutions  that  can  be  obtained  at 
the  carbon  K-edge  and  the  nitrogen  K-edge. 

Assuming  an  average  specimen  diameter  and  as  compact  a  spacing  as  possible  between  the 
sphere  and  specimen,  the  laser  beam  would  have  to  Illuminate  about  an  area  uf  about  to  * 
cm^  This  means  the  free-electron  laser  would  have  to  supply  about  10‘*  J  In  a  single  pulse, 
which  is  quite  reasonable. 


Nickel 

Osmium 

44,7  A 

31.6  A 

e 

30° 

22* 

«t 

43  A 

41  A 

41  A 

63  A 

Table  II.  Spatial  Resolution  at  carbon  and  nitrogen  K-edges. 
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Fig.  4.  Specular  reflectivities  of  various  elemerK^  at  selected  grazing  angles. 


Fig.  5.  Specular  reflectivity  of  nickel  as  a  function  of  grazing  angle. 

8.  CONCLUSIONS 

If  the  proposed  design^'®  can  be  reduced  to  practice,  no  Impediment  is  seen  to  using  a 
free  electron  laser  to  obtain  water  window,  x-ray  holographic  images  of  biological  specimens 
The  best  arrangement  is  to  use  the  Fourier  transform  technique  with  a  nickel  reference  scattering 
sphere  and  a  CCD  to  register  and  record  the  interference  pattern. 
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The  University  of  Illinois  at  Chicago 

Laboratory  for  Atomic.  Molecular,  and  Radiation  Phyeics  (M/C  273) 

Oepadment  of  Physics 

College  of  Liberal  Arts  and  Sciences 

Box  4346,  Chicago,  Illinois  60680 

(312)  996-4866 


Aug.  6,  1990. 

MCR  Technology  Corporation 
P.  O.  Box  10084 
Chicago,  IL  60610-0084 

Dear  Dr.  Rhodes: 

Enclosed  please  find  the  report  for  the  ‘Propagation  Studies  for  High  Brightness 
X-Ray  Sources*  (subcontract  090189-SDI-i). 

The  main  findings  of  this  study  are  (1)  capillary  tubes  are  effective  In  enhancing 
harmonic  conversion  as  well  as  par.iimetric  processes,  (2)  propagation  properties  of 
intensity  laser  beam  are  influenced  by  electrons  produced  from  Ionizations  to  the 
extend  that  the  laser  fluence  inside  the  capillary  tube  falls  off  rapidly. 

One  of  the  aspects  in  this  study,  namely  large  fluctuation  in  the  harmonic  Intensity, 
have  not  been  studied  in  great  detail.  However,  one  can  speculate  that  it  is  directly 
related  to  the  focal  intensity  distribution  of  the  pump  beam.  Recently,  we  have 
developed  a  method  of  examining  the  focal  intensity  distribution  from  any  focussing 
system.  From  these  Investigation,  although  incomplete,  reveals  that  unless  the  laser 
beam  is  properly  aligned  to  the  focussing  optics,  the  focal  intensity  distribution 
can  be  distorted  severely  by  astigmatism  and  spherical  aberration.  This  immediately 
suggests  that  the  performance  of  harmonic  and  parametric  conversion  can  be 
expected  to  improve  substantially  when  these  optical  distortions  can  be  controlled 
judicially.  Furthermore,  availability  of  this  diagnostic  tool  will  allows  us  to  perform 
better  controlled  experiments  to  Investigate  spectral  stability  Issue  of  harmonic 
radiation  in  the  future. 

It  has  been  a  pleasure  working  on  this  contract.  If  there  Is  any  questions  about 
the  report,  please  feel  free  to  contact  me. 

Sincerely  yours. 


Ting  Shan  Luk 
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Abstract 


The  advent  of  subpicosecond  laser  systems  capable  of  achieving  focal 
intensities  exceeding  10**  W/cm^  allow  consideration  of  a  fundamentally  new  mode 
of  electromagnetic  propagation,  namely  channeled  propagation,  In  which  both 
relativistic  electronic  motions  and  electron  expulsion  from  the  region  of  high  laser 
intensity  produces  a  channel  which  could  be  useful  as  a  high-«spect-ratlo 
longitudinally  pumped  x-ray  laser  source.  Two  classes  of  experiments  to  investigate 
longitudinal  pumping  were  performed.  The  first  used  small  bore  capillary  tubes  to 
investigate  the  propagation  of  the  pump  laser  through  an  extended  gas  target  and  to 
examine  the  radiation  omitted  from  this  target.  The  second  class  of  experiments 
Investigated  the  laser  beam  propagation  through  a  static  gas  by  photographing  the 
breakdown  of  the  target  gas  in  the  laser  focal  region  in  order  to  determine  the 
spatial  location  of  the  emitted  radiation. 
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Channeled  Propagation  Simulation  In  Caplllarv  Tube  Targets 


I.  Introduction 

The  Interaction  of  high  intensity  subpicosecond  laser  systems  with  gaseous 
targets  allow  the  production  of  plasmas  that  are  characterized  by  a  low 

temperature,  perhaps  as  low  as  a  few  tens  of  electron-volts,  and  a  density  that  may 
approach  the  critical  value.  A  fundamentally  new  mode  of  electromagnetic  propaga¬ 
tion  is  theorized  to  arise  in  this  species  of  plasma,  namely,  channeled  propagation 
produced  when  the  laser  intensity  exceeds  a  value  approximately  one  half  of  the 

Compton  intensity,^  which  for  a  KrF*  laser  is  4.5  x  10**  W/cm*.  Because  of 

charge  displacement,  the  boundary  of  a  channel  can  be  associated  with  a  sharp 
change  of  index  of  refraction  of  the  target  medium  resulting  from  relativistic 
motions  and  the  expulsion  of  electrons  from  the  region  of  high  laser  Intensity.  If  a 
channel  is  formed,  a  self-contained  cylindrical  region  of  target  material  would  be 
created  which  could  serve  as  a  high-aspect-ratio  source  of  stimulated  radiation, 
Recent  dynamical  calculations^  have  Indicated  the  existance  of  such  a  stable  mode 

of  propagation.  The  ability  to  propagate  high  intensity  radiation  via  channeled 
propagation,  and  thus,  longitudinally  pump  a  target  material  is  the  main  motivation 
for  pursuing  the  Investigations  described  below. 

It  has  been  established  that  atoms  can  be  put  Into  high  charge  states  through 
a  multiphoton  absorption  process  when  gaseous  atomic  targets  are  exposed  to  strong 
electromagnetic  fields  created  by  high  intensity  laser  systems. Experiments  in 
this  laboratory  have  demonstrated  that  the  resulting  ion  can  be  left  in  an  excited 
state,  although  whether  the  excitation  is  a  direct  consequence  of  the  multiphoton 
process  or  is  the  result  of  colHsional  recombination  processes  is  yet  to  be 
determined.  What  is  certain  Is  that,  depending  upon  the  target  material,  copious 
quantities  of  short  wavelength  radiation  are  produced  from  the  small  interaction 
volume  of  target  gas  produced  from  a  pulsed  gas  jet  and  the  focal  volume  of  the 
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laser  system,  the  lighter  gases  up  through  Ne  being  the  noted  exceptions.^ 

The  following  premise  was  proposed  to  investigate  longitudinal  pumping  of  a 
gaseous  target  medium.  Assume  that  the  observed  radiation  from  a  pulsed  gas  jet 
target  is  produced  by  either  a  direct  multiphoton  process  or  by  collisional 
recombination  processes.  Further,  assume  that  this  radiation  can  be  made  to 
originate  from  an  extended  source  rather  than  the  typical  focal  volume  of  the 
laser.  Can  or  will  this  extended  emission  region,  which  would  be  generated  by 

confining  the  target  gas  inside  a  small  bore  capillary  tube,  exhibit  anr.plification? 

Longitudinal  pumping  is  accomplished  by  utilizing  a  small  bore  capillary  tube  to 
extend  the  length  of  the  target  gas  to  a  cylindrical  volume.  In  the  event  that 
recombination  is  the  primary  method  of  producing  excited  state  Ions,  the  capillary 
tube  walls  will  assist  in  rapidly  cooling^  the  plasma  produced.  Small  bore  glass 
capillary  tubes  are  used  because  the  smooth  interior  wall  has  a  high  reflectivity  at 
grazing  Incidence,  a  fact  which  allows  for  the  reflection  of  the  laser  beam  rather 
than  its  absorption,  much  like  the  anticipated  refractive  index  boundary  resulting 
from  channeled  propagation  at  high  laser  intensity  in  a  gaseous  target.  Emission 
will  then  presumably  occur  along  some  fractional  length  of  the  capillary  tube.  Only 
radiation  emitted  in  the  direction  parallel  to  the  capillary  bore,  or  parallel  to  the 

axis  of  propagation  and  not  absorbed  by  residual  gas,  will  be  able  to  escape  the 

capillary  tube  and  be  detected.  Whether  the  detected  radiation  is  the  result  of  a 
direct  multiphoton  absorption  process  yielding  an  excited  state,  and  therefore 

prompt,  or  If  It  is  due  to  collisional  recombination  mechanisms,  amplification  parallel 
to  the  capillary  lube  bore  should  result  if  absorption  from  the  ambient  material  is 
sufficiently  low. 

II,  Capillary  Transmission  Studies 

A.  Low  Power 

Preliminary  studies  considered  the  transmission  of  a  HeNe  laser  beam  focused 
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into  and  through  a  small  bore  capillary  tube.  These  studies  established  the 
methodology  o’  focusing  and  aligning  a  laser  through  a  small  bore  capillary  tube. 
Essentially,  the  capillary  tube  was  found  to  be  correctly  aligned  when  the 
transmitted  beam  produced  a  TMoo  image  on  a  screen  positioned  beyond  the  output 
side  of  the  capillary  tube.  Figure  1  illustrates  the  initial  setup  for  these  studies. 
The  capillary  tube  was  positioned  inside  a  small  vacuum  chamber  equipped  with  CaF, 
input  and  output  windows  and  pumped  to  a  base  pressure  of  a  few  mlllitorr.  The 
exact  base  pressure  was  not  important  as  long  as  it  was  much  lower  than  the  target 
gas  pressure  and  sufficient  not  to  affect  the  laser  beam  propagation,  especially 
when  tho  KrF*  laser  was  used.  Alignment  of  the  HeNe  beam,  and  the  KrF*,  was 
facilitated  by  placing  the  vacuum  chamber  with  the  enclosed  capillary  tube  upon  an 
X-Y-Z-0-4)  translation  stage.  A  HeNe  laser  was  focused  into  a  250  pm  bore 
capillary  lube  which  has  a  length  of  approximately  10  cm.  The  focusing  optic  was  a 
56  cm  focal  length  lens.  Since  the  HeNe  laser  beam  was  expanded  and  collimated  to 
a  1.3  cm  diameter  beam,  an  effective  f/42  focal  system  resulted.  The  transmission 
of  the  HeNe  beam  through  this  capillary  tube  was  measured  using  an  empty  capillary 
tube  as  well  as  the  capillary  tube  with  Nj  flowing  through  it.  Gas  was  supplied  to 
tne  capillary  tube  by  attaching  another  tube  to  a  small  hole  which  had  been  drilled 
into  tho  side  of  and  at  the  center  of  the  capillary  tube.  A  constant  gas  flow  was 
mainta  ed  while  measuring  the  transmission.  The  capillary  tube  was  translated 
along  the  laser  path  to  assure  that  the  entrance  of  the  capillary  tube  was 
positioned  at  the  focal  position  of  the  laser  beam.  The  target  gas  density  was 
varied  by  changing  the  gas  line  pressure  to  the  capillary  tube  from  10  psi  up  to  70 
psi.  The  transmission  of  HeNe  through  the  capillary  tube  for  a  target  gas  line 
pressure  of  50  psi  is  plotted  In  Figure  2  and  compared  to  a  static  air  fill  of  760 
torr  or  less.  Nearly  90%  of  the  HeNe  beam  was  transmitted  through  the  empty 

capillary  tube.  However  almost  97%  was  transmitted  through  the  capillary  tube 
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Vacuum  Gauge 


250  /Ltm  Bore 
Capillary  Tube 


Figure  1.  Initial  experlmentlal  setup  to  measure  the  transmission  of  laser  light 
through  a  small  bore  capillary  tube. 


Transmitted  Energy  (arh.  units) 


lli^  '  =  42 


G AS :  Nj  Lons  FocCil  Lengtfi  ->  56  cm 
Beam  Dfameter  13  cm 


Micrometer  Position  (Relative  Units) 


Figure  2.  A  relatively  high  percentage  of  HeNe  laser  light  was  transmitted 
through  the  250  pm  bore  capillary  tube.  When  gas  was  flowed  through 
the  tube,  the  transmission  Increased. 


68 


when  Nj  wee  flowing  through  it  and  being  pumped  out  of  the  encompassing  vacuum 
chamber.  This  increase  in  transmission  was  attributed  to  a  lensing  effect  of  the 
gas  jet  escaping  out  the  end  of  the  capillary  tube.  The  escaping  gas  served  to 
favorably  deflect  the  laser  beam  into  the  capillary  bore.  Increasing  the  gas  flow 
above  a  characteristic  maximum  value,  typically  40  psi,  but  which  depended  upon  the 
diameter  of  the  capillary  bore,  yielded  no  additional  transmission  due  to  this  lensing 
effect.  Changing  the  bore  diameter  of  the  capillary  tube  to  150  pm  resulted  In  a 
reduced  transmission,  perhaps  because  of  the  focusabillty  of  the  HeNe  laser  beam. 
However,  the  lensing  effect  became  more  pronounced  as  indicated  in  Figure  3  with 
tne  transmission  Increasing  from  72%  to  83%.  Smaller  bore  capillaries  were  not 
considered  because  of  an  inability  to  remove  all  debris  from  the  bore  resulting  from 
the  side  hole  diiliing  process. 

Having  demonstrated  that  a  HeNe  laser  beam  could  be  transmitted  through  a 
small  bora  capillary  tube  of  considerable  length,  propagation  of  the  KrF*  beam 
through  the  capillary  cube  was  considered.  Initially,  a  KrF*  beam  of  approximately 
100  pJ  was  focused  into  the  capillary  tube  using  the  same'  setup  as  described  In 
Figure  1,  Alignment  of  the  capillary  tube  bore  with  respect  to  the  KrF*  laser  beam 
was  found  to  be  critical  as  evident  by  the  reduced  quality  of  the  TM,,(,  mode.  Data 
was  taken  after  the  best  TMo,  mode  image  was  obtained.  Transmission  through  an 
empty  10  cm  long  capillary  tube  with  a  250  pm  bore  diameter  was  determined  to  be 
approximately  40%  to  50%.  The  transmission  dropped  rapidly  as  the  capillary  tube 
was  rotated  away  from  the  maximum  transmission  orientation. 

B,  High  Power 

After  performing  several  experiments  involving  the  acquisition  of  emission 
spectra,  which  are  discussed  below,  the  transmission  of  the  smaii  bor©  capillary 
tubes  was  tested  using  approximately  100  mJ  of  KrF*  energy  from  the  amplifier 
Prometheus.^  The  focusing  optic  used  was  a  10  cm  diameter  lens  with  a  focal 
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Figure  3.  When  the  capillary  bore  was  reduced  to  150  pm,  the  transmitted 

HeNe  energy  was  reduced.  However,  the  effect  of  gas  flowing  through 
the  tube  was  increased. 
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length  of  100  cm.  Since  the  KrF*  laser  beam  from  Prometheus  has  a  rectangular 
2"  by  4"  shape,  the  focal  system  was  f/l0  to  f/20.  Two  sizes  of  capillary  tubes 
were  used,  40  mm  long  with  a  0.75  mm  bore  and  30  mm  long  with  a  0.25  mm  bore. 
Gas  was  admitted  into  the  capillary  tube  through  a  pulsed  valve  gas  jet  assembly 
similar  to  that  shown  in  Figure  4.  Later  modifications  allowed  capillary  tube 
lengths  from  2.5  cm  up  to  23  cm  to  be  tested.  Target  density  in  the  capillary  tube 
could  be  varied  in  two  ways.  One  could  change  the  backing  pressure  to  the  pulsed 
valve  or  one  could  change  the  time  delay  between  the  laser  fire  signal  and  the 

pulsed  valve  open  signal.  The  gas  flow  through  the  capillary  tube  was  monitored 

and  calibrated  using  a  differential  microphone.  Typical  responses  for  two  capillary 
tube  lengths  and  bores  are  illustrated  in  Figure  5, 

For  the  transmission  test  using  the  Prometheus  amplifier,  the  timing  was  set 

so  that  the  capillary  tube  was  completely  filled  with  gas  and  the  gas  flow  was 
nearly  a  maximum  at  the  capillary  tube  exit.  Figures  6-9  Illustrate  the  results.  As 

noted  In  Figures  6  and  7,  some  fraction  of  the  KrF*  laser  beam  could  be 

transmitted  through  all  of  the  test  gases  even  at  high  pulsed  valve  backing  pressure 

when  using  a  750  pm  bore  4  cm  long  capillary  tube.  However,  when  the  capillary 
tube  bore  was  reduced  by  a  factor  of  3,  to  250  pm,  and  the  length  was  decreased  to 
3  cm,  the  transmission  of  the  KrF*  beam  dramatically  changed.  Transmission 
through  the  empty  capillary  tube  dropped.  While  transmission  through  He  and  Me 
remained  high,  even  at  high  pulse  valve  backing  pressure,  transmls.sion  through  the 

heavier  gases  rapidly  dropped  to  zero.  This  effect  cannot  be  attributed  only  to  a 
refractive  effect  due  to  the  abundance  of  electrons  occuring  from  ionizing  Ar  and 

Kr  because  H,  behaved  like  Ar  and  Kr  in  attenuating  the  laser  beam  as  illustrated 
In  Figure  8.  The  transmission  of  the  KrF*  beam  through  an  empty  30  mm  long 
capillary  tube  with  a  bore  of  250  prn  was  measured  to  be  approximately  60%,  or 
roughly  equal  to  the  transmission  through  a  much  longer  capillary  tube  using  100  p,j 
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Figure  5.  Using  a  microphone,  the  gas  flow  through  the  caplllery  tube  could 
bo  monitored.  A  750  pm  bore  23  cm  long  capillary  with  50  psi  backing 
pressure  was  used  In  (a).  A  250  pm  bore  3  cm  long  capillary  with  40  psi 
was  used  in  (b).  Maximum  flow  occures  when  the  curve  crosses  the  x- 
axis.  The  x-axis  time  scale  is  500  psec  per  division. 
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Figure  6.  The  KrF*  beam  could  be  transmitted  through  a  750  pm  bore  4  cm 
long  capillary  tube  even  at  pulsed  valve  backing  pressure  o(  2500  torr  for 
Ne,  Nj  and  Ar.  Note  the  high  transmission  for  Ne. 
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Figure  7.  The  KrF  beam  could  be  transmitted  through  a  750  pm  bore  4  cm 
long  capillary  tube  even  at  pulsed  valve  backing  pressure  of  2500  torr  for 
He  and  Hj.  Note  the  high  transmission  for  He  and  that  the  transmission 
for  Hj  Is  simular  to  that  of  Ar. 
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Figure  8.  When  the  capillary  tube  bore  was  reduced  to  250  pm  and  the  length 
decreased  to  3  cm,  the  KrF  beam  transmission  was  significantly  reduced. 
While  He  and  Ne  still  transmitted  the  beam  Hj,  Ar  and  Kr  did  not. 
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Figure  9.  Schmatic  illustrating  the  laser  being  focused  into  a  capillary 

tube  of  length  I  and  bore  diameter  b.  The  focusing  optic  has  a  diameter 
of  d  and  a  focal  length  f.  The  first  reflection  occu''s  at  position 
r  Into  the  capillary  tube. 


of  energy  and  a  larger  f-number  focal  system.  Clearly,  the  KrF*  beam  is  being 
absorbed  into  the  capillary  tube  walls  at  each  reflection.  Finally,  note  that  addition 
of  gas  does  not  produce  a  lensing  effect.  Rather,  the  laser  transmission  through 
the  capillary  tube  decreases,  for  all  gases,  as  the  gas  flow  through  the  capillary 
tube  is  increased.  Only  by  increasing  the  capillary  tube  bore  diameter  to  750  pm 
does  the  laser  beam  transmission  impove.  A  defocusing  action  due  to  the  electron 
density  could,  at  least  in  part,  account  for  these  observations. 

C.  Capillary  Reflections 

Internal  reflections  are  important  since  the  laser  beam  is  partially  absorbed  at 
each  reflection.  Consider  Figure  9  in  which  a  lens  Is  used  to  focus  laser  light  ’ 
and  at  the  bore  of  a  capillary  tube.  Using  light  rays  to  define  similar  triatic 
the  first  reflection  is  found  to  occur  at  the  position 

r  =  bf/d 

where  b  is  the  capillary  bore  diameter,  d  is  the  laser  beam  diameter  and  f  is  the 
focal  length  of  the  lens.  Additional  reflections  will  occur  at  (2n-1)r  along  the 
length  of  the  capillary  tube  where  n  is  an  integer.  The  total  number  of  reflections, 
N,  in  a  capillary  tube  is  found  from 

(2N  -1)r  =  I 

where  I  is  the  total  length  of  the  capillary  tube.  In  the  two  transmission  examples 
cited  above,  there  will  be  2  reflections  for  the  40  mm  long  capillary  tube  and  4 
reflections  for  tho  30  mm  long  capillary  tube,  the  small  bore  capillary  tube 
experiences  more  reflections.  Following  several  light  rays  through  the  capillary 
tube  indicates  that  beyond  the  point  2r,  the  laser  beam  establishes  an  average 
intensity  over  a  diameter  of  approximately  b/2.  The  angle  of  reflection  Is  given  by 

6  =  tan-‘(d/2f). 

If  the  focusing  optics  result  in  an  angle  of  incidence  which  is  too  large,  absorption 
will  be  severe.  Denote  the  absorption  of  t*^e  laser  beam  at  a  reflection  as  T.  For 
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the  examples  of  Figures  6-8,  the  laser  beam  traveling  ttr-Owtgh  the  250  pm  bore 
capillary  tube  v/ould  experience  considerablly  more  absorption,  than  the  laser 
beam  traveling  through  the  750  pm  bore  capillary  tube. 

Sl?!ce  the  gas  flow  through  the  capillary  tube  will  be  important,  as  indicated 
by  the  strong  attenuation  of  the  KrF*  beam  when  using  Hj,  Ar  and  Kr,  as  well  as 
the  need  to  know  the  target  density,  a  model  was  established  which  assumes  that 
the  gas  flow  Is  viscous.  The  validity  of  this  model  was  checked  using  an 
interferometer  to  measure  the  rata  of  gas  flow  in  the  capillary  tube  and  round  to 
be  within  10%  of  the  calculated  value.  For  details  of  the  interferometer 
measurement  see  Appendix  A.  Further  testing  of  the  gas  flow  was  accomplished 

using  a  microphone  to  monitor  the  flow  through  the  capillary  tube  as  illustrated  in 
Figure  5.  From  the  model,  the  pressure  at  the  end  of  a  tube  is  related  to  the 
pressure  at  the  input  of  the  tube  by 

P-  «  -or/2  +  (or’  +  4Po7/“ 

where 

a  =  16r(V^l./a’, 

r\  is  the  viscosity  of  the  ga.s,  is  the  molecular  velocity  of  the  gas,  L  is  the 
length  of  the  tube,  and  a  is  the  bore  radius.  A  worked  example  is  presented  In 

Apoendix  B  indicating  that  target  densities  as  large  as  10‘*  cm"*  are  encountered. 
With  modest  ionization  of  the  target  gas  along  the  capillary  tube,  mther  large 

electron  densities  could  result. 

III.  Spectral  Studies 

A.  Normal  Incidenca  Spectrometer 

Spectra  from  various  gases  were  obtained  using  the  capillary  target 
geometry.  Two  seperate  data  acquisition  systems  were  employed,  a  1-meter  normal 
incidence  spectrometer  (NIS)  with  a  1200  I/mm  (.rating  and  a  l-meter  grazing 

incidence  spectrometer  (GIS)  with  a  300  I/mm  gratmg.  Each  spectrometer  employed 
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a  microchannel  plate  (MCP)  detector  mounted  at  the  exit  focal  plane  of  the 
spectrometer.  For  the  NIS,  a  charge  coupled  device  (CCD)  camera  was  coupled  to 
the  MCP  allowing  21-nm  regions  of  the  spectrum  to  be  simultaneously  collected. 

For  the  GIS  an  optical  multichannel  analyser  (OMA)  was  coupled  to  an  MCP 

mounted  tangen  'ally  to  the  Rowland  Circle  and  allowed  6-nm  regions  of  the 
spectrum  to  be  collected.  Spectral  resolution  on  the  NIS  was  0.9  nm  due  to  using 
the  capillary  tube  bore  as  the  spectrometer  entrance  slit.  The  GIS  gave  0.2  nm 
resolution.  Alignment  of  the  capillary  tube  bore  with  the  optical  axis  of  the 
spectrometer  was  achieved  using  an  X-Y-Z  manipulator  which  allowed  tilt  motion. 

The  gases  H,,  He,  N^,  Ne,  Ar,  Kr  and  Xe  were  used  as  target  materials. 
Spectra  between  35  nm  and  90  nm  were  obtained  for  al>  gases,  with  a  few  extending 
up  to  130  nm.  The  gases  Ar  and  N,  will  be  discussed  separately  with  the  remaining 
gases  being  summarized  together  for  completeness.  Three  general  areas  of 
investigation  were  examined:  (i)  the  spectrum  as  a  function  of  cap’Hary  tube 

length  and  bore  diameter;  (2)  the  spectrum  as  a  function  of  target  density;  and  (3) 
the  spectrum  as  a  function  of  laser  energy  and  confocal  parameter.  Most  of  the 
experlmential  effort  is  limited  to  the  first  two  areas.  Careful  studies  requiring 
changing  the  laser  energy  and  the  confocal  parameter  were  not  possible,  only  large 

changes  were  attempted  to  clearly  determine  if  an  effect  was  present. 

The  first  experlmential  configuration  utilized  the  NIS.  A  23  cm  long  capillary 
tube  with  a  0.75  mm  bore  diameter  was  positioned  at  the  Image  plane  and  served  as 
the  entrance  aperture  to  the  spectrometer.  Argon  was  selected  as  the  initial  target 
material  because  of  familiarity  with  the  VUV  spectrum.®  Figure  10  Illustrates  the 
typical  spectrum  obtained.  Two  features  are  noted.  First,  the  third,  fifth  and 
“reventh  harmonics  are  evident  although  the  seventh  Is  weak.  These  throe  lines  can 
be  used  as  an  approximate  wavelength  calibration  to  compensate  for  errors 
Introduced  by  using  the  capillary  bore  as  the  spectrometer  entrance  silt.  Second, 
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750  pm  bore  as  the  target  geometry.  The  KrF  laser  energy  was 
approxiTiatei'.  10'  mJ.  Focusing  was  with  an  f/25  to  f/10  optical  system. 
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all  of  the  remaining  spectral  lines  are  unidentified  lines,  and  the  spectrum  Is  not 
nearly  as  rich  as  that  obtained  previously  using  only  a  gas  jet  target.^  The 
strengths  of  the  signals  vary  considerably  with  the  third  harmonic,  at  nearly  2000 

times  stronger  than  the  fifth  harmonic,  being  the  strongest  feature.  No  lines  were 
observable  below  33  nm  because  scattered  light  from  zeroth  order  prevented 
observation  below  this  region. 

Argon  was  searched  at  longer  wavelength.  Figure  11  illustrates  an  Interesting 
feature  observed  at  si  24  nm  which  is  thought  to  be  Ar,*.  Since  this  transition  has 
a  natural  line  width  of  «  10  nm,  the  signal  In  Figure  11a  would  represent  a  gain  of 
«3.5  cm~^,  the  strongest  gain  yet  recorded  for  this  excimer  transition. As 
should  be  expected,  this  line  exhibited  a  strong  dependence  upon  laser  energy  and 
confocal  parameter  as  illustrated  in  Figures  lla-d. 

Figure  12  illustrates  the  lines  grouped  around  the  third  harmonic  In  the  Ar 
spectrum  as  the  target  pressure  is  increased  by  Increasing  the  backing  pressure  of 
the  pulsed  valve.  Figure  13  plots  the  signal  strength  as  a  function  of  backing 
pressure  as  the  pressure  is  increased  and  then  de'  eased  and  illustrates  that  the 
effect  is  reproducible.  Two  features  are  svideni.  First,  the  third  harmonic  has  a 

strong  target  density  dependence  as  evidenced  by  the  rapid  increase  In  signal  when 
the  backing  pressure  Is  increased  from  15  to  25  psi.  Second,  when  the  target 
density  reached  a  critical  value,  the  third  harmonic  signal  drops  rapidly  and  a  new 
signal  at  87  nm  rises.  The  association  between  these  two  spectral  signals  will  be 
discussed  later. 

As  noted  earlier,  the  target  density  wiihin  the  capillary  tube  could  be  varied 

other  than  by  changing  the  pulsed  valve  backing  pressure.  By  changing  trie  timing 
between  the  laser  fire  command  and  the  pulse  valve  open  command,  the  laser  pulse 
could  bo  made  to  arrive  at  any  time  durii'g  trie  gas  flow.  Figure  14  Illustrates  the 

change  In  signal  which  could  be  acfileved  by  altering  the  target  density  in  this 
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Figure  11.  Argon  spectrum  Illustrating  the  laser  energy  dependence  and 
confocal  parameter  dependence  of  the  Ar,  124-nm  line.  (a).  Laser 
energy  <  26  mJ.  Large  confocal  parameter  achieved  by  aperaturlng  down 
the  laser  beam.  (b).  Lasei  energy  =  25  mJ.  Full  laser  beam  aperature. 
(c).  Laser  energy  -  100  mJ,  Full  laser  beam  aperature,  (d).  Laser 
energy  =  140  mJ,  Full  laser  beam  aperature. 
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Figure  12.  Argon  spectrum  illustrating  the  growth  of  lines  near  the  third 
harmonic  a$  the  pulsed  valve  backing  pressure  is  Increased,  (a).  20  psi. 

(b).  30  psI.  (c).  40  psi.  Note  that  as  the  B7-nm  line  grows,  the  third 

harmonic  decreases. 
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Figuro  13.  Plot  of  the  signal  strength  vs  pulsed  valve  backing  pressure  of  the 
third  harmonic  and  ez-not  line  In  argon.  The  backing  pressure  was 
Increased  (dots)  and  then  decreased  (circled  dots)  to  demonstrate  that  the 
result  was  repeatable. 
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Figure  14,  Using  a  23  cm  long  capillary  tube  with  a  750  pm  bore,  the  region 
of  Ar  spectrum  around  the  third  harmonic  changes  as  the  target  density 
Is  changed  by  varying  the  delay  between  the  laser  fire  signal  and  the 
turn  on  of  the  pulsed  valve,  (a).  500  psec.  (b).  1000  psec.  (c).  1500 

psec. 
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fashion.  Note  the  presence  of  the  87  -nm  feature  as  seen  in  Figure  12  and  that 
the  third  harmonic  has  decreased  significantly  while  lines  at  lower  wavelength  have 
appeared  to  replace  it.  Changing  the  timing  was  a  more  accurate  and  reproducible 
method  of  changing  the  target  density  than  changing  the  pulsed  valve  backing 
pressure  on  the  present  setup. 

When  a  9  cm  long  capillary  tube  with  a  0.75  mm  bore  was  used  with  Ar  as  the 
target  gas,  the  basic  spectrum  was  not  changed.  However,  the  variation  witii  target 
density  did  become  more  pronounced.  Figure  15  illustrates  the  variation  in  the 
third  harmonic  region  and  the  87-nm  region  as  the  target  density  was  changed. 
Note  the  sharp  growth  of  the  signal  in  the  87-nm  region  and  its  near  equal 
intensity  with  the  third  harmonic  region  once  it  !s  established.  More  about  this 

variation  with  target  density  will  be  discussed  later  in  association  with  the  observed 
change  in  the  N,  spectrum. 

B.  Grazing  Incidence  Spectrometer 

The  instrumentation  was  changed  to  the  GIS  to  obtain  better  resolution  of  the 
observed  spectrum.  In  order  to  reduce  the  absorption  of  omitted  radiation  along 
the  capillary  tube  axis,  a  40  mm  long  capillary  with  a  0.75  mm  bore  was  used 

Since  no  new  spectral  lines  were  seen,  the  region  around  the  third  harmonic  was 

selected  for  a  more  detailed  study.  Figure  16  illustrates  the  results,  The  laser 
energy  was  -200  mJ  or  greater,  The  target  gas  density  was  changed  by  using  a 
capacitance  manometer  to  monitor  the  pulsed  valve  backing  pressure,  Accurate 

reproducable  low  backing  pressures  to  the  pulsed  valve  were  achleveable  by  thi,$ 
i.iethod.  The  time  delay  between  the  laser  fire  signal  and  the  pulsed  valve  open 
signal  was  held  constant.  It  was  selected  so  that  the  gas  flow  had  reached  the  end 
of  the  capillary  bore,  but  not  yet  reached  a  maximum  value.  Considerable  structure 
was  found  In  the  region  of  the  third  harmonic  as  the  target  gas  density  was 
changed.  None  of  the  lines  could  be  identified  with  the  known  spectrum  of  argon. 
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The  instrumentation  was  changed  to  the  GIS  to  obtain  better  resolution  of  the 
observed  spectrum.  In  order  to  reduce  the  absorption  of  emitted  radiation  along 
the  capillary  tube  axis,  a  40  mm  long  capillary  with  a  0.75  mm  bore  was  used. 
Since  no  new  spectral  lines  were  seen,  the  region  around  the  third  harmonic  was 
selected  for  a  more  detailed  study.  Figure  16  illustrates  the  results.  The  laser 
energy  was  *200  rnj  °r  greater.  The  target  gas  density  was  changed  by  using  a 
capacitance  manometer  to  monitor  the  pulsed  valve  backing  pressure.  Accurate 
reproducable  low  backing  pressures  to  the  pulsed  valve  were  achleveable  by  this 
method.  The  time  delay  between  the  laser  fire  signal  and  the  pulsed  valve  open 
signal  was  held  constant.  It  was  selected  so  that  the  gas  flow  had  reached  the  end 
of  the  capillary  bore,  but  not  yet  reached  a  maximum  value.  Considerable  structu'e 
was  found  in  the  region  of  the  third  harmonic  as  the  target  gas  density  was 
changed.  None  of  the  lines  could  be  identified  with  the  known  spectrum  of  argon. 
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Figure  15.  Using  a  9  cm  long  capillary  tube  with  a  750  pm  bore,  the  region  of 
Ar  spectrum  around  the  third  harmonic  exhibits  simular  structure  as 
Illustrated  In  Figure  14.  (a).  250  psec.  (b).  500  psec.  (c).  750  psec. 

(d).  1000  sec.  (e).  1500  psec.  (t).  2000  psec. 
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Flaure  16  Using  a  4  cm  long  capillary  tube  with  a  750  pm  ^or®,  a  rnoro 

^  detailed  Ar  spectrum  around  the  third  harmonic  Is  obtained  with  the  GIS 

?he  target  cS  s i  y  was  now  changed  by  varying  the  pulsed  valve  backing 

pmssure  (a).  200  torr.  (b).  325  torr.  (c).  550  torr.  (d).  760  torr 

(e).  1975  torr.  (f).  2475  torr. 
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The  capiltpry  tube  was  changed  to  one  which  was  30  mm  long  and  had  a  bore 
diameter  of  0.25  mm.  The  results  for  the  region  around  the  third  harmonic  are 
illustrated  In  Figure  17.  Again  the  structure  around  the  third  harmonic  is  distinct 
as  the  target  gas  pressure  is  increased.  With  the  smaller  bore  capillary  tube,  the 

structure  is  more  quickly  achieved  as  the  target  gas  density  is  increased  and  the 
structure  extends  over  a  greater  wavelength  range.  Once  the  extended  structure  is 
established,  the  observed  lines  exhibit  a  strong  dependence  upon  laser  energy. 
Rarely  do  all  of  the  lines  appear  together  during  a  single  laser  shot,  rather  they 

jump  in  and  out  depending  upon  the  laser  energy.  This  is  evident  in  Figures  17e 
and  17f  In  which  the  target  gas  density  changes  little,  yet  the  spectrum  changed 
noticeably. 

The  most  unusual  feature  found  in  the  Ar  spectrum  Involves  the  63-nm  line, 
see  f-igura  10.  This  line  was  found  to  exhibit  amplification!  While  the  amplification 
!s  small,  only  it  did  seem  to  lend  creedance  to  the  idea  that  this 

capillary  target  goosT.etry  could  be  longitudinally  multiphoton  pumped  to  produce 
lasing.  Not  only  did  this  line  exhibit  gain,  but  it  also  exhibited  line  narrowing. 

Figure  IP  illustrates  the  increase  in  signal  and  the  decrease  in  line  width  as  a 
function  of  pulse  valve  backing  pressure. 

The  difficulty  in  understanding  the  63-nm  line  In  Ar  is  manifested  by  the 

discovery  that  many  of  the  lines  observed  around  the  third  harmonic  and  87  nm 
could  fit  into  a  parametric  process  involving  six  wave  mixing.  For  details,  see 

Appendix  C.  Figure  19  illustrates  the  various  six  wave  mixing  processes  Identified. 

That  six  wave  mixing  had  never  been  reported  in  Ar  was  not  the  problem.  The  63- 
nm  line  playing  a  prominent  role  in  the  six  wave  mixing  process  does  present  a 
problem.  If  the  63-nm  line  Is  the  result  of  six  wave  mixing,  it  can  exhibit  gain, 

However,  it  Is  not  clear  that  such  a  parametric  process  would  exhibit  line 

narrowing.  The  measurement  was  repeated  with  the  same  results,  a  rapid  Increase 
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Figure  17.  Argon  spectrum  around  the  third  harmonic  using  a  3  cm  long 
capillary  tube  with  a  250  pm  bore.  Note  the  increase  In  structure 
compared  to  Figure  16.  The  pulse  valve  backing  pressure  is;  (a)  50  torr; 
(b)  290  torr;  (c)  620  torr;  (d)  1040  torr;  (e)  1950  torr;  (I)  2000  torr. 
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Figure  18,  The  appearent  gain  in  the  63-nm  argon  iine  is  compared  to  the 
known  argon  resonance  line  at  106.7  nm.  During  the  increase  in  63-nm 
signal,  this  line  exhibited  line  narrowing. 
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Figure  19.  Four  argon  energy  level  diagrams  illustrating  six-wave  mixing 
processes  which  could  account  lor  some  of  the  observed  radiation.  The 
dotted  arrow  in  (d)  is  a  predicted  but  unconfirmed  transition  at  5.78  pm. 
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in  signal  to  a  saturation  point  as  the  target  gas  density  is  increased  and  a 
corresponding  linewidth  decrease.  These  aspects  of  the  behavior  remain  to  be 
completely  understood. 

Nitrogen  was  an  important  target  gas  to  consider  for  these  capillary  tube 
target  experiments  because  of  the  results  obtained  from  muitiphoton  Ionization 
experiments  as  well  as  the  radiation  experiment  results  obtained  from  using  a  gas 
jet  target.''^”^^  Previous  measurements  on  molecules  suggested  that  molecular 
species  would  exhibit  interesting  properties  in  regard  to  the  production  of  excited 
states.  Moreover,  the  previous  Nj  work  implied  that  more  complicated  molecular 
species  containing  high  Z  material  symmetrically  surrounded  by  lighter  material 
could  be  better  candidates  for  target  gases. 

When  N,  was  used  as  the  target  gas  in  a  23  cm  long  capillary  tuba  with  a  0.75 
mm  bore,  qualitatively  similar  results  were  obtained.  The  spectrum  Illustrated  In 
Figure  20  shows  a  strong  third  harmonic,  a  weak  fifth  harmonic  and  a  few  weak 
unidentified  lines.  Regions  of  the  spectrum  not  illustrated  were  devoid  of  spectral 
lines.  Figure  21a-b  and  21c-e  illustrate  the  rapid  change  in  line  intensity  as  the 
target  density  is  changed  by  increasing  the  delay  between  the  laser  fire  signal  and 
the  pulsed  valve  open  signal.  The  results  are  the  same  as  observed  in  Ar, 
especially  as  pertains  to  the  region  around  the  third  harmonic. 

Later  a  9  cm  long  capillary  tube  with  a  0.75  mm  bore  was  used  as  the  target 
geometry.  The  typical  Na  spectrum  obtained  is  illustrated  in  Figure  22.  Basically 
there  is  no  change  in  the  spectrum  from  using  the  much  longer  23  cm  capillary  tube 
except  for  the  few  additional  longer  wavelengths.  Using  the  shorter  capillary  tube 
did  result  in  a  significant  change  in  the  observed  spectrum  when  the  target  density 
was  changed  by  altering  the  time  delay  between  the  laser  fire  command  and  the 
pulse  valve  open  command.  Figures  23  and  24  illustrate  the  point.  Figures  23a-b 
show  a  dramatic  change  In  the  third  harmonic  with  it  being  completely  replaced  by 
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Figure  20.  Nitrogen  spectrum  obtained  using  a  23  cm  long  capillary  tube  with 
a  750  pm  bore. 
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Figure  21.  Nitrogen  spectrum  obtained  using  a  23  cm  long  capillary  tube  with 
a  750  pm  bore.  Note  the  rapid  change  in  the  line  structure  as  the 
target  density  is  changed  by  altering  the  time  delay  between  the  laser 
fire  signal  and  the  pulsed  valve  open  signal,  (a).  500  psec. 

(b).  1500  psec.  (c).  500  psec.  (d).  750  psec.  (e).  1000  psec. 
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Figure  22.  Nitrogen  spectrum  obtained  using  a  9  cm  long  capillary  tube  with 
a  750  pm  bore.  Note  the  longer  wavelengths. 
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Figure  23,  Changing  the  N,  target  density  by  varying  the  time  delay 
slgnUicatly  altered  the  spectrum,  (a).  500  psec.  (b).  1500  psec.  (c). 

750  psec.  (d).  1000  psec. 
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Figure  24.  Changing  the  N,  target  density  by  varying  the  time  deiay 
significantly  altered  the  spectrum,  (a).  500  ^sec.  (b).  1000  psec. 

(c).  1500  (isec.  (d).  2000  psec.  (e).  3000  ^sec.  (f).  4000  psec. 
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a  new  line  of  equal  strength  while  Figures  23c-d  show  an  equally  dramatic  variation 
concerning  the  line  at  ^95  nm.  The  variation  of  the  spectral  line  intensities 
observed  In  Figure  24  is  difficult  to  explain  on  the  basis  of  only  a  changing  target 
density.  The  pulsed  valve  is  turned  on  for  approximately  250  psec.  During  this 
time  gas  is  flowing  into  the  small  reservoir  seperating  the  pulsed  valve  body  and 
the  capillary  bore.  The  pulsed  valve  closed  before  the  gas  reached  the  end  of  the 
capillary  bore  at  approximately  500  psec.  Changes  in  the  spectrum  are  expected  to 
occur  as  the  target  gas  density  increases  and  the  capillary  tube  is  filling  with  gas. 
This  change  is  clearly  evident  in  Figures  24a-b.  When  the  capillary  tube  was  filled 
with  gas  and  the  gas  density  was  maximized  throughout  this  volume,  the  spectrum 
should  exhibit  little  change  for  some  range  of  timing  change.  This  is  clearly 

shown  in  Figures  24b-c.  Beyond  this  time  delay  the  target  gas  density  is 

decreasing.  As  the  time  delay  is  Increased,  the  resulting  spectrum  should 

reproduce  the  spectrum  seen  with  a  small  time  delay.  Figures  24d-f  show  that  this 
Is  clearly  not  happening.  What  are  the  differences  in  the  target  between  Figures 

24a  and  24f  which  could  account  for  this  variation  in  observed  spectrum?  They  are 
three.  First  the  target  gas  density  really  is  different.  In  Figure  24f  the  target  gas 
density  is  the  result  of  residual  gas  remaining  in  the  capillary  tube.  In  Figure  24a 
the  gas  is  beginning  to  flow  into  the  capillary  tube  but  has  not  yet  reached  the  end 
of  the  tube.  Second,  the  resulting  lengths  of  the  target  gas  are  different  for  the 

reason  stated  above.  The  target  gas  density  integrated  along  the  interaction  length 
should  be  much  less  In  Figure  24f  than  in  Figure  24a,  because  of  the  much  lower 
gas  density.  Third,  the  target  gas  density  gradient  in  Figure  24f  should  be  small, 
since  the  gas  flow  has  ceased  whereas  in  Figure  24a  it  should  be  large  since  the 
gas  flow  is  )ust  beginning,  see  Figure  5.  This  target  gas  density  gradient  is 
probably  the  dominating  factor  in  the  observed  spectrum  because  the  Increased 
density  gradient  would  have  a  strong  effect  upon  the  propagation  of  the  laser 
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beam  through  the  capillary  tuba. 

Gas  mixtures  of  N3  and  Ne  were  also  tried.  Neon  was  added  because  Ne  Is 
the  best  harmonic  generator  found  to  date  for  KrF*.®  Mixtures  ranging  from  10% 
Ne  to  90%  Ne  were  used.  The  resulting  spectrum  was  either  the  same  or  the 
spectrum  was  weaker  depending  upon  the  spectral  region.  Variation  in  the  shot  to 
shot  laser  energy  could  have  had  an  effect,  but  it  was  difficult  to  accurately 
monitor  such  behavior.  The  spectrum  changed  with  each  laser  shot,  but  generally 
the  time  averaged  spectrum  exhibited  little  variation. 

When  the  GIS  was  employed  for  greater  resolution,  emphasis  was  place  upon 
the  third  harmonic  region.  Using  a  30  mm  long  capillary  tube  with  a  0.25  mm  bore 
produced  a  spectrum  strongly  dependent  upon  target  density,  confocal  parameter  and 
shot  to  shot  variation  in  laser  energy.  The  results  are  summarized  In  Figure  25. 
Figures  25a-b  illustrate  the  effect  produced  by  increasing  the  confocal  parameter 
and.  hence,  reducing  the  laser  energy.  Figure  25c-d  illustrate  the  large  variation  In 
signal  strength  even  when  10  laser  shot  averages  were  considered.  Increasing  the 
confocal  parameter  by  aperturing  down  the  laser  beam  had  the  effect  of  producing 
slightly  shorter  wavelengths  when  the  target  gas  density  was  increased  as  Illustrated 
in  Figure  26. 

A  quick  look  was  taken  at  the  fifth  harmonic.  Similar  structure  was  observed 
which  varied  over  the  accumalated  10  shot  averages  as  illustrated  in  Figure  27. 

The  other  target  gases  used,  H,,  He.  Ne,  Kr,  produced  similar  modulation  of 
the  third  and  fifth  harmonic  as  a  function  of  target  gas  density  and  laser  energy. 
These  results  are  summarized  in  Figures  28-31.  A  partial  spectrum  of  Xe  Is 
Included  in  Figure  32.  Little  was  done  with  either  Kr  or  Xe  because  of  the  ease  in 
ionizing  these  gases.  Ionization  effectively  prevents  the  laser  beam  propagation  in 
the  capillary  tube  as  shown  in  Figure  9.  The  lines  which  were  appearing  in  all  of 
the  gases  around  the  harmonics  were  related  to  the  fundamental  or  the  harmonics 
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Counts/Shof 


Figure  25.  Nitrogen  spectrum  near  the  third  harmonic  using  a  3  cm  long 
capillary  tube  with  a  250  pm  bore.  The  spectrum  changes  with  target 
density,  confocal  parameter  as  well  as  varied  from  shot  to  shot. 

(a).  Pulsed  valve  backing  pressure  =  20  torr.  Beam  aperatured  to 
1*  diameter,  (b).  Pulsed  valve  backing  pressure  »  50  torr.  Full 
beam  aperature.  <c).  Pulsed  valve  backing  pressure  =  100  torr.  Full 
beam.  (d).  Pulsed  valve  backing  pressure  =  100  torr.  Full  beam. 
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Figure  26.  Nitrogen  spectrum  as  in  Figure  25  with  the  laser  beam  apertured 
to  1".  The  pulsed  valve  backing  pressure  was  changed,  (a).  20  torr. 

(b).  465  torr.  (c).  700  torr. 
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Figure  27.  Nitrogen  spectrum  near  the  fifth  harmonic  showing  that  even  here 
the  spectrum  changed  as  the  target  density  was  changed  and  that  there 
varations  over  10  shot  averages,  (a).  20  torr.  (b).  50  torr. 

(c).  50  torr. 
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Figure  28.  Hydrogen  spectrum  near  the  third  harmonic  showing  the  modulation 
obtained  by  varying  the  target  density.  Each  spectrum  is  a  ten  shot 
average.  (a).  30  torr.  (b).  65  torr.  (c).  105  torr.  (d).  175  lorr 

(e).  250  torr. 
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Figure  29.  Helium  spectrum  near  the  third  and  fifth  harmonic, 
strong  variation  In  tha  single  shot  plots,  (a),  58  torr. 

(b).  750  torr,  10  shots,  (c).  1400  torr.  i  shot,  (d), 

1  shot.  (e).  850  torr,  i  shot. 
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Figure  31.  Krypton  spectrum  near  the  third  harmonic.  Single  shot  spectra, 
(a).  85  lorr.  Laser  energy  =  223  mj.  (b).  85  torr,  Laser  energy 

=  252  mJ.  (c).  375  torr.  Laser  energy  =  260  mJ.  (d).  375  torr. 

Laser  energy  =  193  mJ, 
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In  some  unknown  way.  It  is  not  clear  at  present  if  they  are  due  to  a  wave  mixing 

process  as  discussed  for  the  case  of  Ar,  although  a  six-wave  mixing  scheme  does 

exit  In  Xe. 

IV.  Photographic  Study  of  Propagation  in  Static  Gases 

In  addition  to  the  longitudinal  pumping  .studies  using  the  capillary  tube  to 
obtain  spectra,  beam  propagation  was  investigated  by  photographically  recording  the 
focal  region  in  Hj,  He  and  Air  for  different  gas  densities.  Initially  an  f/10  lens 
system  was  used  to  focus  dSO  millijoules  of  energy  In  air.  Figure  33  shows  the 
typical  results.  The  convergence  of  the  incoming  laser  beam  is  observable  as  a  blue 
glow  due  to  excitation  of  the  air.  Once  the  laser  intensity  becomes  sufficient  to 
Ionize  the  air,  the  observed  emission  becomes  white.  The  resulting  bright  fireball 
at  the  focal  zone  of  the  tens  system  Is  the  result  of  the  explosion  of  the  ionized 

gas.  At  this  point  the  laser  Intensity  is  sufficient  to  significantly  ionize  N,  and  O3 

producing  many  free  electrons.  These  electrons  should  change  the  refractive  index 
of  the  surrounding  space.  Hence,  the  resulting  plasma  should  have  either  absorbed 
the  laser  beam,  deflected  it  or  caused  it  to  diverge.  What  the  photograph  clearly 
shows  leaving  the  breakdown  region  Is  a  beam  with  a  cone  angle  less  than  that  of 
the  incoming  laser  beam.  Either  this  emerging  beam  is  the  result  of  third  harmonic 
excitation  of  the  surrounding  air,  or  some  other  process  allowed  the  KrF*  beam  to 
propagate  through  the  breakdown  region  and  emerge  with  a  smaller  divergence. 

The  results  of  focusing  the  laser  beam  In  air  lead  to  a  more  careful  study  of 
the  laser  beam  propagation  through  static  gas  targets  of  various  density.  Now  the 
KrF*  beam  was  focused  using  a  30  cm  focal  length  aspheric  mirror.  With  up  to  200 
mJ  of  energy  available  from  the  Prometheus  amplifier,  the  maximum  focal  intensity 
Is  expected  to  reach  10**  W/cm^.  It  is  assumed  that  the  laser  beam  was  focusible 
to  a  spot  size  of  2  micrometers.^ 

Figures  34-38  illustrate  some  of  the  salient  features  of  the  laser  beam 
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Figure  33.  Laser  beam  focused  in  air  using  f/iO  optics.  In  (a)  note  that 
the  beam  leaving  the  focal  breakdown  region  has  a  smaller  divergence 
than  the  Incoming  beam.  In  (b)  note  that  the  focus  region  has  two 
lobes,  the  smaller  one  appearing  just  before  the  main  breakdown. 
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Figure  34.  Laser  beam  focused  Into  Xe  at  0,59  torr  pressure,  Photo  Is 
single  shot.  Note  the  center  dark  region.  This  could  be  caused  by 
either  a  lack  of  radiating  target  material  In  the  focal  volume  or 
by  the  radiation  In  the  focal  volume  being  of  higher  energy. 
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Figure  35.  Laser  beam  focused  Into  air  at  0.15  torr,  0.35  torr,  1.23  torr 

and  3  torr.  Photos  are  of  10  to  20  shots.  The  photos  are  generally  of 
two  colors,  a  blue  or  white  center  with  long  red  tails  along  the 
laser  beam  direction.  Note  that  as  the  target  density  increases  the 
dark  centra!  region  vanishes. 


Figure  36.  Laser  beam  focused  into  He  at  0.4  torr,  1.2  torr,  3.3  torr 

and  10  torr.  Photos  are  of  10  to  30  shots.  Photos  appear  simular 
to  those  of  air  except  at  the  lowest  target  density  where  now  a 
bright  central  region  appears. 


Figure  37.  Laser  beam  focused  into  He  at  too  torr  pressure.  Photos  are 
singio  shots,  (a).  Basically  a  bright  breakdown  with  small  red 
tails,  (b).  Filter  is  used  to  block  out  the  red.  The  focusing 
lobes  are  now  clearly  visible. 
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Figure  33.  Laser  beam  focused  into  H,  at  0.36  torr.  1,1  torr,  3  torr, 

10  torr,  30  torr  and  100  torr.  Photos  range  from  40  shots  at  low 
pressure  to  1  shot  at  high  pressure.  Note  the  formation  of  the 
bright  focusing  lobes  In  (a)  to  (e)  and  the  collision  of  these  lobes 
in  (f)  as  evidenced  by  the  small  jets  perpendicular  to  the  laser 
beam  direction. 
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propagation  through  the  focal  region  observed  In  Xe,  air,  He,  and  H,.  All  of  them 
show  basically  the  same  phenomena,  a  breakdown  symmetrical  about  the  point  of 
focus.  This  is  in  contrast  to  the  asymmetric  breakdown  illustrated  in  Figure  34.  In 
all  cases,  the  laser  beam  intensity  became  sufficient  to  produce  breakdown  quite 
some  distance  before  the  focal  region.  The  figures  show  that  the  width  of  the 
breakdown  region  increases  as  the  target  gas  density  increases. 

Using  He,  the  spatial  extent  of  the  observed  blue  and  red  radiation  was 

mapped.  The  focal  region  was  mapped  onto  the  entrance  slit  of  a  0.3  meter 
spectrometer  using  a  Nikon  105  mm  lens.  The  imaging  experienced  a  reduction  of 
approximately  1.7.  Ion  lines  were  seen  to  originate  from  the  center  of  the  focal 

region  where  the  laser  intensity  would  be  greatest.  The  neutral  lines  are  spread 

out  over  the  entire  breakdown  region  and  experience  a  dip  at  the  center  where 

neutrals  would  tend  to  be  more  scarce. 

There  are  perhaps  many  factors  necessary  to  explain  the  propagation  of  a  laser 
beam  through  the  gaseous  targets  illustrated  in  Figures  33-38.  But  in  each  case  it 
is  clear  that  a  large  portion  of  the  laser  energy  is  deposited  into  the  region  leading 
up  to  the  focal  area  as  evidenced  by  the  observed  breakdown.  This  must  have  some 
effect  upon  beam  propagation  because  the  abundance  of  free  electrons  changes  the 
refractive  index  of  the  medium.  Figure  33  clearly  indicates  that  breakdown  can 
become  so  catasrophic  that  the  laser  beam  propagation  is  strongly  affected.  The 
results  of  Figures  6-9  give  an  indication  when  this  occurs  for  the  different  gases. 

V.  Conclusions 

Three  main  conclusions  are  drawn  from  the  experiments  described  above. 
First,  the  laser  beam  propagation  is  strongly  affected  by  ionization  of  the  target 
gas  and  by  any  density  gradient  in  the  target  gas  that  the  laser  beam  has  to  pass 
through  In  order  to  reach  the  focal  point.  Second,  as  evidenced  by  the  spectrum 
obtained  from  the  various  target  gases,  the  lengths  of  the  capillary  tubes  used  in 
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these  studies  were  probably  too  long  for  longitudinal  pumping  of  the  target  gases 
because  a)  once  breakdown  occurs  in  the  capillary  tube  the  capillary  tube  wall 

absorbs  the  laser  energy,  presumably  in  a  nonlinear  process,  rather  than  reflecting 
it  back  into  the  target  gas,  and  b)  the  resulting  laser  intensity  along  the  capillary 
tube  was  too  low  to  produce  sufficient  ionization  for  the  surrounding  gas  to  be 
transparent  to  the  short  wavelength  radiation  produced.  A  somewhat  greater 
intensity  is  required.  Third,  longitudinal  pumping  of  the  target  gas  in  the  capillary 
tuue  generates  a  very  strong  third  harmonic  at  82.8  nm.  This  signal  could  have 

enough  power  to  be  experimentally  quite  useful.  For  example,  using  Ne  as  the 

target  gas  could  yield  a  third  harmonic  with  as  much  as  1%  energy  conversion. 
This  would  be  a  sufficient  conversion  to  serve  many  important  application,  such  as 

holography. 
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VII.  Appendices 


Appendix  A 

Measurement  of  the  gas  target  pressure  in  the  capillary  tube  is  determined 
uc.ng  a  HeNe  laser  interferometer.  The  experimental  setup  is  illustrated  In  Figure 
A1.  The  index  of  air  at  a  pressure  of  760  torr  is  given  by 

n  -  1  =  0.000293.  A1 

Since  the  optical  path  in  the  capillary  tube  Is  given  by 

Op  =  (n  -  1)  d  P/760  A2 

where  d  is  the  length  in  the  capillary  tube  and  P  is  the  average  gas  pressure  in 
torr  Inside  the  capillary  tube,  the  number  of  fringes  for  the  HeNe  laser  is  given  by 

f  -  (n  -  1)  d  P/760X  =  0.00609  A3 

where  f  is  in  units  of  fringes  per  cm-torr.  The  average  pressure  in  the  capillary 

tube  is  given  by  combining  Equations  A2  and  A3. 

As  an  example,  the  flow  through  a  10  cm  long  capillary  tube  with  a  160  pm 
bore  was  measured  by  placing  it  inside  a  vacuum  chamber  as  illustrated  in  Figure 
Alb.  Port  B  was  opened  to  vacuum  and  the  entire  chamber  was  pulled  to  less  than 
1  mTorr  pressure.  Using  1  atm  backing  pressure  at  port  A,  gas  flow  was 
established  through  the  capillary  tube  with  port  B  open  to  vacuum.  When  port  A 
was  closed,  29  fringes  were  counted.  Likewise,  when  port  A  was  reopened,  29 
fringes  were  counted.  When  the  capillary  tube  was  not  in  position,  so  that  air 

leaked  into  the  vacuum  chamber  through  port  A,  the  fringe  shift  was  approximately 
1/2.  From  equations  A2  and  A3  this  gives  the  average  pressure  In  the  capillary 

tube  as 

Pg  =  28.5  fringes-cm-torr/(l0  cm  x  0.00609  fringes) 

Pg  =  468  torr. 

The  rate  of  flow  through  the  capillary  tube  can  be  found  by  closing  port  B  and 

opening  port  A  to  atmosphere.  When  this  is  done,  32  fringes  per  30  seconds  is 
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Figure  A1.  The  target  gas  pressure  Is  determined  by  Interferometric  means 
as  Illustrated  In  (a).  The  capillary  tube  was  positioned  Inside  a  vacuum 
chamber  as  Illustrated  in  (b)  and  gas  flowed  In  and  out  of  ports  A 
and  B  producing  a  change  In  the  fringe  pattern  observed  on  the  screen. 
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recorded. 


Since  the  volume  of  the  vacuum  chamber  Is  1.7  liters,  the  flow  rate  Is 


found  to  be 


Q  =  32  f  r  I  nges  x  760  torr  x  1. 
30  sec  296  frng 


Q  =  4.65  torr-llter/sec. 


7  liters 
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Appendix  B 


Assume  viscous  flow  through  a  capillary  tube  of  length  I  and  bore  radius  a. 
The  conductance,  Q,  through  the  capillary  tube  at  each  point  x  along  the  tube  Is 
given  by 

Q  *  ira**  P  ^  B1 

8ti  dx 

where  P  Is  the  pressure  along  the  tube  and  r|  Is  the  viscosity  of  the  gas.  The  gas 
flow  through  the  capillary  is  subject  to  several  assumptions.  First,  the  flow  speed 
is  equal  to  the  speed  of  sound,  or  the  molecular  speed,  v^.  Second,  steady  state 
flow  is  assumed  so  that 

Q  =  na^P'v^  =  constant  B2 

where  P'  is  the  pressure  at  the  output  end  of  the  capillary  tube.  The  boundary 
condition  is 

P|x=l  =  P’-  B3 

Putting  these  conditions  together  gives 

P’  -  -leqP'v^x/a*  +  Pi  B4 


where  P^  is  the  pressure  at  the  input  end  of  the  capillary  tube.  Applying  the 
boundary  condition  again  results  In  a  quadratic  equation  for  the  output  pressure  In 
terms  of  the  input  pressure  as 

2P'  =  -a  +  (a’  +  4P’)‘/*  B5 


where 


«  =  l6nVfpL/a’.  Be 

Both  equations  B5  and  B6  are  In  units  of  pbars.  The  average  pressure  in  the 
capillary  tube  can  be  obtained  by  integrating  Equation  B4  to  give 

Pa  “  (VO  /i(Po  -  Pax/ljV^’dx 

=  (2/3P'a)(Po  -  P'*J  B7 
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where  i  is  the  length  from  port  A  to  the  end  of  the  capillary  tube. 

Consider  Figure  Bl.  The  basis  structure  of  the  gas  flow  path  is  laid  out  with 
the  respective  length  and  dimensions  indicated.  For  a  given  input  pressure  of  P,  = 
760  torr,  what  is  the  output  pressure  P,  when  the  gas  is  argon?  The  molecular 
velocity  is  given  by 

=  (2kT/3m)V*  B8 

where  k  in  Bcitzman's  constant,  T  Is  the  gas  temperature  and  m  Is  the  molecular 
mass.  For  argon  at  room  temperature.  300  ®K,  the  molecular  velocity  is  3.4  x  10“ 
cm/sec.  The  viscosity  is  1.8  x  10~“  dyne-sec/cm^  For  an  effective  length  between 
Po  and  P,  of  1.55  cm  and  between  P^  and  P,  of  5  cm,  equations  65  and  B6  give 

CKi  =  74.5  torr 
Pi  =  723.7  torr 
a,  =  6840.3  torr 
Pj  =  75.7  torr. 

The  average  pressure  in  the  capillary  tube  can  be  calculated  from  Equation  B7.  It 
is  467.4  torr.  This  1$  to  be  compaired  to  the  experimental  value  found  in  Appendix 
A  of  468  torr.  From  the  difference  of  only  4%,  it  would  appear  that  the  viscous 
flow  model  of  Appendix  A  is  sufficient  to  describe  the  average  flow  through  the 
capillary  tube.  With  this  average  output  pressure  in  a  0.15  mm  bore  capillary  tube, 
the  average  target  density  wouid  be  1.7  x  10'*  cm"*. 
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Figure  Bi.  Typical  schmatic  of  the  gas  flow  path  in  the  pulsed  valve 
capillary  tube  assembly. 
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Abstract— V/t  develop  an  analytic  (heury  of  charKe-difiplacemenl  sclf- 
channclinR:  a  mcchanbm  tl,at  can  dynamically  leap  a  vhorl  intcnac 
pulic  of  HkIiI.  We  foeux  our  allenlion  on  Ihe  case  of  moii  inieresi  for 
X-ray  ampllHcalion:  llie  sIrnnRiy  saturated  channel,  for  which  all  free 
electrons  are  expelled  from  the  channel  core  and  the  channel  walls  arc 
uverdense.  S4ime  curious  results  are:  1)  the  intensity  at  the  channel 
walls  is  independent  of  the  total  laser  power,  2)  the  radius  of  Ihe  chan¬ 
nel  increases  very  slowly  with  laser  power,  asymptotically  as  the  fourth 
root,  and  3)  the  power  in  the  channel  wall  is  a  con.slanl.  The  channel 
is  cncrRelically  stable  in  the  sense  that  a  bifurcation  will  cause  a  net 
increase  in  the  electrostatic  potential  cnerRy,  but  is  only  marRinally 
stable  against  relativistic  lilamcntalion  in  Ihe  walls.  We  also  And  that 
Ihe  channel  is  an  eHective  waveguide  for  all  .secondary  radiation.  .Scal¬ 
ing  studies  show  a  substantial  advantage  exists  by  using  the  highest 
frequency  driving  laser  available.  .Since  very  large  energy-deposition 
rales  from  mullipholon  coupling  are  expected,  Ihe  channel  is  ideal  for 
generating  coherent  short-wavelength  radiation,  perhaps  well  into  the 
X-ray  range. 


I.  Introduction 

Amplification  in  the  x-ray  region  require.s  pro¬ 
digious  energy  deposition  rates  |l|  .spatially  orga¬ 
nized  in  a  high-aspcct-ratio  volume  of  material.  We  will 
show  that  the  u.se  of  extremely  intense,  shon-pulse  radia¬ 
tion  may  be  able  to  produce  the.se  conditions  by  combin¬ 
ing  the  energy  deposition  [2),  (3)  arising  from  high-ordcr 
multiphoton  proce.sses  with  a  mode  of  channeled  propa¬ 
gation  involving  the  complete  expulsion  of  free  electrons 
from  the  central  region  of  the  channel.  It  is  significant  that 
the  conditions  needed  for  the  strong  multiphoton  coupling 
arc  identical  to  those  found  required  for  the  confined 
propagation.  Previous  work  (2)  di.scus.ses  the  energy  dep¬ 
osition  rates.  Recent  work  suggests  a  paradigm  for  statc- 
•selectivc  excitation  by  the  ob.servation  of  inner-orbital 
molecular  transitions  produced  by  multiphoton  coupling 
14).  The  present  discu.ssion  concentrates  on  the  propaga¬ 
tion. 

When  a  very  intense  pul.se  of  radiation  enters  a  cold 
medium,  multiphoton  ionization  produces  a  local  plasma. 
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which  leads  to  the  following  qualitative  behavior.  The 
conversion  of  electrons  from  bound  to  free  states  by  the 
ionization  induces  a  strong  reduction  in  the  refractive  re¬ 
sponse  of  the  atomic  material.  The  resulting  plasma  ex¬ 
hibits  a  lurther  decrease  in  index  from  the  free-clectron 
component  collaterally  produced.  If  the  ambient  electron 
density  Nq  is  such  that  the  plasma  frequency  = 
sAxwoe*/ m  is  below  the  frequency  of  the  radiation  w,  an 
interesting  mode  of  channeled  propagation  could  develop. 
For  a  sufficiently  .short  pulse  (  -  100  fs).  the  ions  remain 
spatially  fixed  while  the  relatively  mobile  electrons  are 
expelled  by  the  ponderomotive  potential  from  the  high- 
intensity  region  of  the  beam.  A  state  of  equilibrium  can 
then  be  established  between  the  ponderomotive  and  the 
electrostatic  force  densities  owing  to  the  net  charge  dis¬ 
placement.  Since  the  electrons,  which  embody  a  negative 
contribution  to  the  index,  are  expelled,  sn  on-axis  region 
of  relatively  high  refractive  index  is  formed  which  can 
support  channeling  of  the  laser  beam. 

For  the  optimum  condition  conducive  to  la.ser  action, 
we  would  like  the  density  of  excited  states  to  be  as  high 
as  possible,  and  consequently,  a  high  atom  density  is  de¬ 
sired.  However,  we  may  not  have  w  <  or  the  laser 
beam  may  not  initially  propagate  to  form  the  channel. 
While  this  constraint  may  be  .somewhat  mitigated  in  the 
relativistic  regime,  in  general  it  will  be  most  propitious  to 
.set  the  atom  density  such  that  at  the  anticipated  ionization 
level,  the  plasma  frequc.ncy  is  ju.st  below  the  laser  fre¬ 
quency. 

In  this  paper,  1)  we  review  .self-channeling  in  a  low- 
density  plasma  and  present  some  simple  physical  models 
for  calculating  the  critical  intensity  for  channeling,  the 
channel  radius,  and  the  power  at  which  all  electrons  are 
expelled;  2)  we  di.scuss  the  ionization  state  achieved  in 
initially  cold  material,  and  show  that  it  is  a  good  approx¬ 
imation  to  assume  a  common  ionization  state  for  all  atoms 
involved  in  the  channel;  3)  we  show  that  the  index  of  re¬ 
fraction  is  dominated  by  the  free-electron  component;  4) 
we  develop  an  analytic  solution  for  the  strongly  .saturated 
channel  with  overdense  walls,  which  is  the  case  of  great¬ 
est  interest  for  X-ray  amplification;  5)  we  show  that  the 
motionles.s-ion  approximation  is  valid  in  the  regime  of  in¬ 
terest;  6)  we  demon.strale  energetic  stability  of  the  chan¬ 
nel,  although  it  is  only  marginally  stable  against  relativ¬ 
istic  lilumcntation;  7)  we  discu.ss  channeling  of  secondary 
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radiation  and  energy  deposition  rates;  and  8)  we  derive 
scaling  relationships  that  show  a  substantial  advantage 
with  the  use  of  high-frequency  driving  lasers. 


II.  Review  of  Self-Channelinc  in  Low-Density 
Plasma:  Simple  Analytic  Models 
Some  aspects  of  the  self-channeling  behavior  have  been 
calculated  numerically  in  the  regime  of  low-plasma  den¬ 
sity  (S).  To  get  our  arms  around  the  problem,  we  review 
those  results  and  show  how  they  can  be  obtained  from 
simple  analytic  models.  We  then  proceed  to  the  case  of 
an  initially  cold  material  subjected  to  exceedingly  high 
intensities. 

The  steady-state  force  balance  between  the  radially  out¬ 
ward  ponderomotive  force  and  the  oppositely  directed 
electron-ion  attraction  for  a  completely  ionized  tenuous 
( w  »  Uj,)  plasma  is  given  by 


2re' 

mu'c 


W(r)  =  e 


(no  -  n(r))  (r  -  r‘ ) 

- r: — ^ 


(I) 


For  a  cylindrically  symmetric  Gaussian  intensity  distri¬ 
bution  /(r),  the  electron  density  n(r)  is 


nlr)  =  no  - 


2/(0)^ 

moj'crl 


in  which  is  the  Gaussian  radius.  We  can  estimate  the 
condition  neces.sary  for  self-focusing  by  describing  the 
charge  displacement  by  two  regions  such  that  n(r  < 
To  s/2  )  =  n  <  0 )  and  n ( r  a  To  s/2  )  -  and  equat¬ 

ing  the  angle  of  total  intemiil  reflection  to  the  angle  cor¬ 
responding  to  the  first  minimum  of  dih'raction.  This  yields 
a  critical  intensity 


,  rmu'e  ,  „  .^-<1  fnu'c 

=  - - =  4.3  X  10  - W  •  cm  ' . 

64(I-Fe')r, 


(3) 

where  r,  =  e'/mc^  denotes  the  classical  electron  radius. 
Significantly,  this  intensity  is  independent  of  the  ambient 
electron  density  jiq  and  contrasts  with  the  power  threshold 
(7)  normally  arising  from  induced  index  changes  in  trunri- 
pareni  dielectrics.  We  note  that  if  the  critical  intensity  is 
not  achieved,  ionization  will  generally  result  in  u  refrac¬ 
tive  index  profile  that  will  defocus  the  beam  17).  The  abil¬ 
ity  to  focus  and  form  a  self-trapped  beam  will  al.so  depend 
on  the  initial  intensity  distribution  and  the  ri.se  lime. 

If  the  intensity  exceeds  7^.  the  beam  will  tend  to  form 
a  channel,  but  the  intensity  range  over  which  (2)  is  a  good 
description  of  the  electron  distribution  is  quite  narrow. 
The  more  common  and  more  interesting  situation  is  de¬ 
scribed  by /(O)  >  \mu'cnijrl.  Under  such  conditions,  all 
free  electrons  arc  expelled  from  the  core  ot  the  chiiimcl. 
which  is  said  to  be  saturated  or  cavitating.  The  radius  of 
the  channel  will  not  necessarily  be  the  same  as  the  laser 
focal  spot,  and  v/ill,  in  general,  depend  on  the  electron 
density. 


Here  is  a  very  simple  and  intuitive  model.  We  estimate 
of  the  radius  of  the  saturated  region  («(/■  <  r^)  =  0] 
by  treating  the  channel  as  a  dielectric  waveguide.  We  as¬ 
sume  that  the  expelled  electrons  form  a  cylindrical  sheath 
around  the  saturated  region  with  a  thickness  comparable 
to  so  that  n(r^.  <  r  <  2r„)  =  j/ic,  and  n(r  >  2r„ )  = 
no-  Consequently,  the  charge  density  n(r)  changes  from 
»o  to  Vi)  iit  the  interface.  The  cutoff  frequency  w,  of  the 
dielectric  waveguide  |8]  can  be  written  as 

1 _ m _ 

~  Ke  \4T(n(r  >  Th)  -  n(r  <  r,.))’ 

where  Xqi  =  2.405  is  the  first  root  of  the  Bessel  function 
7o(jr),  i.e.,  7o(->^oi)  =  0-  Wc  ignore  the  effect  of  the  re¬ 
gion  r  >  2r„.  The  threshold  for  propagation  occurs  at  ui 
=  for  which  we  find  /•„.  =  (jtoic  s/3/4  )/w,,  =  le  joii,. 
This  agrees  with  numerical  calculations  (5),  which  also 
.show  =  2c ! oip  for  the  unsaturated  case  and  thai  (2)  is 
u  good  description  of  the  charge  distribution  with  r,j  =  /•„ . 
The  power  P,  for  the  onset  of  saturation  can  be  found 
from  the  condition  7(0)  =  ;mu'cnoro,  which  gives 


1.74  X  10"’ 


W, 


which  is  also  in  good  agreement  with  numerical  calcula¬ 
tions  |5). 

111.  Ionization  State  and  Refractive  Properties 

The  .state  of  ionization  and  its  radial  profile  are  impor¬ 
tant  aspects  of  the  analysis.  For  these  estimates,  we  use 
the  formulation  of  Keldysh  (9)  corresponding  to  the  situ¬ 
ation  in  which  tunneling  dominates  (I0|,  |ll|.  and  have 
included  u  correction  for  the  effect  of  the  Coulomb  lield 
in  the  final  state  1 1 1 J.  The  influence  of  atomic  shell  struc¬ 
ture  is  accounted  for  through  the  use  of  computed  ion¬ 
ization  potentials  for  multiply-charged  ions  1121,  ,\n  ex¬ 
ample  involving  holmium  illustrates  the  outcome  of  this 
procedure.  It  is  found  that  for  a  pulse  of  100  fs  duration, 
nickel-like  holmium  (Ho^’’’^)  will  be  produced  for  inten¬ 
sities  spanning  the  1.3-4. 8  X  lO'"  W  •  cm  '  range 
Therefore,  if  the  spatial  beam  profile  were  a  Gaussian  with 
a  radius  and  a  peak  intensity  /„  =  4,8  x  I0‘"  W  • 
cm'*,  the  nickel  state  would  exist  from  r  =  0  to  r  = 
l.3r„,  nearly  the  entire  region  of  intere.st. 

The  radial  dependence  of  the  refractive  index  A'lrj 
governs  the  condition  for  channeling.  Since  the  ions  are 
incrtially  confined  for  the  short  time  (  -  100  fs)  consid¬ 
ered,  the  expulsion  of  the  free  electrons  from  the  channel 
results  in  a  reduced  plasma  frequency  or,  equivalently,  an 
increase  of  ihe  refractive  index  in  the  central  region. 

The  contribution  to  the  index  from  the  ions  under  the 
strong-lield  condition  involves  two  opposing  effects.  I  hcy 
are:  I)  the  nonlinear  contribution  proportional  to  ,V,. 
which  is  enhanced  at  more  elevated  field  strengths,  and 
2)  the  tendency  to  produce  high  charge  staic.s,  corre 
sponding  to  low  polarizabilitic.',  and  depressed  v.ilncs  id 
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Ni.  llie  resulting  ionic  component  of  the  index,  however, 
is  overwhelmed  by  the  free-electron  contribution.  This 
fact  is  demonstrated  by  another  example  involving  Ho  at 
an  intensity  of  3.7  x  lO”  W  ’  cm"^,  the  value  at  which 
Ho  becomes  krypcon*Iike  (Ho^''^).  The  refractive  index, 
through  the  first  nonlinear  term,  is  =  1  Nq 
where  No  =  IvN^ao  and  Nj  -  (r /Z)  NaOt^,  with  rep¬ 
resenting  the  ion  density.  We  can  relate  these  suscepti¬ 
bilities  for  Ho^’^  to  known  values  (13]  for  neutral  Kr, 
namely,  ao  =  24.8  x  lO"^*esu  and  aj  =  137.7  x  10'^* 
csu.  Since  Ho^'^  has  a  radius  of  -0.032  that  of  Kr,  hy- 
drogenic  scaling  yields  oro  —  10~^'  esu  and  =:  10~^‘ 
esu  for  Ho^'^.  At  the  given  intensity,  the  ionic  contribu¬ 
tion  then  is  {N  -!)„,=  10~^*A/a.  hor  a  collisionless 
(Ho^''^)  plasma  with  an  electron  density  of  SlA/^,  we  ob¬ 
tain  the  free-electron  contribution  to  the  index  at  248  nm 
as(/V  ~  1),  =  -(u)^/2«-)  =  -8.5  X  10 so  that 
(N  -  \)r/(N  -  1  )„,  =  -  lO’.  The  free-electron  contri¬ 
bution  is  far  larger  than  the  ionic  contribution.  Therefore, 
the  ionic  contribution  can  be  .safely  neglected. 

IV.  The  Strongly  Saturated  Channel  with 

OvERDENSE  WaLLS 

We  now  turn  to  the  case  of  greatest  interest  for  X-ray 
amplification;  the  strongly  saturated  channel.  In  Section 
II,  we  reviewed  the  case  of  a  tenuous  plasma  (ca  »  Wp). 
Here,  no  such  restriction  applies,  but  we  assume  a  density 
that  allows  enough  initial  propagation  for  the  channel  to 
form.  We  are  seeking  the  steady-state  conditions  of  prop¬ 
agation:  the  electrons  are  completely  expelled  out  to  a  ra¬ 
dius  and  the  channel  is  bounded  by  overdense  walls 
with  a  very  rapid  rise  in  electron  density.  Our  objective 
is  to  obtain  an  analytic  approximation. 

To  this  end,  we  make  the  following  assumptions;  I)  all 
electrons  are  expelled  from  the  center  of  the  channel  out 
to  a  radius  r„.;  2)  beyond  r^,  the  electron  density  rises  to 
a  overdense  value  /io».  in  a  radial  distance  that  is  very  small 
compared  to  r„,;  3)  the  density  remains  at  for  a  radial 
distance  6  representing  one  e  folding  of  the  high-fre¬ 
quency  field  and  then  drops  to  the  ambient  value  Oq  « 
now't  and  4)  the  relativistic  shift  of  the  plasma  frequency 
is  determined  by  the  high-frequency  field  at  Fh- 

A.  Analytic  Solution 

Within  the  cylinder  r  <  r^,  the  fields  for  circularly  po¬ 
larized  light  are  described  by  the  equation  [6] 

V\^p  +  Cl^l>  =  0,  (ri) 

where  is  defined  so  the  intensity  /  =  and  q  i.s  an 
eigenvalue  established  by  the  boundary  condilion.s.  The 
amplitude  of  the  high-frequency  field  will  drop  oft' expo¬ 
nentially  [14]  in  the  overdense  channel  wall  with  a  pen¬ 
etration  depth  given  by 

r  c 

b  -  (7) 

—  W"  “’/m 


where  Up„  is  the  relativistically  corrected  [15]  plasma  fre¬ 
quency, 


^pw 


VpH-O 


1  + 


rn^c'ui^ 


(8) 


and  Up^o  =  'J^rno^e~/m  is  the  nonrelativistic  plasma  fre¬ 
quency  in  the  channel  wall.  We  assume  the  same  ioniza¬ 
tion  state  throughout  the  channel  wall  as  was  Justified  in 
Section  111.  The  amplitude  of  the  high-frequency  field  for 
F  >  is  given  by 


Hr)^  (9) 


A.ssuming  that  all  the  electrons  from  the  saturated  region 
are  moved  to  a  thin  cylindrical  shell  with  inner  radius 
and  outer  radius  +  6  where  h  «  because  » 
hq,  by  using  (7)  and  (8)  we  can  approximate 


6  = 


2vnor,r^ 


(10) 


where ■  4irr,/m<j3~c. 

The  solution  to  (6),  which  represents  the  saturated  re¬ 
gion  r  <  r^,  is 


t^(r)  =  tt(0)7o(qr),  (II) 


an  expression  we  want  to  join  to  (9)  at  f  =  r„.  We  note 
that  this  treatment  is  substantially  different  from  the  treat¬ 
ment  used  by  Barnes  et  ai  [16]  which  assumed  u  Gauss¬ 
ian  solution  near  the  beam  axis. 

The  boundary  conditions  are  continuity  of  ^  and  d^/</F, 
which  give 


^(0)  7o(qFj  =  ^(fJ  (12) 


1  _  ^qI^Fh,)  2TnQF, 

•7,(^f„)  ^/T+~2^p(r^)' 


We  now  use  local  force  balance  to  obtain  the  wall  ra- 
dill.s.  The  radial  ponderomolive  force  on  an  electron  at  r, 

IS 


Fpir.) 


2ire^  d 
mw‘c  dr 

4irf  ’ 
mut'e 


^^i.0)q  Jo{qrr,)Ji(qr^). 


(  I.M 


Equation  (15)  is  nonrelativistic.  As  the  intensity  ap 
proaches  the  relativistic  regime,  the  quiver  motion  of  the 
electrons  changes  from  a  linear  displacement  in  the  dec 
tric-fidd  direction  to  u  figure  S  pattern  [17]  extending  m 
the  direction  of  propagation.  Linder  extremely  rdaiivisiK 
coiidilions,  the  ponderomotivc  force  will  be  mainly  in  ihc 
propagation  direction,  although  the  radial  gradient  ol  m 
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tensity  will  still  contribute.  We  use  (IS)  as  an  approxi¬ 
mation  to  the  force  and  we  shall  show  that  the  intensity 
at  Th,  is  independent  of  the  total  laser  power.  Thus,  the  use 
of  (IS)  is  equally  good  for  all  laser  powers. 

From  Gauss’  law,  we  obtain  the  electrostatic  force 

„  ,  ,  2  4t  T’ 

K{rw)  =  e\—\  pdp 
fw  •'O 

=  2re^nor„.  (16) 

Balancing  the  forces  =  -Ff(r^)  gives 

g  fiomu  c 

Combining  (14)  and  (17),  we  obtain  directly 

w'(o)  Ugr.)  =  --} 

and.  using  (12),  we  find 

=  Vl  +  2&yp\r^).  '19) 

The  solution  to  (19)  is  then  simply 

1  +  yfi 

<A'(0  -  — (20) 


W'(0)  Ugr^)  = 


0  20  40  60  60  too 

PyJ(O) 

Fig.  1 .  The  Bessel  function  pirameler  ai  ihe  channel  wall  as  a  function  of 
ceniral  inicnsity. 


40  60 

Pyi(O) 


The  intensity  at  the  channel  wall  Is  independent  of  laser 
power.  Furthermore, 

Wp- =  ®  0.414«,y,o,  (21) 

always  holds.  From  (12)  and  (20),  we  have 


Fig.  2.  The  channel  radius  as  a  function  of  central  intensity 


tion  of  a  laser  beam  through  a  charge-displaced  channel. 
The  channel  walls  are  overdense.  The  treatment  is  rela¬ 
tivistic  and  appropriate  for  very  high  intensity.  The  radial 
dependence  of  t!  j  laser  intensity  is 


/Vl  -F  y/2\ 

V  xlriO)^  )’ 

(22)  7(r)  =  1 

where  /({’’[•/of-*)!  =  >/o[-fo''(^)]  =  -t  and  we  have  re¬ 
stricted  the  functions  to  principal  values.  Fig.  1  shows  qr, 
as  a  function  of  the  channel  center  intensity  /(O),  as  given 
by  (22).  Notice  that  gr^  -♦  0  as  /(O)  -*  ( 1  +  and 

qr*.  -♦  Jfoi  ft.'.  /(O)  -»  00.  From  (17)  and  (22),  we  have  the 
result 


1  -F  ^  ,  /VTT^N 

r  /  Vi  +  V2\i 

''r"A7(^).' 


J\  J'S 


Fig.  2  shows  the  channel  wall  radius  r„  as  a  function  of 
the  channel  center  intensity  /(O),  as  given  by  (23).  Notice 
that  r*  -»  0  as  7(0)  -♦  ( 1  +  -v^  )//3,  as  expected,  and 
increases  indefinitely  with  7(0). 

B.  Summary 

To  summarize  results  so  far,  we  have  derived  a  solution 
to  an  approximate  model  for  the  z-independent  propaga- 


and  the  radial  dependence  of  the  electron  density  is 
'0,  r  s  Th,; 

n{r)  =  n^r^llhy  r„,  <  r  <  fy  -F  6,  I  2.‘>  I 
^/Iq.  /■  >  r„,  -F 

where  6  =  ( 1  -F  9^  )/(2xno'',r',„),  as  can  be  seen  by  %ub 
stituting  (20)  into  (10).  Fig.  3  shows  the  intensity  as  a 
function  of  radius  as  given  by  (24)  for  some  reprcsenij- 
tive  values  of  the  intensity.  Fig.  4  shows  electron  Jensny 
di.stributions  as  given  by  (25)  together  with  the  coire 
.sponding  intensity  distributions. 

C.  Asymptotic  Forms 

A  Taylor  expansion  at  jcoi  gives  Jq(x)  ==  y|(.r„i )  ' 

-  X),  which  can  be  used  to  find  the  asymptotic  bchavnT 
of  (22)  as  7(0)  -»  oo. 


gr^  -  A'oi  - 


'l  +  v/2 


y,(xoi)V^ 


=  -TOI  ~ 


s/^7(0) 


=  I  Hfi  I 
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f,V5iinjr, 


Fig.  3.  The  ndial  disiribulion  of  intensity  for  various  values  of  (he  central 
intensity. 


Fig.  5.  Comparison  of  the  asymptotic  solution  obtained  by  Taylor  expan¬ 
sion,  (27),  to  the  exact  solution,  (23),  for  the  channel  radius  as  a  function 
of  central  intensity.  Dashed  curve  is  the  asymptotic  solution. 


(b) 

Fig.  4.  The  electron  density  distribution  superimposed  on  the  density  dis¬ 
tribution  for  various  values  of  central  intensity,  (a)  •  20.  (b) 

^^'(0)  -  100. 


A  Taylor  expansion  also  gi\esJo(x)J,(j;)  =  /kxoiXxoi 
~  jt),  which  can  be  combined  with  (26)  to  find  the  asymp¬ 
totic  behavior  of  (23),  specifically, 


Wi(-tQ,)  V(1  +  V2)/3/(0)  ~  1  -  y/2 

2xnor, 

1.94Vd/(^  -  2.41 


(27) 


Fig.  5  compares  the  asymptotic  solution  given  in  (27)  to 
the  exact  solution  given  in  (23).  Remarkably,  the  asymp¬ 
totic  solution  agrees  with  the  exact  solution  within  a  few 
percent  over  the  entire  domain  of  definition:  it  is  an  ex¬ 
cellent  approximation.  Equation  (27)  shows  that  in  the 


limit  of  very  high  power,  the  channel  radius  increases  as 
the  fourth  root  of  intensity  at  the  channel  center  or  about 
as  the  fourth  root  of  total  laser  power. 

V.  Validity  of  the  Motionless-Ion  Approximation 

The  motionless-ion  approximation  is  quite  good  over 
the  intensity  and  pulse-length  regime  of  interest.  The  pon- 
deromotive  force  on  an  ion  is  proportional  to  its  charge 
and  inversely  proportional  to  its  mass  M.  For  Ho  at  Iq  ~ 
4.8  X  I0^°W  *  cm'^,  the  pondeiomotive  force  moves  an 
ion  a  distance  -4  x  10' times  that  of  the  electron, 
which  is  completely  negligible.  A  far  larger  force  on  the 
ions  results  from  the  charge  displacement  itself  since  the 
ions’  mutual  electrostatic  repulsion  will  expel  them  from 
the  beam.  However,  the  ion  motion  due  to  this  force  is 
also  small  for  a  sufficiently  short  pulse.  Assume  that  the 
saturated  channel  has  been  formed  at  /  -  0.  With  the  elec¬ 
trostatic  force  on  an  ion  at  radius  r  <  given  as  F  = 
2xr(l^e)^na  Gauss’  law,  an  ion  with  initial  position  rg 
has  its  motion  described  by  r  =  rg  exp  (  y^)  where 

which  leads  to  the  approximate  result  for  small  displace¬ 
ments: 

^  -  '•q  ^  ^2  ^29) 

fg  M 

At  /g  =  4.8  X  10^°  W  •  cm"^,  the  ion  migration  gives  a 
relative  change  in  its  position  of  about  9  percent  after  lOU 
fs. 

VI.  Energy  Loss  Mechanisms 

The  two  principal  mechanisms  leading  to  energy  loss 
in  the  channel  are:  1)  ionization,  and  2)  the  field  energy 
associated  with  the  ch.^rge  displacement. 

A.  hniiation 

The  energy  lost  to  ionization  can  be  found  by  summing 
the  ionization  potentials  [12]  of  electrons  successively 
stripped  from  the  atom  up  to  the  ionization  level  produced 
by  the  field.  For  the  Ho”^  case,  with  0^  =  2.3  x  lO" 
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cm'^  and  a  channel  radius  ~  M***.  io**-  while  in  the  region  outside 

ization  produces  a  loss  of  -~4.3  mJ  •  cm"‘.  I  +  V2  )  (““ 

B.  Charge  Displacement  0  Jr. 

The  radially  directed  electric  field  is  given  by  j  /  ■  +  J2)  [  r 


exp  [2(r,  -  ryb]rdr 


£(r)  =  2ire«o 


fw  +  r^{r^  -  r)/5, 

/•„  <  r  <  r^,  +  6; 
0,  r  2:  +  5 


where  we  have  corrected  for  the  approximation  in  (25). 
By  direct  integration,  the  charge-displacement  field  en¬ 
ergy  per  unit  length  of  the  channel  is  given  by 

^  _  1  f*  £2/;.^^;.  -  +  6)  (3r,„  +  5) 

flh  4  Jo  *  ^  ^  ~  12 

2  2  2  4  2 

t  €  rtorj.  me 


•  Ijro,y,(xo,)  V(T  +  yfl)  0/(0)  -  1  -  V^f. 

(31) 

For  the  example  above,  this  would  be  -44  ■  cm~'. 

VII.  Stability  of  the  Channel 
Two  questions  concerning  the  stability  of  such  channels 
have  been  examined.  The  first  involves  energetic  stability 
a.ssociated  with  the  splitting  of  a  single  channel  into  two 
or  more  channels,  each  having  a  fraction  of  the  total  elec¬ 
trostatic  energy.  The  second  considers  stability  agaiast  re¬ 
lativistic  self-iWu.sing  in  the  walls  of  the  channel,  a  mode 
observed  in  numerical  simulations  [18]. 

A.  Energetic  Stability 

For  appraisal  of  the  case  based  on  energetics,  using  the 
approximation 

^^jI(x')x' dx'  =Xoii?(-roi)  1  -  ^1 


which  is  accurate  to  within  10  percent  over  the  0  <  x  < 
Xfli  range,  for  a  total  power  P  =  Fq  4  F, ,  we  find  that 
inside  the  saturated  region 


Po  =  2t/(0)  Jl{qr)rdr 

^  C  J*M0l\O) 
v2/lftr-fxni*/ 1  f  JTni  )  \  /3/( 


_  VI  ^  >/2 

2nar,[x0iJx  (xqi  )  >^/(0)  -  J 

r  r  V’y 


2ir(  1  +  V2  )  r  ^ 

0  L  2  4 

_  mu)^c(  1  4-  V2  )^  { 1  +  +  4jrnor,ri;) 

Slx^nlr^rl 

•«> 

Combining  (31)  and  (33)  and  taking  the  limit  /(O)  -♦  00, 
we  find  further 


1  [(1  +  y/l)  rxoiVr^eC'noP  ]'^' 


Xoi7|(xoi)  \0/{b) 


a  result  which  shows  that  the  electrostatic  energy  e,  per 
unit  length  /  scales  as  (laser  power)^'^\  Therefore,  this 
energy  would  be  increased  by  splitting  into  two  channels, 
each  of  which  contained  half  the  laser  power.  Further  in¬ 
spection  of  (31)  and  (33)  shows  that  an  even  more  stable 
situation  exists  for  powers  less  than  the  asymptotic  limit. 

B.  Stability  Against  Relativistic  Filameniation 

The  toul  power  propagating  in  the  channel  walls  is 
given  by  (34).  The  power  threshold  for  relativistic  seif- 
focusing  (19),  [20]  is 

4T/Jog 

in  which  A"  =  1.5  x  10'°  W,  so  that  the  condition  P,  > 
Pi  barring  relativistic  .self-focusing  in  the  walls  is 

£  >  ^(1  4-  v^)’  =  2.6  X  10'®  W, 

Although  this  stability  criterion  is  not  met,  h  is  unly 
slightly  surpassed.  The  assumption  that  the  relaiivisiic 
shift  of  the  plasma  frequency  is  given  by  the  inicn.siiy  at 
the  channel  wall  slightly  overestimates  the  total  power 
flowing  in  the  walls.  Thus,  we  can  only  conclude  that  the 
stability  against  relativistic  filamentation  is  marginal 

Furthermore,  although  the  power  in  the  wall  ot  the 
strongly  saturated  channel  is  always  the  same,  ihe  time 
(or  distance)  for  the  instability  to  develop  fully  depends 
on  the  geometry  of  the  channel.  A  thick  annulus  will  til- 
ament  more  quickly  than  a  thin  annulus.  The  more  power 
flowing  down  the  channel,  the  longer  it  will  take  lor  (he 
instability  to  develop.  The  question  of  rclativistu  nla- 
mentation  can  be  addressed  by  detailed  numerical  simu¬ 
lations,  but  ultimately  only  by  experiment  [21]  Both  arc 
beyond  the  scope  of  this  paper. 

C.  Linear  Stability 

If  the  channel  is  to  be  used  as  the  excited  vniurnc  'nr 
an  X-ray  laser,  it  is  desirable  to  have  the  channel  pfi’i'j 
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gate  in  a  straight  line.  Linearity  of  the  channel  can  be 
ensured  by  laying  down  a  track  of  higher  refractive  index 
that  the  pulse  will  follow.  For  a  gas,  this  might  be  do.ne 
with  a  laser  prepulse.  In  a  solid,  microfabrication  tech¬ 
niques  might  be  used. 

VIll.  Tov.'ard  X-Ray  Amplification 

Several  additional  considerations  pertain  to  the  use  of  a 
charge-displaced  channel  for  obtaining  X-ray  amplifica- 


A.  Channeling  of  Secondary  Radiation 

Using  the  dielectric  waveguide  as  a  model  for  the  chan¬ 
nel  reveals  an  interesting  aspect  of  this  mode  of  propa¬ 
gation.  We  see  that  if  the  frequency  w  exceeds  the  cutoff 
in  (4),  a  condition  essential  for  the  channel  to  form,  then 
all  higher  frequencies  will  also  be  guided.  This  occurs 
because  the  index  is  dominated  by  the  free  electrons  or, 
equivalently,  the  angle  of  total  internal  reflection  is  pro¬ 
portional  to  the  angle  of  the  first  minimum  of  diffraction. 
Since  all  of  the  radiation  which  falls  within  the  acceptance 
angle  of  the  waveguide  will  propagate  as  a  guided  mode, 
these  channels  can  serve  naturally  as  bright  directed 
sources  of  energetic  radiation. 

We  note  that  linear  stability  is  only  a  practical  require¬ 
ment;  we  need  to  know  which  way  the  laser  is  poir^ting. 
The  channel  acts  as  an  optical  fiber  for  all  secondary  ra¬ 
diation. 

B.  Deposition  Rates 

Energy  deposition  rates  associated  with  such  channels 
are  expected  to  be  extremely  high.  Extrapolation  (21  of 
experimental  energy-transfer  rates  sugge.sts  that  the  mul¬ 
tiphoton  absorption  cross  section  will  limit  to  a  fre¬ 
quency-independent  value  of  -  10*‘“  cm^  for  heavy  ele¬ 
ments  at  intensities  above  lO'^  W  •  cm"^,  a  value  near 
the  Compton  intensity  in  the  ultraviolet.  For  the  holmium 
case  cited  above,  with  an  electron  density  /Iq  ~  8.9  x 
10'“^  cm~^,  the  rate  is  in  exce,ss  of  10'“*  W  ■  cm’’,  a  value 
far  above  that  needed  for  amplification  of  radiation  in  the 
kilovolt  region  (I). 

C.  Population  of  Active  Atoms 

From  (27),  we  find  that  the  number  of  atoms  per  uni. 
length  in  the  saturated  core  of  the  channel  is 

WiUoi)  Vo  -F  -  1  -  >/2 


where  f  is  the  degree  of  ionization.  Using  248  nm  radia¬ 
tion  at  4.8  X  10™  W  ■  cnr’.  we  ti.id  ^  1  x  lu'" 

cm  " '  or,  recalling  the  example  of  hoimium,  n,^^^  1  ,H  x 
lO"  cm"'.  This  is  a  rather  small  linear  density,  although 
volume  density  is  quite  high.  It  v  ill  be  necessary  to  u.se 
the  atoms  very  efficiently.  Perhaps  inner-shell  inultipho- 
ton  ionization  is  the  test  mechanism. 


D.  Scaling  with  Driving-Laser  Frequency 

The  deposition  rate  has  a  strong  dependence  on  primaiy 
laser  frequency,  which  gives  a  substantial  advantage  to 
the  use  of  higher  frequencies. 

When  a  tenuous-plasma  channel  is  un.saturated  or  only 
slightly  saturated,  its  radius  is  ~2c/<i>p.  The  power  prop¬ 
agating  down  the  channel  is  roughly  proportional  to  its 
central  intensity  times  its  square  of  the  radius  (!ij.  To 
compare  the  effectiveness  of  different  driving  lasers  for 
creating  channels  that  might  become  lasers  themselves,  it 
is  appropriate  to  compare  systems  with  the  same  values 
of  w/wp.  So  fora  given  central  intensity,  the  driving-laser 
power  scales  as  to"^. 

For  the  strongly  saturated  channel,  which  is  the  main 
.subject  of  this  paper,  driving-laser  power  scales  even  more 
rapidly  with  frequency.  From  (33),  we  obtain 

_  Vgf'(Q)(i  +  '^) 

^  "  ^0  -  2no'-,JCoi 

^  4i?(xo,)  /tVV/^'(0)7i  +  >/2) 

in  the  high-power  limit.  Forcnannels  with  the  same 
the  driving-laser  power  scales  as  Furthermore,  highei 
frequencies  generally  allow  shorter  transform-limited 
pulse  lengths,  enabling  the  ions  to  remain  motionless  at 
higher  driving-laser  power.  The  ion  migration  .scales  ap¬ 
proximately  as  the  square  of  the  pulse  length.  If  pulse 
length  then  scales  as  w"',  the  migration  .scales  as  u  ’. 
Finally,  the  charge-displacement  energy  given  in  (31)  be¬ 
comes 

^  -ro|y)(.roi)  (1  4  yfi}  131(0) me' 
dx  I6r, 

x4iy?U,)(l  +  V^)/3/(0)c 


in  the  high-power  limit.  So  b.e  charge-displacement  en¬ 
ergy  .scales  as  w'^,  while  the  corresponding  iooization  en 
ergy  goes  asoj”^. 

IX.  C0NCLU.S10NS 

The  analysis  we  have  presented  shows  thi/t  nonlinear 
absorption  anc  channeled  propagation  combine  coopera¬ 
tively  in  producing  conditions  favcrable  to  the  u.se  o!  in¬ 
tense  short-pul.se  high-frequency  radiation  for  the  produc 
tier  of  X-ray  amplification.  NumericaS  calculations  will 
be  rtsjUired  to  obtain  a  detailed  de.scriptian  of  the  dynum 
ic.s  of  channel  formation  and  .stability;  liiowever,  the  pres¬ 
ent  technology  i;:  clcse  to  being  able  to  answer  rnunv 
questions  in  the  laboratory  121J.  The  ability  to  channel 
both  the  energy  deposition  and  resulting  emission  can 
enormously  increase  the  high-frequency  gein  and  provuli- 
a  natural  mecha  iism  for  generating  low-divergence  high 
brightnc.ss  X-ray  .sources. 
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The  characterittics  of  relativiatic  propagation  in  plaamas  of  lubpicosecond  ultraviolet  (248  nm)  radia¬ 
tion  arc  studied  for  both  spatially  homogeneous  plasmas  and  plasma  columns.  It  is  established  for  the 
first  time  that  the  defocusing  properties  of  the  interaction  can  represent  a  dynamical  mechanism  stabiliz¬ 
ing  the  mode  of  propagation  against  radial  oscillations.  The  dependence  of  both  quasistabilized  modes 
and  pulsing  waveguide  regimes  on  the  initial  transverse  intensity  distribution  is  examined  and,  for  the 
latter,  the  locus  of  the  fint  focus  produced  in  a  homt^eneous  plasma  is  calculated. 


PACS  numben;  S2.40.Ob.  42.10. -s.  42.65.Jt 

The  dynamics  of  propagation  of  extremely  intense 
subpicosecond  pulses  of  radiation  in  self-generated  plas¬ 
mas  is  a  rapidly  developing  area  of  study.  In  this  Letter 
certain  aspects  of  this  question  are  treated  theoretically 
for  the  following  ranges  of  physical  parameters:  peak  in¬ 
tensity  /as  I0"-10^*  W/cm^  pulse  duration  raO.S-l.O 
ps,  initial  radial  aperture  roSal-B  >tm,  wavelength  X 
-248  nm,  and  plasma  elearon  density  — 10^-10^' 
cm  A  key  finding  of  these  calculations  is  insight  into 
the  physical  processes  enabling  the  formation  at  a 
quasistable  mode  of  propagation  of  such  high-intensity 
radiation  in  plasmas. 

For  pulses  having  a  peak  intensity  in  the  10'^*I0^'- 
W/cm^  range,  rapid  (Arai  i  ft)  processes  of  multiphoton 
ionization''^  occurring  on  the  .leading  edge 
W/cm^)  of  the  wave  form  will  remove  several  (~6-l0) 
atomic  electrons.  Thereby,  for  a  focused  beam  with  a 
given  radial  aperture  ro,  a  plasma  column  is  rapidly 
formed  in  which  the  central  temporal  zone  of  the  pulse 
propagates.  The  spatial  and  temporal  dynamics  govern¬ 
ing  the  resulting  propagation  are  determined  by  ihe  com¬ 
petitive  interaction  of  diffraction,  defocusing,  and  self- 
focusing.  We  note  that  in  the  early  portion  of  the  pulse, 
some  defocusing  is  expected  to  occur  while  the  ionization 
is  commencing,  since  the  contribution  to  the  refractive 
index  from  the  free  electrons  tends  to  locally  reduce  iu 


value  in  the  central  high-intensity  region.  In  addition,  as 
shown  below,  the  process  of  defocusing  can  give  rise  to  a 
previously  unknown  mechanism  of  stabilized  propaga¬ 
tion. 

For  subpicosecond  pulses  and  the  radial  dimensions 
considered,  grou  plasma  motion  involving  the  ions  is 
negligible.^  Similarly,  for  full  ionization,  noninertial 
Kerr  self-focusing^  ia  insignificant.  In  the  present  case, 
both  the  relativistic  increase  in  the  electron  mau  arising 
from  induced  oscillations  in  the  intense  field  and  the 
charge  displacement  resulting  from  expulsion  of  free 
electrons  from  the  ionized  column  by  the  ponderomotive 
force can  produce  strong  self-focuaing  action. 

We  now  describe  the  specific  equations  governing  the 
propagation  and  present  calculated  solutions  illustrating 
the  development  of  self-focuaing  and  channeled  propaga¬ 
tion.  The  origin  of  the  stability  exhibited  by  the  propa¬ 
gating  radiation  is  of  primary  impcrtance.  In  consider¬ 
ing  the  central  temporal  region  of  the  pulse,  we  neglect 
the  energy  loases  associated  with  further  Ionization  end 
spatial  broadening  of  the  plums  c  si.-cr  these 

ePecU  can  be  quite  smu^L  '  Wm  i  ,  .>•  condition,  the 
S|/av*  bchs-'.-r  the  i/tdnviio.i  field  in  the  column, 
including  relativistic-'Chsrge-displacement  self-focuaing. 
diffraction,  and  defocusing,  is  governed  by  the  nonlinear 
j  Sebrodinger  equation 


Cl  «/  dz 


2ko  2»k^  2«/», 


(I) 


where  E(i,z,r)  is  the  slowly  varying  (l,z)  complex  amplitude  of  the  wnvc,  r  is  the  transverse  coordinate. 
Aj.’"Wr9r)rd/dr,  4o“2irA,  ci  is  the  group  velocity  in  the  plasma,  and  c  is  the  vacuum  speed  of  light. 

The  noiilincr  term  Sep{r,\E\^)  describing  the  combined  relativistic'charge-displacement  lelf-focusing  and  defocut- 
ing  of  the  pulse  is  expressed  as 


"Sitr  «R“i -«//(«*+ v^).  SR.-l . 

—  {4xe^N,l/]/Irw,(l  +  3//4)  (»|,“4wr^yV,.o/rj„ 


/,-Wa’cV4xf^  /-(c/2x)|£|*. 
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In  thac  expressions,  Uf  is  the  relativistically  shifted 
pUisnia  frequency,  oi  is  tite  laser  angular  frequency,  v  is 
the  electroP'iun  ooUiskmal  frequency,  !,  is  the  relativistic 
intensity/'^  and  Nt\J\  is  the  electron  density  arising 
frciii  the  quasistatic  balance  of  the  ponderomotive  and 
the  electrostatic  forces  densities.^  We  assume  that 
Ni{r,z)^Ntfif(.r,z)  is  the  static  ion  density,  in  which 
Ni.o  is  a  constant  denoting  the  maximal  ion  density  and 
/(r.z)  describes  the  spatial  distribution  of  the  ions. 
Hence, /(r,z)>0,  max/(r)*“l  I/(r,z)«l  for  homo¬ 
geneous  plasma].  For  quasineutrality,  this  gives 
A,,o/(/',z)  as  the  initial  electron  density  with  Ntja 
and  Z  the  ionic  charge. 

The  expreuion  for  A,[/},  which  describes  the  pertur¬ 
bation  of  the  plasma  from  local  neutrality  by  both  the 
ponderomotive  and  relativistic  mechanism,  is  given  by 

/V,I/l"Af,,omax  0,/(r,z)+»fdiv  -  )  1/2 

where 

K  “  (2ca>^m,,yV,,o)  “ ' ,  <3) 

and  m,  denotes  the  electron  mass.  In  deriving  Eq.  (2), 
the  electronic  pressure  has  been  neglected,  since  it  is 
small  for  the  beam  and  plasma  paramcten  examined 
here.  The  expression  maxiO,— )  simply  provides  for 

Calculations  of  the  propagation  have  been  performed 
corresponding  to  the  following  parameters  for  the  radia¬ 
tion  field  and  the  plasma;  11*248  nm,  //*1.34xlO” 
W/cm^  /o-f/,ai3.0xlO”  W/cm^  ro-3  nm,  N1.0 
-7.5x10” cm"’,  and Z- 10, 

Normalization  of  r  and  /  by  the  initial  beam  radius  ro 
and  the  maximal  value  of  the  initial  intensity  /o,  respec¬ 
tively,  gives  «f| —(cVeT’*  1.4x10“’,  the  dimensionless 
analog  of  the  coefficient  k  defined  by  Eq.  (3).  Therefore, 
for  this  set  of  conditions,  the  charge  displacement  ’  is  ex¬ 
pected  to  be  relatively  small  to  that  its  influence  is 
confined  to  the  structure  of  the  focal  regions  which  arise 
from  the  relativistic  mechanism.  Because  of  this,  we  will 
utilize  the  approximation  that  x— 0  in  the  current  work. 
It  it  important  to  note,  however,  that  other  choices  of  the 
parameters  can  lead  to  a  generally  strong  role  for  the 
charge-displacement  proccu  and  future  study  is  being 
directed  to  that  case. 

The  collitiona!  absorption  of  light  in  plasma  is  given 
by  */  —  -  V—  — p  ~/ko,  where  /i  ”  is  the  absorp¬ 

tion  coefficient.  Estimates’  of  n  ~,  for  the  values  of  the 
parameters  used  in  this  work,  show  that  the  collisional 
absorption  is  extremely  weak  and  docs  not  influence  the 
numerically  obtained  solutions.  Therefore,  the  results 
presented  ^low  have  been  calculated  with  n  ~  —0. 

The  nonlinear  Schrodinger  equation  stated  in  Eq.  (I) 
is  only  applicable  to  the  description  of  the  dynamics  of 
the  field  in  a  small  focal  zone,  if  its  longitudinal  scale  / 
satisfies  the  condition  The  numerical  experiments 


showed,  for  the  range  of  parameters  studied,  that  / 
>  50X.  In  this  connection,  we  note  that  the  spatial  limi- 
utions  on  the  minimal  dimensions  of  the  focal  zones,  and 
consequently  the  bounds  on  the  intensities  developed  in 
those  regions,  arise  from  saturation  of  the  nonlinearity 
expressed  by  6eitU,\E\^). 

The  dynamics  of  propagation  vith  relativistic  self- 
focusing  have  been  examined  usinf  model  pulses  i(z,r,t) 
having  spatial  and  temporal  Gaussiai,  or  hyper-Gaussian 
incident  beam  profiles  given  by 

/(0,r,»)»/o(r’,l)  — /*,  exp{  —  O/t)^'  —  (r/ro)  ,  (4) 

NiZ:2,  Ni>2 

for  /*—/,»/• —10'*  W/cm’,  the  latter  being  the 
threshold  intensity  above  which  rapid  multiphoton  ion¬ 
ization  occurs  to  the  particles  in  the  medium.  Further¬ 
more.  we  have  concentrated  on  the  central  temporal  tone 
of  the  pulse,  |/|  ^  lo,  where  to  satisfies  the  expression 

/o</o)r/o(0./o)-/i«exp(-0o/r)^')»/* ,  (5) 

Ihe  transverse  profile  of  the  plasma  column  created  uy 
the  front  of  the  pulse  is  simulated  by  the  hyper-Gauuian 
functiun 

/(r)  — exp{  -  (r/r^ )  ^*1 ,  2: 2 ,  ( 6 ) 

with  the  aperture  of  the  plasma  column  r,  defined  by 

/o(r*,/o)«yo(fo)wpl~(r'*//’o)^’)  — /*  .  ‘'I 

For  a  Oauttian  transverse  intensity  distribution  (,V . 
an  aperture  defined  with  /*— lO”  W/cm’,  and  ukmg 
we  find  r«ai2,68ro.  Consequently,  ihe 
homogeneous  pluma  approximation  /(r)*  I  is  valid 

The  beam  profiles  of  the  solutions  exhibit  vrverii 
salient  features.  Among  them  are  the  producuon  of  mul¬ 
tiple  foci  and,  under  appropriate  conditions,  the  develop¬ 
ment  of  a  quasistable  confined  mode  of  propagation  lack¬ 
ing  the  generation  of  sharply  focused  regions  Imtsvr 
tantly  for  the  latter,  the  calculations  reveal  the  nature  A 
the  mechanism  leading  to  the  stabilization.  The  re^ult» 
of  computations  presented  in  this  work  are  given  m  iiw 
frame  G,z,r),  j-r-fz/ci),  connected  with  the  •«*« 
front  of  the  beam.  The  data  in  the  figures  below  iiiu« 
tratc  the  propagation  of  the  radiation  along  the  :  ito  '.>r 
j  —  const. 

Figure  1(a)  illustrates  the  propagation  that  develop* 
for  the  physical  parameters  stated  above  correspi>mling 
to  a  pulse  with  an  entrance  intensity  loU)  —  ]  having 
a  transverse  distribution  given  by  /V}-8.  In  accord  wuh 
an  earlier  study,*  as  the  nonlinearity  dcgflfp)  1*  va 
turated  at  a  sufficiently  low  level,  the  behavior  of  put «es 
propagating  in  a  homogeneous  plasma  without  abworp- 
tion  can  generate  a  mode  sim'lar  to  a  putting 
waveguide.*'”  The  result  shown  in  Fig.  1(a)  illutiriirv 
Ihe  formation  of  such  a  regime.  Interestingly,  the  pu  1 
ing  waveguide  represents  an  nlternatio.-i  of  ring  mii, 
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FIG.  I.  The  (r,x)  diitributioM  for  reUlivutic  lelMocuaing 
of  the  propaBatini  inUtuily  for  «  peek  initiel  intensity  It  •"it, 
hevinf  a  trensverM  distribHlion  given  by  in  (4)  for 
the  beam  and  piuma  parameters  defined  in  the  teat,  (a)  Spa> 
tially  homogeneous  fuliy  stripped  plasma,  (b)  Axially  lym* 
metric  plasma  column  with  r«**r«  whose  transverse  profile 
corresponds  to  N)  *8  in  Eq.  (6). 

tures  and  focal  regions  along  the  axis  of  propagation. 
Significantly,  the  power  trapped  in  it  corveiponda  to  ap- 
proximately  90%  of  the  total  initial  power  of  the  beam. 

From  F/q.  (7),  it  follows  that,  for  placeaulike  initial  in* 
tensity  distributions,  the  aperture  r«  of  the  plasma 
column  becomes  comparable  with  the  beam  aperture 
For  example,  in  the  case  treated  above,  we  have 
r»  » l.28ro.  As  the  apertures  of  the  laser  beam  and  the 
plasma  column  tend  toward  coincidence,  however,  dc" 
focusing  becomes  relatively  more  significant.  Therefore, 
this  aspect  of  the  propagation  must  be  carefully  taken 
into  account  when  the  evolution  of  pulses  having  plateau' 
like  initial  intensity  distributions  was  studied. 

Figure  I  (b)  shows  the  propagation  of  a  pulse  in  a  plat* 
ma  column  for  Ni  ‘~8  corresponding  to  the  same  physi¬ 
cal  parameters  as  those  pertaining  to  the  illustration  in 
Fig.  1  (a).  In  this  case,  the  transverse  profile  of  the  elec¬ 
tron  density  in  the  plasma  column  is  described  by/(r)  in 
Eq,  (6)  for  Nj  “8  and  r,  —ro,  The  comparison  of  Figs. 
1  (a)  and  1  (b)  demonstrates  the  strong  effect  of  defocut- 
ing  on  the  spatial  dynamica  of  propagation  when  the 
column  hai  a  radial  dimension  close  to  thst  of  the  beam. 
Dcfocusing  causes  a  fraction  of  the  beam  to  spread  away 
from  the  column  while  the  remaining  energy  adjusta  to 
balance  the  relativistic  lelf-focusing,  dcfocusing,  and 


FIG.  2.  The  locus  (i,z)  of  the  first  focal  zone  for  a  puUe 
with  a  Gauiaian  initial  spatial  and  temporal  intensity  distribu¬ 
tion  for  the  parsmeters  defined  in  the  text.  The  solution  per- 
uins  u>  relativisiic  Mlf-foctuing  in  a  spatially  bomogencoui 
fuliy  stripped  pUsma. 

diffractioa.  In  thia  way,  the  dynamica  of  the  interaction 
tend  toward  the  formation  of  a  spatially  stable  mode  of 
propagation.  In  addition,  as  evi^nt  from  Fig.  Kb),  a 
substantial  fraction,  nominally  23%  of  total  initial  beam 
power,  can  be  confined  in  such  a  mode.  Importantly,  the 
further  study  of  such  cases  has  also  indicated  that  such 
stable  patterns!  can  be  maintained  for  a  length  of  at  least 
twenty  normal  diffraction  lengths. 

The  solutions  of  Eq.  (I),  with  the  initial  conditions 
given  by  Eq.  (4),  satisfy  important  similarity  laws. 
Specifically,  the  results  shown  in  Figs.  1  (a)  and  1  (b)  are 
valid  for  any  other  set  of  physical  parameters  »atisfyina 
the  relationships  a|Ss((Aroro)VsRolw^«/<»^~'2.486x lo' 
and  *”2/9.  Furthermore,  if  full  account  of  both 

the  relativistic  and  charge-displacement  mechanisms 
(xisuO)  is  made,  these  similarity  statements  are  pre¬ 
served  nwf  r?  •»oj/4ai. 

Following  earlier  studies'*’’"  we  have  analyzed  the 
evolution  of  plane  waves  having  small  perturbations  in 
order  to  determine  tlie  dependence  of  the  self-focusing 
length  on  the  initial  pulse  parameten.  The  results  of  this 
analysis  showed  that  the  Largest  growth  rate  of  the  per¬ 
turbations  and,  consequently,  the  minimum  seif-focusing 
length  occur  at  an  entrance  intensity  on  the  beam  axis  of 
IqU)"'}!,.  Specifically,  the  conipiitatiotu  have  shown 
that  this  inference  remains  true  for  the  relativistic  self- 
focusing  of  pulses  with  Gaussian  initial  intensity  distri¬ 
butions. 

The  dynamic  motion  of  the  foci  illustrated  in  Fig.  2 
exhibits  important  characteristics.  For  example,  the  cal¬ 
culations  show  that  the  locus  of  the  first  focus  in  the 
(f,r)  plane.^  for  the  initial  condition  detemiined  by  Eo. 
(4)  with  yV,-2.  yV2“2.  7* -6L as 8.04x10*’  W/cmi 
r*”500  fs,  and  ro**3  pm,  reaches  the  rriinimal  z  for 
tgO)  "•  ft, .  If  the  extreme  intensity  on  the  beam  axis  is 
Imi-tr.  this  locus  hss  three  reversal  points.  Two  of 
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them,  comspooding  to  the  sanw  z,  are  due  to 
toil)  1 1,.  The  third  and  oentrai  one,  corresponding  to 
a  greater  z,  occurs  at  In  principle,  this 

feature  establishes  a  clear  diagnostic  signature  for  rela¬ 
tivistic  self-focusing  that  enables  it  to  be  distinguished 
from  the  Kerr  nonlinearity*  which  produces  a  locus  with 
a  tingle  point  of  reversal. 

The  spatial  character  of  propagation  of  sufficiently  in¬ 
tense  subpioosecond  ultraviolet  pulses  in  self-generated 
plasmas  it  significantly  influenced  by  a  relativistic- 
charge-displacement  mechanism  and  two  important 
classes  of  propagation  have  been  distinguished.  In  a  spa¬ 
tially  homogeneous  plasma,  the  interaction  can  produce 
self-focusing  in  the  f^orm  of  a  pulsing  waveguide  involv¬ 
ing  a  structure  of  alternating  rings  and  focal  zones.  The 
location,  formation,  and  nature  of  the  focal  regions,  in¬ 
cluding  the  maximal  intensities  occurring  in  them,  de¬ 
pend  essentially  on  the  properties  of  the  initial  transverse 
intensity  distribution.  The  locus  of  the  first  focus  in  the 
(/,z)  plane  for  pulses  with  Gaussian  initial  intensity  dis¬ 
tributions  and  maximum  intensities  has  three  re¬ 
versal  points.  For  this  class  of  pulses,  the  relativistic 
self-focusing  length  is  minimal,  if  the  intensity  of  the  ini¬ 
tial  pulse  on  the  axis  equals  f//. 

A  mechanism  leading  to  a  spatially  stabilized  mode  of 
propagation  is  also  revealed.  In  this  case,  the  defocusing 
that  occurs  when  the  plasma  column  has  a  radius  close 
to  that  of  the  beam  is  particularly  significant.  The 
effective  action  of  the  defocusing,  which  serves  as  a  spa¬ 
tially  distributed  nonuniform  energy  loss,  causes  the  rap¬ 
id  formation  of  intensity  profiles  that  become  quasistabi- 
lized  along  the  beam  axis.  Therefore,  it  is  accurate  to 
state  that  the  process  of  defocusing  contributes  directly 
to  the  dynamical  development  of  stabilization.  Further 


computational  resulu  using  the  procedures  described 
thovt  are  available.'^ 
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ABSTRACT 

Calculations  of  the  dynamics  of  propagation  in  plasma  of  intense  axisymmetric 
laser  pulses  incorporating  both  relativistic  and  charge-displacement  mechanisms 
are  presented.  It  is  shown  that  the  combined  action  of  these  two  processes  can 
lead  to  stabilized  confined  modes  of  propagation  whose  asymptotic  properties 
correspond  well  to  the  lowest  z-independent  eigenmodes  of  the  governing  nonlinear 
Schrodinger  equation.  The  effect  of  the  charge-displacement  is  large  and  cavitation 
of  the  electron  distribution  is  a  general  feature  of  the  solutions  for  the  range  of 
parameters  studied.  The  results  of  a  specific  example  show  that  approximately 
one  half  of  the  incident  power  (~  4  TW)  can  be  trapped  in  the  channels  which, 
for  a  wavelength  of  248  nm,  an  electronic  plasma  density  of  ~  7,5  x  10^“  cm'*, 
and  an  initial  beam  radius  of  -  3  pm,  develop  propagating  intensities  of  ~  10^‘  W/cm*. 
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The  dynamics  of  propagation  of  intense  subpicosecond  laser  pulses  in  plasmas 
in  the  strongly  relativistic  regime  is  currently  undergoing  vigorous  theoretical 
analysis^"®.  Of  primary  significance  is  the  possible  formation  of  quasistable  self- 

0  Q 

channeled  modes  of  propagation.  Separately,  the  analyses  of  both  relativistic 
and  charge-displacement^  mechanisms  give  indications  that  confined  modes  of 
propagation  can  exist.  In  this  Letter  we  examine  theoretically  the  consequences  of 
the  combined  action  of  these  two  mechanisms  with  particular  attention  to  the 
development  and  character  of  highly  stable  confined  modes  of  propagation.  For 
radiation  with  a  wavelength  of  248  nm,  the  range  of  intensity  of  interest  is  I  > 
10‘*  W/cm^  since  both  relativistic®  and  charge-displacement^  mechanisms  are 
significant  in  that  region. 

The  basic  formalism  and  procedures  describing  these  calculations  have  been 
described  previously®.  Assuming  collisionless  and  lossless  propagation®  (p  =  v  = 

*  0),  as  in  the  earlier  analysis,  the  complex  field  amplitude  E(t,2,r)  is  governed 
by  the  nonlinear  Schrodinger  equation 

1  3  3  \  I  ik, 

- +  —  p  ^  -  A  E  t-  -  depE  »  0,  (1) 

Cj  9t  9z  J  2kj) 

where  r  is  the  transverse  coordinate,  Aj.  =  3V3r*  +  r'*9/3r,  kj,  =  2ti/X,  c.  = 

Is  the  group  velocity  in  the  plasma,  and  c  is  the  vacuum  speed  of  light.  The 

nonlinear  term  i5ep(r,  |E|*)  of  key  significance,  which  describes  the  combined 

action  of  the  relativistic  and  charge-displacement  mechanisms,  is  given  by 

1  -  (1  +  3I/Ir)  -‘/"max  C,  f{r)  +  (cVwp^)Aa.(1  +  3l/Ir)‘/’jj. 

(2) 

in  which  1^  ~  3mg6)’cV4it8’  is  the  relativisiic  intensity®  ®,  is  the  plasma  frequency, 
u  is  the  laser  angular  frequency,  and  f(r)  denotes  the  spatial  distribution  o(  the 
initially  unperturbed  electron  density.  The  last  term  in  the  inner  brackets  of  Eq, 
(2),  (cV^'ipg)Ax(l  +  3I/If)'/^,  describes  the  charge-displacement  process  and  was 
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neglected  in  the  earlier  studies^  of  purely  relativistic  propagation.  The  form  of 
Eq.  (2)  assumes  that  the  propagating  radiation  is  circularly  polarized,  since  the 
relativistic  T-factor  can  be  written  as  (1  +  31/1^)'/®  in  that  case.  Futhermore, 
the  structure  of  the  terms  in  Eq.  (2)  clearly  reveals  that  the  relativistic  and 
charge-displacement  processes  are  fundamentally  connected  and  cannot  be  considered 
as  truly  separate  and  independent  mechanisms. 

In  consideration  of  the  charge-displacement,  the  motion  of  the  heavy  ions  is 
neglected,  in  accord  with  previous  treatments^,  since  the  time  scale  of  the  laser 
pulse  is  assumed  to  be  sufficiently  short  (~  100  fs).  In  this  picture,  it  is  known 
that  the  electronic  component  of  the  plasma  can  be  substantially  perturbed  by  the 
ponderomotive  potential  generated  by  the  propagating  radiation For  the  present 
calculation,  the  electron  distribution  is  determined  by  the  force  balance  between 
the  radially  outward  ponderomotive  force  and  the  oppositely  directed  electron-ion 
attraction^.  This  ponderomotively  driven  process  is  capable  of  completely  expelling 
electrons  from  certain  spatial  regions,  a  condition  described  as  electronic  cavitation"^ 
The  calculations  show  that  the  charge-displacement  plays  a  strong  role  in  determining 
the  spatial  character  of  ti.e  propagating  energy  and,  particularly,  in  stabilizing  the 
confined  high-intensity  modes. 

The  essential  finding  of  these  calculations  is  insight  into  the  formation  and 
stabilization  of  spatially  confined  modes  of  propagation  arising  from  the  cooperative 
action  of  the  relativistic  and  charge-displacement  processes.  The  numerical 
simulations  were  performed  with  model  pulses  I(z.r,t)  having  spatial  and  temporal 
Gaussian  or  hyper-Gaussian  incident  beam  profiles  and  initial  homogeneous  or 
hyper-Gaussian  unperturbed  electron  densities,  in  accord  with  the  procedures 
described  previously®.  Using  that  earlier  notation®,  the  relevant  parameters  are  X 
248  nm,  If  =  1,34  X  10’°  W/crn’,  !«  =  3  x  10‘’  W/cm’,  r^  =  3  pm,  Ng  q  = 


Ul 


7.5  X 


10*®  cm'*.  Na  =  2.  and  f(r)  ■  1  (spatially  uniform  plasma).  Moreover,  the  results 
described  below  are  valid®  for  any  other  set  of  physical  parameters  satisfying  the 

(Tgiilp  \  2 

- ^1  =  2.486  X  10*  and  a,  ■  lo/I^  =  2/9.  The  results  of 

computations  presented  in  this  work  are  given  In  the  frame  (q.z.r),  q  =  t  -(z/c,), 
connected  with  the  wave  front  of  the  beam.  The  graphical  data  represented  below 
illustrate  the  propagation  of  the  radiation  along  the  z  axis  for  q  =  constant. 

The  solutions  for  the  spatial  profiles  of  the  electromagnetic  energy  and  the 
corresponding  electron  density  exhibit  several  important  characteristics  which  are 
shown  In  Fig.  (1)  and  Fig.  (2),  respectively.  Among  them  are  (1)  the  generation 
of  intense  focal  regions,  (2)  the  development  of  a  stable  confined  channel  of 
propagation,  (3)  the  trapping  of  a  substantial  fraction  of  the  incident  power  in 

the  channel,  and  (4)  strong  cavitation  in  the  electron  density. 

Figure  (la)  Illustrates  the  overall  behavior  of  the  normalized  intensity  distribution 
I(r.2)/Io  and  Fig.  (lb)  magnifies  the  region  between  the  first  two  peaks  appearing 
in  Fig.  (la).  The  corresponding  normalized  electron  density  N(r,2)/No  is  shown  in 
Fig.  (2a).  Initially,  although  a  small  amount  of  energy  can  be  seen  diffracting 

away,  a  large  fraction  of  the  incident  radiation  collapses  to  form  both  an  axial 
focal  zone  (-45%)  and  a  high-intensity  annular  region.  This  early  stage  of  evolution 
involves  very  little  charge-displacement  and  resembles  rather  closely  the  dynamics 
seen  In  the  purely  relativistic  case®.  However,  in  the  present  situation,  the 
development  of  these  high-intensity  zones  leads  rapidly  to  strong  perturbations  in 
the  electron  distribution  and  the  formatiort  of  two  regions  of  cavitation.  This 
modified  charge  distribution,  the  cross  section  of  which  Is  shown  In  Fig  (2b)  for 
the  axial  position  z  =  95.4  pm,  clearly  illustrates  simultaneous  cavitation  along  the 
axis  (r  <  0.25  pm)  and  in  a  relatively  narrow  annular  region  located  at  r  -  1.7  pm. 
Subsequently,  as  energy  is  exchanged  between  these  two  structures,  the  annulus 
transfers  some  of  its  energy  to  the  axial  region  and  a  smooth  stable  confined 
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mode  of  propagation  on  the  z-axis  containing  ~  46%  of  the  incident  power  develops 
(z  ?  600  pm).  The  electron  density  evolves  in  a  corresponding  manner  and  a 
single  cavitated  zone  on  the  axis  is  formed  (z  S  600  pm).  The  peak  intensity 
associated  with  the  channeled  propagation  is  exceptionally  high,  reaching  values  of 
~  10*‘  W/cm®  in  this  example. 

The  relativistic  influence  on  the  electron  mass  and  the  ponderomotively 
driven  displacement  of  charge  both  encourage  the  formation  of  channeled  propagation. 
On  the  basis  of  the  dynamical  picture  revealed  above,  the  combined  effect  of 
these  two  mechanisms  can  be,  in  reasonable  approximation,  summarized  in  the 
following  simple  way.  The  relativistic  effect  leads  to  the  initial  concentration  of 
the  radiation  and  the  resulting  displacement  of  electronic  charge  reinforces  this 
tendency  and  stabilizes  the  confinement.  The  cooperative  nature  of  this  action 
appears  to  lead  to  highly  stable  conditions  of  propagation. 

The  behavior  of  the  computed  solutions  for  large  z  can  be  compared  with 
the  z-independent  eigenmodes^-^^  of  the  nonlinear  Schrodinger  equation  (1).  In 
this  comparison,  it  has  been  shown  that  the  field  amplitude  tends  asymptotically 
to  the  lowest  elgonmode  Vg(r)  of  this  equation.  V^{t)  is  a  real-valued  positive 
monotonIc  function  of  r,  carries  a  real  dimensionless  index  s,  vanishes  as  r 
and  is  related  to  the  electric  field  amplitude  by  the  relation 

e(r.z)  =  Vs(r)  expl-iz(l  -  s)wjy2k^,c"].  (3) 

The  computations  have  shown,  for  the  case  studied,  that  the  field  distribution 
for  large  z  tends  to  the  amplitude  \/5(r)  lor  the  index  s  =  0,554.  The  appropriate 
intensity  distribution  Is  then  given  by  I^ft)  ■  Vg*(r},  Figure  (2c)  exhibits  the 
normalized  asymptotic  field  amplitude  (Is(r)/Io)*'^’  and  the  corresponding  asymptotic 
electron  density  distribution  Ng(r)/Ng.  Significantly,  these  asymptotic  profiles, 

v/hen  compared  to  the  corresponding  distributions  arising  from  tlie  numerical  dynamical 
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calculation  of  the  propagation  for  z  ^  900  pm,  differ  by  less  than  one  percent. 
Therefore,  the  limiting  behavior  of  the  dynamically  derived  propagating  solutions 
seems  well  described  by  the  simpler  z-independent  analysis. 

Earlier  work^  concerning  the  charge-displacement  mechanism  with  a  steady- 
state  picture  showed  that  the  electrostatic  energy  per  unit  length  associated  with 
a  channel  arising  from  the  charge-displacement  scaled  as  P*/®,  with  P  denoting  the 
laser  power.  Consequently,  for  a  fixed  total  power,  this  electrostatic  energy  will 
necessarily  be  a  minimum  for  a  single  channel.  Therefore,  the  collapse  of  the  annular 
structure  shown  In  Fig.  (1)  is  consistent  with  the  system  dynamically  seeking  to 
minimize  this  electrostatic  energy.  Indeed,  for  the  electron  density  distributions 
shown  In  Fig.  (2),  direct  calculation  shows  that  the  electrostatic  energy  per  unit 
length  of  the  large  z  solution  (asymptotic)  is  approximately  one  half  that  associated 
with  the  region  containing  the  cavitated  annulus  [Fig.  (2b)].  In  addition,  it 
should  be  noted  that  the  electrostatic  energy  per  unit  length  in  this  example  is  far 
less  (<  0.2%)  than  the  corresponding  energy  per  unit  length  in  the  radiation  field. 

The  results  of  the  dynamical  calculations  shown  in  Fig,  (1)  and  Fig.  (2)  can 
be  directly  compared  in  three  ways  to  the  simple  analytical  picture  previously 
used  to  estimate  the  properties  of  the  propagation  arising  solely  from  the  mechanism 
of  charge-displacement^.  (1)  The  analytical  treatment  found  that  the  intensity  at 
the  wall  was  independent  of  laser  power  and  equal  to  lyy  =  (1  +  /2)mai’c/4irrg 
=  1.1  X  10®®  W/cm®.  The  dynamical  result  presently  calculated  gives  good  agreement 
with  lyy  ~  1.2  X  10®®  W/cm®.  (2)  The  peak  intensity  developed  at  the  center  of 
the  channel  In  the  asymptotic  regime  is  also  found  to  agree  within  -  10%.  (o) 
Finally,  the  channel  radius  shown  in  Fig.  (2c)  is  ~  0.4  pm,  while  the  analytic 
treatment^  gives  a  corresponding  value  of  ~  0.7  pm.  This  difference  is  at  least 
partially  understood  by  the  fact  that  the  analytic  theory^  did  not  take  into 
account  the  relativistic  increase  in  electron  mass  in  the  equation  for  the  force 
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balance.  The  use  of  a  higher  mass  would  naturally  reduce  the  effective  ponderomotive 
force,  a  correction  that  would  translate  into  a  reduced  channel  radius  as  indicated 
by  the  relativistic  caiculation.  Overaii,  however,  we  can  state  that  the  fuil 
numerical  treatment  of  the  dynamics  and  the  approximate  analytical  theory  produce 
results  that  are  in  rather  close  agreement  for  these  important  properties  of  the 
asymptotic  behavior. 

Several  other  cases  have  also  been  considered.  In  particular,  calculations 
examining  the  dynamics  of  propagation  have  been  studied  for  initial  intensity 
distributions  other  than  Gaussian®  (e.g.  Nj  =  8)  for  both  homogeneous  piasmas 
(f(r)  ■  1]  and  simuiated  plasma  columns  If(r)  =  exp  (-(r/ro)'].  The  behavior  was 
found  to  be  qualitatively  the  same  as  that  shown  in  Fig.  (l)  with  asymptotically 
trapped  fractions  of  the  initial  laser  power  ranging  from  34%  to  77%. 

The  dynamical  behavior  of  channeled  propagation  arising  from  both  relativistic 
and  charge-displacement  mechanisms  has  been  studied  theoretically.  Two  principal 
findings  have  been  established.  (1)  These  two  processes  can  cooperatively  reinforce 
one  another  and  lead  to  stable  confined  high-intensity  modes  of  propagation  in 
plasmas  which  are  capable  of  trapping  a  substantial  fraction  (~  50%)  of  the  incident 
power,  (2)  The  spatial  profiles  of  both  the  intensities  and  electron  distributions 
tend  asymptotically  to  those  derived  as  lowest  eigenmodes  of  a  z-independent 
analysis. 
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Fig.  (1);  (a)  The  normalized  intensity  distribution  l(r,z)/I(,  calculated  for  an 


initially  homogeneous  plasma  [f(r)  “1).  lo  =  3  x  10‘*  W/cm“,  r^,  =  3 
pm,  X  =  248  nm,  and  =  7.5  x  W/cm’.  (b)  Enlargement  of  the 
transition  zone  between  the  first  two  peaks  shown  in  Fig.  (la). 

Fig.  (2):  (a)  The  normalized  electron  density  distribution  N(r,z)/N|,  calculated  for 

an  initially  homogeneous  plasma  (f(r)  «  1]  and  corresponding  to  the 
data  shown  in  Fig.  (1).  N,,  =»  Nqq  =  7.5  x  lO’"  cm'*,  (b)  Radial  cross 

section  of  the  normalized  electron  density  shown  in  Fig.  (2a)  for  axial 
distance  z  =  95.4  pm.  Cavitation  is  present  along  the  central  axis  (r  2 
0.25  pm)  and  in  an  annulus  located  at  r  ~  1.7  pm.  (c)  Radial  dependence 
of  the  asymptotic  solutions  for  the  normalized  amplitude  (is(r)/lt,j‘/’ 
and  the  normalized  electron  density  Ng(r)/tio  for  the  index  s  =  0.554. 
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ABSTRACT 

Experimental  studies  examining  a  new  relativistic  regime  of  high-intensily 
short-pulse  propagation  in  plasmas  have  been  performed  which  present  evidence  for 
the  formation  of  a  stable  mode  of  spatially  confined  (channeled)  propagation.  For 
an  electron  density  of  ~  1,35  x  10^‘  cm’’  and  a  power  of  ~  3  x  10“  W,  the 
re-sult.s  indicate  a  channel  radius  <  i  pm  and  a  peak  intensity  -  10^’  W/cm’. 
Comparison  of  these  findings  with  a  dynamical  theory  yield  agreement  for  both 
the  longitudinal  structure  and  the  radial  extent  of  the  propagation  observed. 
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A  fundamentally  new  regime  of  electronnagnetic  propagation  is  expected  to 

arise  in  plasmas  for  short-pulse  radiation  at  sufficiently  high  intensity.  Calculations 

of  the  propagation  in  plasmas,  incorporating  both  relativistic^ and  charge- 

3—7 

displacement  mechanisms,  indicate  that  the  combined  action  of  these  processes 
can  lead  to  a  new  stable  form  of  spatially  confined  (channeled)  propagation. 
This  Letter  reports  (1)  the  results  of  the  first  experimental  study  probing  the  physical 
regime  relevant  to  the  observation  of  relativistic/charge-displacement  self-channeling 
and  (2)  presents  the  initial  comparison  of  these  experimental  findings  with  matching 
theoretical  calculations  performed  with  the  computational  procedures  described  in 
Ref.  7. 

The  experimental  arrangement  used  in  these  studies  is  illustrated  in  Fig.  i(a). 
The  source  of  radiation  was  a  subpicosecond  KrF*  (X  =  248  nm)  laser  that  has 
been  described  elsewhere.®  It  delivered  a  linearly  polarized  power  of  -  3  x  10“ 
W  (-  150  mj,  pulse  duration  ~  500  fs)  in  a  beam  with  a  diameter  of  ~  42  mm. 
When  this  radiation  was  focused  into  the  chamber  with  lens  LI  (f/7),  a  focal 
radius  Tq  -  3.5  pm  was  measured,  giving  a  maximum  intensity  Iq  ~  8.6  x  10‘^ 
W/cm^  The  medium  was  provided  by  filling  the  chamber  statically  with  gas  (He, 
Ne,  Ar,  Kr,  Xe,  Nj,  COj,  or  a  mixture  of  Xe  (4%)  and  (96%)]  up  to  a  maximum 
density  of  -  1.89  x  10’®  cm'^. 

The  diffracted  248  nm  radiation  was  measured  as  a  function  of  the  angle  (0) 
with  respect  to  the  direction  of  the  incident  radiation.  The  incident  laser  beam 
was  blocked  by  a  metal  disk  on  the  output  window  of  the  chamber  and  lens  L2 
imaged  the  region  near  the  focal  zone  on  a  fluorescent  screen  S.  The  diaphragm 
D  in  front  of  lens  L2  restricted  the  collection  of  the  diffracted  light  to  a  solid 
angle  of  -  5®  while  simultaneously  increasing  the  depth  of  field.  The  angle 
between  the  axis  of  the  lens  L2  and  the  axis  of  the  incoming  laser  radiation 
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could  be  readily  varied  up  to  a  maximum  arigle  of  0  -  15*.  Two  flat  mirrors  coated 

for  high  reflection  99%)  at  248  nm,  both  having  a  spectral  bandwidth  of  -  10 

nm,  served  in  reflection  as  spectral  filters  (F)  for  the  diffracted  laser  radiation 
so  that  only  the  scattered  248  nm  radiation  could  illuminate  the  screen.  An 

attenuator  A  was  employed  to  adjust  the  intensity  on  the  screen  and  the  images 
formed  were  recorded  through  the  visible  fluorescence  produced  with  a  microscope 
and  a  CCD  camera. 

The  characteristic  behavior  observed  is  well  illustrated  by  the  data  recorded 

with  Na.  The  measured  result,  shown  in  Fig.  (ib),  corresponds  to  a  density  ' 

1.35  X  10®°  cm'®.  To  the  left,  in  the  photographic  inset,  a  relatively  large  cone 
of  light  Rayleigh  scattered  from  the  plasma  is  visible  at  all  angles  as  the  energy 
propagates  toward  the  focal  point  of  the  lens,  while,  in  the  region  to  the  right  of 
the  conical  apex,  a  narrow  filament  developed.  The  diameter  of  this  filament  is 
not  greater  than  10  pm,  the  measured  spatial  resolution  of  the  imaging  system 
The  distribution  of  intensity  observed  along  the  filament  exhibited  several  bright 

features  attributed  to  diffraction  because  they  could  not  be  seen  for  0  >  20“ 

Since  the  axis  of  the  imaging  lens  corresponded  to  an  angle  0  =  7.5®,  the  scale 
along  the  abscissa  of  the  photographic  data  is  reduced  by  -  8-fold,  giving  the 
maximum  length  of  the  filament  as  ~  1  mm.  The  graph  in  Fig.  (1b)  represents  a 
one-dimensional  axial  profile,  taken  along  the  direction  of  propagation  (z),  of  the 
observed  intensity  pattern  (inset).  Three  peaks  (a,  p.  T)  are  visible  with  a  spatial 

separation  of  6  =  200  ±  20  pm.  The  normal  Rayleigh  range  for  the  focal  geometry 
used  was  ~  200  pm. 

The  diameter  of  the  filamentary  channel  is  an  important  dynamical  variable 

which  we  estimated  by  measuring  the  maximum  angular  deviation  of  the  diffracted 
light.  The  experimental  value  (f)  of  this  diffracted  cone  was  4>  ~  20®,  a  magnitude 
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indicating  a  radius  -  0.9  pm  though  the  relation  41  =  1.22  X/r^.  Filaments  of 
this  general  nature  were  observed  at  densities  above  -  1.35  x  io”  cm'*  in  Nj, 
Ne,  Ar,  Kr.  COj  and  a  mixture  of  Xe  (4%)  and  Nj  (96%).  but  not  in  He  and  Xe, 
two  materials  discussed  further  below. 

Two  mechanisms  exist  that  could  modify  the  refractive  index  of  the  medium 

and  lead  to  the  observed  behavior.  They  are  (1)  the  Kerr-effect  stemming  from 

the  ions  and  (2)  the  relativistic/charge-displacement  process.^  Since  the  pulse 

duration  is  very  short  (-  50C  fs),  the  motion  of  the  ions  is  negligible,^  and  no 

contribution  can  arise  from  expulsion  of  the  plasma  from  the  high-intensity  zone. 

An  implication  of  the  estimate  of  the  channel  radius  (r^  ~  0.9  pm)  is  that  the 

observed  propagation  is  associated  with  intensities  in  the  10‘*~10^*  W/cm*  range. 

9  10 

Under  tnese  conditions,  available  experimental  evidence  '  on  muliiphoton  ionization 

indicates  that  Ho  should  be  fully  ionized  and  the  C,  N,  and  0  atoms  constituting 

the  molecular  materials  would  retain,  at  most,  only  is  electrons. 

9  10 

Consider  explicitly  the  case  of  N,.  which  has  estimated  ’  threshold  intensities 
for  the  production  of  N®*.  N**,  and  of  1.6  x  10‘*  W/cm*,  6.4  x  10‘’  VV/cm*. 
and  1.3  x  10^*  W/cm^,  respectively.  Hence,  the  volume  of  the  channel  would  be 
largely  ionized  to  N**,  with  certain  localized  high-intensity  regions  contributing 
some  N**.  Two  consequences  of  this  pattern  of  ionization  follow,  namely,  (ot)  the 
Kerr-effect  arising  from  the  ions  is  small,  since  the  polarizabilities  of  the  remaining 
Is  electrons  are  low,  and  (P)  the  electron  density  (ng)  initially  produced  in  the 
focal  region  is  nearly  uniform.  Therefore,  ng  =  1.35  x  10*'  cm'*  for  the  data  on 
Nj  shown  in  Fig.  1(b). 

A  critical  power  for  self-channeling,  arising  from  the  relativislic/charge-- 
displacement  mechanism,  can  be  defined^’  as 
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Per  *  (m|cV0*)  f  go(Orclr  (u/wp  q)’  =  1.62xl0‘*  ]  W.  (1) 

0  '  ng/ 

where  mg,  e,  and  c  have  their  customary  identifications,  u  is  the  laser  angular 
frequency,  Up  q  is  the  plasma  frequency  for  the  uniform  unperturbed  plasma  with 
electron  density  ng,  is  the  critical  electron  density  (  for  X  =  248  nm,  n^^  ^ 
1.82  X  10”  cm'*),  and  gQ(r)  is  the  Townes  mode.^^ 

The  critical  powers  associated  with  the  experimental  conditions,  for  He  and 
Nj  at  a  medium  density  p  =  1.35  x  10*®  cm'*,  are  1,08  x  10‘*  W  and  2.19  x  lO" 
W,  respectively.  Therefore,  since  the  incident  power  was  P  =  3  x  10^‘  W,  no 
filament  was  expected  in  He,  a  prediction  conforming  with  the  observation  of 
none.  Moreover,  the  diffracted  cone  of  radiation  was  also  atsent  with  He.  In 

contrast,  P/Per  =  f  for  N,,  a  condition  that  held  generally  (P/Per  >  1)  for  all 
materials  which  exhibited  evidence  for  channel  formation.  We  note,  however,  that 
some  contribution  from  the  Kerr-effect  may  be  present,  even  for  the  light  materials 
(Ne,  Nj,  and  CO,),  in  the  early  stage  of  channel  formation  prior  to  the  development 
of  a  substantial  level  of  ionization,  and  that  the  heavier  gases,  (Ar,  Kr,  and  Xe) 
may  involve  a  more  significant  influence  from  the  Kerr  process.  A  specific 

estimate  of  the  nonlinear  index  change  arising  from  both  N**  and  N**  at  an 
intensity  of  -  10‘*  W/cm*  indicates  that  their  contribution  is  less  than  10'*  that 
of  the  free  electrons,  hence  the  ionic  contribution  can  be  neglected  in  N,  for  the 

conditions  studied. 

A  direct  comparison  will  now  be  made  between  the  theoretical  analysis,  fully 
described  in  Ref,  (7),  and  the  experimental  findings  for  N,.  This  comparison  can 

be  accomplished  for  both  the  longitudinal  intensity  profile  and  the  radial  extent 

of  the  channel.  Figure  2(a)  illustrates  the  intensity  profile  I(r,z)/Io  calculated 
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with  physical  parameters  corresponding  to  those  of  the  experiment  (i.e.  P  =  3  x 
10‘*  W.  Tq  »  3.5  pm,  ng  =  1.35  X  cm'*,  and  P/Pgr  =  1.37).  Importantly,  all 
of  these  parameters  are  based  on  independent  measurements  of  (1)  the  laser  pulse 
involving  determinations  of  the  energy  and  power  (P),  (2)  the  focal  radius  (Fq)  of 
the  incident  radiation,  and  (3)  the  characteristics  of  the  multiphoton  ionization^ ° 
generating  the  electron  density  (ng).  Therefore,  this  comparison  does  ngt  involve 
a  fit  to  a  free  parameter.  The  normalized  electron  density  calculated  is  presented 
in  Fig.  2(b)  from  which  it  is  seen  that  electronic  cavitation  occurs  only  near  the 
positions  of  the  maxima  in  the  intensity  profile  (Fig.  2(a)].  Curve  (A)  in  Fig.  2(cj 
repre'-.ants  the  one-dimensional  axial  intensity  profile  I(0,z)/Iq  corresponding  to 
the  calculated  distribution  shown  in  Fig.  2(a).  The  spacing  (6)  of  the  maxima  is 
seen  to  be  6  -  185  pm,  a  value  in  close  agreement  with  the  experimental  figure  (5 
=5  200  ±  20  pm)  illustrated  in  Fig.  1(b).  Furthermore,  analysis  has  shown  that  the 
spacing  <S  is  quite  sensitive  to  the  power  P^  and  electron  density  ng,  particularly 
in  the  region  close  to  the  threshold  (see  Eq.  (1)].  With  respect  to  the  results 

illustrated  in  Fig.  2(c),  an  increase  in  Og  by  less  than  5  percent  causes  a  reduction 
in  the  spacing  (5  by  approximately  25  percent.  Therefore,  substantial  changes  in  tne 
physical  parameters  would  grossly  alter  the  comparison  of  the  experimental  ana 

theoretical  results. 

3—7  1 1 

Theoretical  studies  ’  indicate  that  the  charge-displacement  plays  a  ver, 
important  dynamical  role.  In  order  to  test  this  hypothesis,  identical  calculations 

were  made  for  for  the  purely  relativistic  case^  which  explicitly  neglects  the  charge- 
displacement  term,  namely,  elimination  of  the  term  (cV(“p,o'’o))  *  I/Ir)"  ‘ 

in  Eq.  (24)  of  Ref.  (7).  Significantly,  the  resulting  axial  profile  [Curve  B  in  Fig 
2(c)]  exhibits  only  a  single  relatively  weak  maximum,  for  0  <  z  4  600  pm,  an  outcome 

sharply  at  variance  with  both  the  full  theoretical  analysis  and  the  experimental 
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observation.  Although  the  expected  charge-displacement  is  highly  localized  [Fig. 

3  7  11 

2(b)l,  this  comparison  reveals  the  strong  influence  it  has  on  the  propagation. 

At  a  greater  incident  power  (P/Pcr  ^  10).  a  continuous  channel  in  the  electron 
distribution  is  expected  to  develop.^ 

The  measurements  indicated  an  approximate  value  of  r^  =  0.9  pm  for  the 
radial  extent  of  the  channel,  a  result  that  can  be  compared  with  the  corresponding 
theoretical  figure.  Figure  2(d)  illustrates  five  radial  intensity  profiles  I(r,Zj  )/Io 
of  the  distribution  pictured  in  Fig.  2(a).  Since  the  measurement  of  this  angularly 
scattered  radiation  did  not  correspond  to  a  known  longitudinal  position,  this  comparison 
can  only  be  qualitative,  but  the  radial  distributions  shown  indicate  that  the  expected 
value  lies  in  the  interval  0.5  <  r  <  1.0  pm,  a  range  that  comfortably  includes  the 
experimental  value  r^. 

The  results  observed  with  Xe  deserve  additional  discussion,  since  those 
experiments  did  not  give  evidence  for  the  formation  of  a  channel,  Jn  significant 

contrast  to  the  case  involving  N,.  the  electron  density  ng  produced  by  the  multiphoton 
ionization^®  in  Xe  is  expected  to  be  very  nonuniform  spatially.  For  intensities 

spanning  10“  -  10“  W/cm’,  the  corresponding  density  would  vary  by  over  a 
factor  of  two.  Since  this  nonuniformity  woulo  tend  to  reduce  the  refractive  index 
locally  in  the  central  high-intensity  region,  a  significant  defocusing  action  is 

expected  which  could  suppress  the  channel  formation. 

Finally,  we  note  (1)  that  the  intensity  distribution  is  not  expected  to  depeno 

13  14 

Strongly  on  the  state  of  polarization  and  (2)  that  tosses  to  the  plasma  may  be 
significant,  particularly  at  electron  derisities  close  to  n^,f. 

The  first  experiments  examining  a  new  relativistic  regime  of  high-intensiiy 
pulse  propagation  in  plasmas  have  been  performed  and  the  findings  indicate  the 
formation  of  a  channeled  mode  of  propagation  over  a  length  considerably  greater 
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than  the  Rayleigh  range.  Specific  comparisons  of  the  experimental  observations 

with  a  dynamical  theory,  which  expiicitiy  includes  both  the  influence  of  the 

relativistic  mass  shift  and  the  displacement  of  the  electronic  component  of  the 

plasma,  produce  excellent  agreement  for  both  the  longitudinal  structure  of  the 

intensity  profile  and  the  radial  extent  of  the  channel.  While  the  present  channel 
contains  several  foci,  a  continuous  channel  is  predicted  to  develop  at  higher 
power.  Finally,  the  intrinsically  very  high  concentration  of  power  associated  with 
this  mechanism  of  channeled  propagation  provides  an  efficient  and  general  method 
for  the  production  of  conditions  necessary  for  x-ray  amplification.^® 
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Figure  Captions 

f  ^y,.  1:  (a)  Experimental  apparatus  use  in  studies  of  propagation.  See  text  for  i 

description.  (b)  Data  concerning  the  pattern  of  propagation  observed 
with  a  single  pulse  in  N,  at  a  density  of  ~  1.35  x  10^°  cm'*.  The  maximum 

intensity  is  half  the  detector  (CCD)  saturation.  The  radiation  is  incident  i 

1 

from  the  left.  The  inset  shows  the  photographic  data  with  a  vertical  i 

spatial  resolution  of  --  10  pm.  The  graph  illustrates  the  one-dimensional  ; 

axial  profile  taken  along  the  direction  ol  propagation  (z)  of  the  photographic  | 

data  (inset).  A  spacing  of  the  maxima  6  =  200  ±  20  pm  is  indicated. 

Fig.  2;  Calculations  for  N,  with  P  =  3  x  10‘‘  V/,  r,  =  3.5  pm,  Og  =  1  35  x 

cm'*,  and  Iq  =  8.6  x  10‘’  W/cm’;  (a)  Normalized  intensity  I(r,z)/Io. 

(b)  Normalized  electron  density  N(r,z)/No  for  N,  with  N<,  =  Og.  (c) 

Normalized  one-dimensional  axial  intensity  profiles  I(0,z)/I<|.  (A)  Full 

theory  for  data  in  panel  (a),  5  =  185  pm.  (B)  Calculation  with  charge- 
displacement  term  neglected.  (d)  Normalized  radial  intensity  profiles 
I(r.2j)/lo  corresponding  to  panel  (a).  Longitudinal  positions  Zj  =  172 
pm.  Zj  =  245  pm,  Zj  =  358  pm.  z,.  =  441  pm.  and  z*  =  559  pm  and  r^  = 

0.9  pm. 
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